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AN ALGORITHM TO GENERATE INPUT DATA FROM METEOROLOGICAL AND 


SPACE SHUTTLE OBSERVATIONS TO VALIDATE A CHj^-CO MODEL 

ABSTRACT 

This research addressed objective procedures to analyze data from meteoro- 
logical and Space Shuttle observations to validate a three-dimensional model 
which describes the transport and chemistry of carbon monoxide and methane in 
the troposphere. These efforts will ultimately provide the knowledge needed 
to understand the chemical dynamics of these species in the troposphere. There 
have been four aspects of this research: (1) Detailed evaluation of the varia- 
tional calculus procedure, with the equation of continuity as a strong con- 
straint, for adjustment of global tropospheric wind fields; (2) Reduction of 
the National Meteorological Center (NMC) Data Tapes for data input to the OSTA- 
1/MAPS Experiment: (3) Interpolation of the NMC Data for input to the CH|^-C0 
model; and (4) Temporal and spatial interpolation procedures of the CO measure- 
ments from the 0STA-1/MAPS Experiment to generate usable contours of the data. 
Each of these aspects is discussed in this report. 

4 
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INTRODUCTION 

Combined transport/chemistry models can be employed to analyze the circu- 
lation of pollutants from their sources to sinks. Depending on the time scale 
of the phenomena involved, these models can be broadly classified into urban, 
regional, or global scale. Sulfur dioxide and particulates are examples of 
pollutants that are important on urban scales, whereas the SO^-sulfate conver- 
sion is representative of a problem of regional importance. Finally, carbon 
monoxide is a species that has significant sources and sinks on a global scale. 

This report describes research to develop an algorithm to process 
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meteorologlcal and CO observational data for model simulations that would be 
compared with CO concentrations obtained from measurements on Space Shuttle 
missions. In the future, data from National Heteoroiogicai Center tapes and 
CO observations from Space Shuttle experiments will be used to generate input 
and comparison data for model simulations based on the procedures discussed in 
this report. 


OBJECTIVES OF RESEARCH 


The overall objectives of this research can be summarized as follows: 

a) Development of suitable spatial and temporal interpolation 
and extrapolation procedures to provide representative 
approximations of the discrete wind-field data which are 
necessary for the model simulations; and 

b) Development of representative three-dimensional CO data 
from two-dimensional observations, and interpolation of 
meteorological data for use in tropospheric species 
transport models. 

Most of the effort was expended on the first objective which, due to the size 
of the system, proved to be an unwieldy problem. Nevertheless, some sugges- 
tions relative to the second objective will be made. 


DISCUSSION OF RESEARCH 

The greatest emphasis during this research has dealt with the velocity 
field. This is most important since a mass-conservative wind field is essen- 
tial to minimization of systematic numerical errors during integration of 
turbulent transport equations. The approaches which have been evaluated in 
the past to guarantee conservative wind fields will be briefly discussed. This 
actually involves two separate problems: (a) Establishing a conservative wind 
field from the observations generally available at 0:00 and 12:00 GMT; and 
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(b) Interpolating between these observation periods (wind fields are required 
at approximately 1/2-*hour intervals for many models) to establish representa* 
tive "best wind fields" at intermediate times that are also conservative. The 
thrust of our research In these areas will be discussed separately. The pre- 
sentation has been divided into the following four sub-areas: (a) Spatial 
Adjustment of Data; (b) Temporal Interpolation of Meteorological Parameters; 

(c) Decoding the NMC Data Tape; and (d) Temporal and Spatial Interpolation of 
CO Data. 

Spatial Adjustment of Data 

A wind field that is free of errors and satisfies a particular set of 
general relations is a very important Input parameter for all kinds of meteoro- 
logical models, whether these models are for weather forecasting or for atmos- 
pheric transport studies. Unfortunately, such a wind field is not easily 
provided. The observed wind field data, as available mainly from rawinsonde 
measurements, contain observational errors and small scale effects, and as a 
result, do not satisfy general physical relations, say mass conservation or 
other diagnostic relations applicable to a particular atmospheric model. Such 
a wind field when input to model computations can lead to highly inaccurate 
results with initial errors multiplied many fold. In order to overcome this, 
it becomes an important and essential task to treat and adjust the observed 
wind field data in such a way that the observed data are minimally adjusted. 

But at the same time, the adjusted values should satisfy the physical relations 
applicable to a particular atmospheric model. For example, a mass consistent 
wind field is a must in the transport and diffusion atmospheric model. Thus, 
for such models the observed wind fields are adjusted so as to satisfy the 
mass conservation condition, i.e., the continuity equation. 

Several methods of varying degrees of complexity have been investigated 


-i, 


to accomplish this mass conservation. The first and simplest Is establishing 
a vertical velocity component by nunerically differencing the equation of con- 
t i nu i ty* 


3u 3v 
3x 3y 



( 1 ). 


Thus, from Equation (1) the vertical velocity is immediately calculated taking 
w > 0 at z > 0 to initiate the calculation. This procedure has been used with 
completed initial simulations (Peters and Jouvanis, 1979 and Carmichael and 
Peters, 1981). 

This procedure is very simple to use. However, due to errors in the re- 
ported observed wind field, the calculated vertical velocity can be more repre- 
sentative of "noisy data" than of the actual vertical motions. As a result, 
representativeness of the wind field may be sacrificed and not attained by 
this simple method. The procedures discussed next recognize that there are 
observational errors and make allowar.ee for such. 

A technique, employed in the past (Endlich, 1967 and Goodlin, 1976) uses a 
point iterative method to adjust the observed wind field data so as to satisfy 
the desired physical constraints. Goodin et al. (1980) recently used a hybrid 
of the point iterative procedure with numerical satisfaction of the continuity 
equation. 

A seemingly more general and powerful technique is that based on varia- 
tional calculus principles. This technique, which was first proposed by 
Sasaki (1958), allows the observed wind field data to be changed by a minimum 
amount in a least squares sense, while at the same time, adjusted values satis- 
fy the imposed physical conditions within numerical approximation errors. 

Subsidiary conditio ns that are to be satisfied exactly are termed strong 

^ All of the development presented in this section will assume constant density 
for simplicity. However, analysis using a varying density is an obvious and 
straightforward extension and is being included in the actual analysis. 
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constraints, while conditions being imposed approximately are weak constraints. 
This technique has been extended to a number of situations with various forms 
of constraints. Washington and Duquet (1963) used the technique with the geo- 
strophic approximation while Sasaki (19S8) and Stephens (1970) applied the 
balance equation as the dynamic constraint. Achtemeier (1975) extended the 
technique with a number of constraints such as the momentum, energy, continuity, 
and hydrostatic equations. Later work seems to lead to a more complex network 
of non-iinear equations that must be linearized before solving them. Recently, 
Dickerson (1978) and Sherman (1978) used a variational technique to adjust the 
observed wind field in the San Francisco Bay area under the strong constraint 
of the continuity equation. 

Generally, the mathematical formulations of the variationai technique may 
be described as follows. Let u**, v**, and w° be the observed velocity components 
of a wind field, and u, v, and w be the corresponding modified adjusted values. 
The difference between the observed and adjusted field may be expressed in a 
least squares sense as 

a,^u - u**)^ + 02^(v - v»)2 + a^^w - w«)^ (2). 

« 

« 

where 02 » and are weighting functions defined as the reciprocal of twice 
the variance of the errors of observations for the respective observed compo- 
nents. The simplification of » 02 has been used in the past on the basis 
that the errors in the horizontal components should be comparable. In order to 
obtain the modified values, u, v, and w, it is required that the adjustment 
over space and time, subject to the satisfaction of the constraint equations, 
be a minimum. An adjustment functional, therefore, may be expressed as 

I » / / [a^^(u - u®)^ + 02 ^(v - v®)^ + Oj^(w - w®) + yG^ + XG 2 ] dtdl' (3), 
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where is a weak constraint which is introduced into the functional through 
the constant y The strong constraint G 2 is introduced through the Lagrangian 
muttipiier X. Equation (3) defines a variational problem in which a minimum 
vaiue of the integrand under the given constraints G^ and G^ is being sought; 
i.e, the first variation on I must vanish, 51 « 0. As a result, the associated 
Euler-Lagrange equations along with the so*cai1ed natural boundary conditions 
being inherited from the vu'*iationai formulation must be satisfied. For ex- 
ample, a functional such as given in Equation ( 3 ) with the continuity equation 
as the only strong constraint and no weak constraint yields the following 
associated Euler-Lagrange equations along with natural boundary conditions 
(Sasaki, 1958). 
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The natural boundary conditions are 


n X5u 

X 


n X5v 

y 


n X5w “ 0 
z 


( 8 ), 


where 5u, 5v, and 5w are the first variation of velocity components, and n^, 
n^, and n^ are the outward positive unit normals in the subscripted directions. 
From Equations (h) through (7) • an expression for X may be obtained as 


* Note that » 


02 has been assumed. 
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3^X ^ ^ ,®1 . 3^x _ , 2 /3u* . 3v* . 3w% 


(9). 


Equation (9) along with boundary conditions as given in Equation (8) can now be 
solved fo^ X(x,y,z), which in turn is used in Equations (4), (5) > and (6) to 
estimate the adjusted values of the velocity components. 

The adjustment procedures under the continuity condition as a strong con- 
straint, as described above, lead to working equations which are in principle 
easy to apply. These procedures with Equations (3) through (8) have been 
evaluated for small regions, such as the San Francisco Bay area, by Dickerson 
( 1978 ) and Sherman (1978). For the global wind field adjustment problem, the 
validity of these procedures has yet to be completely evaluated. On the global 
scale, it may be required that the wind field not only satisfy the continuity 
relation but also retain other observed kinematic properties. One of the im- 
portant kinematic properties of the wind is Its vorticity. In such cases 
another constraint defining vorticity of the wind may be introduced into Equa- 
tion ( 3 ) through a second Lagrange multiplier X^. Consider that it is required 
to retain only the vertical component of vorticity, c’* estimated from the 
observed velocity components. Then the second constraint may be expressed as 

( 10 ). 
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The functional as defined by Equation (3) with the continuity equation and 
vorticity. Equation (10), as strong constraints leads to the following Euler- 
Lagrange equations 
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with natural boundary conditions 


n^X^6u - ® 

"x^2®'' + nyX^6v ■ 0 

n^X^Sw ■ 0 
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X^ and X 2 are the Lagrange multipliers for the continuity and vorticity con** 
straints, respectively, and the following equations for X^ and X 2 may be ob- 
tained from Equations (11) through (15). 


2 2 2 2 
3^X^ 3^X^ 3^X^ 

5” + »- + (■“j) 2* 
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3^X, 


( 20 ). 


3x 3y 

Equations (19) and (20) along with the boundary conditions (Equations l6 - I8) 
can be solved for X^ and X2 which can In turn be substituted In Equations (11) 
to (13) to estimate the adjusted values of the velocity components. 

This variational formulation Is applicable to those cases In which only 
the observed vertical component of the vorticity Is to be retained. Similarly, 
cases constraining all three components of vorticity may be formulated. Such 
a formulation not only leads to more complicated equations which are more diffi* 
cult to solve, but it also involves vorticity terms which require observed 
values of the vertical velocity. In practive, the vertical velocities are 
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rarely measured at meteorological stations* Therefore* it leads to considerable 
uncertainty and difficulty In calculating the observed x- and y-direction vor- 
ticities. 

Any or ail of the fundamental equations describing the state of the atmo* 
sphere could be used as constraints in a variational sense. However, only 
Achtemeier (1972) has attempted to treat the observed data so they satisfy 
f the complete set of equations comprising the mathematical model for the atmo- 

f sphere, and he did that based on the variational method. This Is probably the 

best approach in treating the limited experimental data obtained every twelve 
hours, and if applied, it would yield sets of “data" which, in terms of the 
mathematical model used, consistently describe the state of the atmosphere at 
discrete points in time. 

I 

I To summarize, the four procedures represented by (a) Equation (l) alone; 

(b) Equations (k) through (9); (c) Equations (ll) through (20); and (d) the 
completely constrained formulation represent a comprehensive analysis of global 
wind field adjustment procedures. Only Procedure (b) has been studied in this 
^ research. The computer program to solve that set of equations is presented in 

the Appendix. ‘.The program is written in SL/1 for use on the NASA LaRC STAR 

I 

I computer. 

I Temporal Interpolation of Meteorological Parameters 

The preceding analyses deal with the meteorological configuration based 
on observations at a specific time. Generally, observations are only available 
every 12 hours while data describing the state of the atmosphere may be re- 
quired at the intermediate time periods. As mentioned previously, transport 
modeling efforts frequently require such information at 30 minute intervals. 

In this section, a procedure to achieve this will be described. It is an ex- 
tension of the mass-consistent wind field methodology presented in the previous 
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sectton. (These temporal interpolation studies may also be applicable to the 
data reduction procedure for the MAPS instrument which requires the temperature 
(to within 2*K) and water vapor profiles to determine the CO concentration. 
Again, since the data are only available at 12-hour intervals, the projection 
of these meteorological variables at the intermediate times at locations under 
the Space Shuttle track are required. The procedure that is suggested could 
also be applied to those requirements.) 

The variational formulation for spatial adjustment Is applicable to the 
observed data available at one time level. However, a similar formulation may 
be extended to any time level, and interpolation of data between two time 
levels may be obtained in such a way that the interpolated data satisfy the 
given constraint ($). Consider that observed values of wind field data are given 


for time levels t^ and It is required to estimate interpolated values 

at time t where t < t < t . If linear interpolation between time levels is 
n n+i 

assumed, then the functional in this case may be written as 


” ' 9 9 7 7 7 7 

I* • / / (a. (U - U*)^ + o/(V - V*) + CL (W - W*)‘ 

Vt ’ 




where 

u - u u . , - u 

U-r ^ and ^ (22), 

‘ ‘n ‘n+1 *n 

with similar definitions for V and W. It can be shown that a simple transfor- 
mation on the velocity components can reduce Equation (21) to a functional si- 
milar to that of Equation (3) from which Euler-Lagrange equations and boundary 
conditions are obtained which adjust the interpolated wind field so that the 
continuity equation is automatically satisfied. Thus, the wind field proceeds 
smoothly from the configuration at t^ to that at t^^^. Equation (21) and the 
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associated conditions apply to linear temporal interpolation, but non-linear 
interpolation schemes are readily handled. Finally, the additional constraint 
of vorticity can also be incorporated. 

Decoding the MC Data Tape 

The National Meteorological Center Data tape includes considerable data 
that are not necessary for interpretation of the OSTA-I/MAPS experimental data. 
In addition, separate tapes are used for the Soutern and Northerm Hemispheres 
and these data must be merged. Finally, the data must be decoded for use by 
the MAPS Science Team. A substantial effort was expended in this phase of the 
project. 

Weather data tapes produced by the National Meteorological Center have 
been coded in a binary format not easily utilized by most computer programning 
languages. For this reason, a preliminary procedure must be performed on the 
data, to transform it into a standard EBCDIC format. (See the Appendix for a 
description of the data on the NMC tape.) DECODE is a PL/I program which con- 
verts the NMC binary data to their EBCDIC equivalents. The tape format Is also 
modified so that only certain desired fields are processed and retained. A 

e 

program listing of DECODE is also contained in the Appendix and has been de- 
signed for use on the IBM 370 at the University of Kentucky. 

DECODE works by manipulating PL/1 data structures. Data are read into a 
structure of bit strings. Then, using a subroutine, the binary data are con- 
verted, first to fixed binary numbers, then to EBCDIC characters for output to 
tape. Conversions are done based on the equations given In Office Notes 8b and 
l8b from the National Meteorological Center (see the Appendix for copies of 
these). A selection is made so that only the desired data are processed. All 
included data are of type 8. Data are Included for the following conditions: 

1. Q ■ 1, height with respect to mean sea level geopotentlai 
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(at 12 levels), 

2. Q ■ l6, atmospheric temperature in degrees Kelvin (12 
levels), 

3« Q ■ 48, u'component of wind (12 levels), 

4. Q ■ 49, v'component of wind (12 levels), and 

5« Q " 88, relative humidity (6 levels). 

A sort is then done to sequence the data in the proper order for merging 
into longer records. The input records contain all the data for 5,36$ points 
on a southern or northern hemispheric grid of l45 x 37* Two Input records are 
required to store all 10,730 o* the global values for each time period of a 
two time-period day, with 15 days making up a complete tape. Alternate data 
values have been omitted, leaving 2 grids of 73 x 19* which corresponds to a 
5” X 5* gridding. The output records contain data for 2,774 global points in 
a west to east, south pole to north pole order, with the sort field values in- 
cluded for Identification purposes. In this manner, selected data from a full 
time period can be contained in a single record* Data fields (and their for- 
mats) included In the output records are 


OATA-TYPE 

« 

F(I4) 

type of surface; 

Cl 

F(I4) 

value of surface levels; 

Y 

F(l4) 

year; 

N 

F(l4) 

month; 

0 

F(i4) 

day; 

1 

F(14) 

time period; 

DATA (1337) 

F(9) 

data values from the southern grid; and 

0AT2(1387) 

F(9) 

data values from the northern grid. 

The records are sorted by year, 

month, day, time period, type of surface, and 


value of surface levels in descending order. The type of data is sorted in the 
following order: 48, 49, 1, 16. 88. 
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DECOOE has been designed for maximum ease of utility by the end-user. 

Only one modification must be made prior to submission of the job with a new 
tape. JCL cards describing the input and output datasets must be prepared for 
each run of the program. Input (the NMC tape) is defined as SYSUT1 , illus- 
trated by the following example: 

//SYSUTl 00 UNIT - 3l»00-5, VOL = SER » tape-number, OISP « (OLD, PASS), 

// OSN « name.of. dataset. on. tape, LABEL « (f Me-number,SL) 

Output records may be sent to either disk or tape. However, the JCL for 
tape is simpler to prepare since the output records are longer than a track 
and therefore would require V6S records on disk. Be sure to specify a large 
enough SPACE parameter when using disk. A sample JCL definition for output 
on tape follows. 

“ //GO. TAPE DO UNIT - 3l»00-5, VOL « SER - tape-number, DiSP = (,PASS), 

// OCB = (RECFM - F,LRECL - 25050, BLKS IZE » 25050), 

// DSN - name. of. dataset. on. tape, LABEL « (f i 1e-number,SL) 

Based on partial tape runs done on the University of Kentucky Computing 
Center's IBM Model It computer running OS/MVT, a projected cost esti- 

mate may be made with some degree of accuracy. Since a load module will be 
used for the production runs, compilation costs for this projection can be 
ignored. Processing only the 2 time periods of the first day took 55 CPU- 
seconds of which 3 seconds can be attributed to compile time. Using the re- 
maining 52 seconds of execution time, the approximate cost is $7*52 for de- 
coding the first day. This is based on a $525 per hour charge (the discounted 
rate for deferred priority). Multiplying by the 15 days on each tape makes 
each tape cost about $113 to decode. It should be understood that this is only 
an approximation, and the actual cost will depend on various circumstances at 




the time the Job is run. Other computer Installations use different charging 
algorithms, so costs will vary even more if run on a different machine. 

Interpolation of the NMC Data for 0STA~1/MAPS 

The NMC data as decoded provide horizontal winds, temperature, and dew 
point, all at the pressure levels of 1000, 850, 700, 500, IjOO, 300 , 250 , 200, 150, 
100, 70, and 50 mb. An interpolation scheme was developed to represent these 
data in an alternate form. This scheme provides interpolated values for a 
5** Longitude x 5” Latitude grid and a non-dimensional vertical coordinate 
which has p ■ 0 at the surface and p • 1 at the tropopause. The data are 
provided at p - 0, 0.0333, 0.0666, 0.10, 0.20, 0.40, 0.60, 0.80, and 1.0* This 
gridding provides a finer structure within the planetary boundary layer. 

With the data unpacked from the NMC tapes consistent with the requirements, 
the data can then be interpolated for a geometric gridding. The mixing layer 
depth will also be calculated at each grid point and will use the 1-dimensional 
planetary boundary layer model of Yamada and Mel lor (1975) • Those researchers 
have provided a version of their program which is being modified for the pre- 
sent purposes.^ 

« 

Tempot\zl and Spatial Intcr^x>lation of 00 Data 

The CO measurements from the OSTA-l/MAPS experiment on board the second 
Space Shuttle flight will correspond to mid-troposphere concentrations below 
the orbit path. These data will be from 38°S Latitude to 38°N Latitude and 
represent variations in both space and time. Thus, it is necessary to Inter- 
polate the concentration-time-position data at the desired times. The following 
treatment would accomplish this. 

1) Transform the individual 1-second signals to correspond to 
non-overlapping, discrete orbiter track lengths of 10.37 
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2 ) 


3) 


s. 


nautical milas. 

The information from Item (1) represents a concentration 
at some time t|^, and some position Tj. The position is, 
of course, expressible in terms of latitude and lonqitude. 
In order to construct an entire concentration field at 
some time t^, interpolation in both space and time is 
required using as much of the data as seems reasonable. 
Define the following. 

c? = Concentration at time t and position r, 

l( 

fjj = Weighting factor for temporal and spatial 
interpolation scheme 

r. , = Distance between positions r, and r, 
fj -i -j 


^kn ~ ‘‘n 

q H Total number of locations which are considered 
to influence interpolated value 
s = Total number of time intervals which are con- 
sidered to influence the interpolated value. 

For the present case, s * q. 

Then, the interpolated concentration is simply 

c, - E E f, , c, 

‘ k-1 J-1 J 

If one considers ail of the data, this could be simplified to 


c"- ? fJ c-^ 

k) Now the crucial question is, What should be used for f!^j7 
There are a few papers on this topic for aerometric data, 
and one scheme used is an inverse square dependence for 


(23) 


( 24 ) 
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spatial Interpolation. Nobody has apparently reported both 
spatial and temporal interpolation. Nevertheless, one plan 
to start with is the inverse square dependence for both space 
and time; i.e.. 


r’2 t*2 

.k _ j 'kn 
U “ q .2 s 
^ •':T ^ K 

j=i k=i 


(25) 


In this manner, is of course normalized. Other weighting 
schemes could be tried. 

5) The procedure in Items (l) - (4) provides interpolated values 
over a 10.37 nautical mile region. Since that is too much 
data to handle reasonably, these values could be averaged over 
1" X 1®, 2.5® X 2.5°, or 5° x 5° longitude- latitude intervals. 

Then these data could be presented in digitized form as well 
as contours. Contours will be most helpful to see large 
scale plumes. 

This procedure is easily coded, and the weighting function could be varied 
to evaluate its effect. Furthermore, several times could be done at once, and 

q -2 

one may be able to take advantage of the data format, and calculate E r , 


just once. 


j»1 
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CONCLUDING REMARKS 

The procedures discussed in this final report will be used to assist in 
the analysis and Interpretation of the CO observations obtained from the OSTA- 
1/MAPS Experiment which will be aboard the second Space Shuttle flight 
scheduled for November, 1981. Specifically, DECODE will be used to provide 
meteorological data from the NMC data tapes In a format useable by the MAPS 
Science Team. This will include u, v, T, and water vapor concentration data 


for the flight period in either a pressure coordinate or a dimensionless 
terrain/tropopause following coordinate. In addition, an interpolation pro- 
cedure for the 0STA-1/MAPS CO observation data has been proposed that includes 
a spatial and temporal weighting function. 

A considerable effort has been expended on objective analysis procedures 
to establish global mass-conservative wind fields from the NMC data. The 
variational calculus procedure with the equation of continuity as a strong 
constraint has been coded and Is in the final stages of evaluation. Novel ad- 
justment procedures which incorporate vorticity preservation as an additional 
strong constraint and temporal interpolation during adjustment have also been 
presented. However, these procedures have not been evaluated. Research is 
continuing in this overall area of wind field adjustment techniques (Kitada 
et al., 1981). 
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FILE PAMLAN* 


-/g<at , p;^tniam 

./cxipvphf ( paml 9m4 

A Monui. F-WIWnAn. lS( UNITl=CriPMTRX»UNIT?=f;ini. MTRX, UNIT3=- 


-AOJMTRX U-INI T/>=OUTPUT-^ l.lN-1 TJ5= I NPUT-)-BFG I N ~t*IA IN ; 


-iXRFFi- 


-jU-ULF-^WINDF rELD-MODULE-^=_PU/RAR I fi-r^FP-TFMPFR -^-,-1 981- 
->Dt SOURCF-s-P AMI . AM * / 



BATA_s_l.AMDATA_ 
JCI— 55-FAMCTL . 


-*JL 


./.*..THIS-PROGRAM_AD.JUS-TS_OBSFRVFa.UTNn-F:iB.n-TO-MAKF- 
_I.I_:MASS=CONSERVATXVE_H.I NH -HI Fl.D . AH JUST I NG -T-EOHN I OUE- 
.IS-BASFD . ON- VAR T T I ONAL -FORMUI .ATI ON . --:5 1 MULT ANFOUS 


-EQUATIONS -FOR I AMOA -TERMS - ARF-SOLVFD-USING -SUGCFSST VF 

-OVER- RELAXATION iSr-O-rRi-I-TERATION METHOn.-F-W FS- NAMETl -AS 


XOFMTRX.vF:OI.MTRX.,.ANin.AnJMTRX — HAVF-NESSASARV- 


.-EST.IMATFD-JC:GFFF-IC;IENTS, - THESE-PTLFS- WU.L-BF -RFAD-TN-THIS 

-PROGRAM. *J 

./* ^THF- EOU.OWr NO-SECTION NFFOS- TO PF- CHANGFP-FOR -niFFFRFNI 

-GRID SYSTEM. XTRAl-J-Q-XTRAT -AND-ISMt TO.ISM7- ARE- 


PAGING..PARAMETERS.-.-ALl_XS..WE.I.GHXING.-FACTOR-POR_U-VELOC;r.XJiL 

... COMPONENT. , 


ORIGINAL PAGE IS 
or POOR QUALm* 


* Anyone interested in this program should contact L. K. Peters at the 
University of Kentucky. 


r 


21 


-COMMFNTINO 
h--/* F1LF-<S- IS 


/* MOniFICATION HISTORY 
FFRRI.IARY— 18. 1^81 




RFOARDINO FS 
NOW OUTPUT 


IMPLICIT I/O 


»/ 

*/- 

*/- 

•*/ 


-/*- 


-*/- 


-/♦— BFGININO OF THF- SPOT I ON 

-l.-TTFRAI.LY-''5,0" DYY» 

-LTTFRAI. L Y -•'■7?" -NUM. FWl 

-I.ITFRAI.I.-Y -*'37'' - NUM.NSI 

-l-ITFRALL Y- ”9''NUM ,Al.T« 

-IrlTFRALI Y-"10''NUM_Al.TPH -/«— NUM. 

- L I TFRAI J.Y " FW .. NS .ALT < 

-L I THR Al .1. Y-“ 3357/. " - X TR A 1 1 

-t-TTFRAl.l Y— " 15992" - XTRA2I 

- L-I TFRAI, I ..Y J 7584 "-XTRASt 

-I.TTFRALI.Y ••24512" XTRA4? 

- LITFRALIrY- "51 1 52" -XTRA5« 

-4 T TFRAI l.Y " 1 1 192" - XTRA7; 

„i_l -TFRAL LY 7/.S •• - 1 SM 1 ? 


-*/- 


,ALTP1“ -NUM..AI..T+1 */ 


l-rl TFRAI. IrY-^-^ 640 ••- I SM2 1 

1 1TFRALLY-" 128^^ ISM3? 

LITFRAI.L-Y-^256^^ - ISM4; - 

L I TFR Al. I . Y •* 256 •• - 1 .3M5 S 

L T-TFRAI. I . Y - •• 256^^- - T SM7 ? 

FF-4 . T TFRAL LY " 0 . 0 1 " Al. 11 


-I FF»/ 


— L-I-TFRAL.LY i!3. 141592653589793"PI 1 

_/«- NIIM FU=NUMBFR OF ORIOS IN FAST-WFST OIRFCTION. 

-NIJM..NS=NUMBFR.OF .i:iRTDS--IN NORTH-SOUTH DIRFCTION. 

—NUM. ALT=NI.IMBFR- OF- GRIDS IN ALT I TUDF- DIRFCTION. - — 

- FW .NS . AI.T*TOTAL NI.IMRFR OF GRIDS. 

_JCILFRANCF=T01..RANCF IN ITERATIONS. 

—RLXFACsRFI.AXATTON FACTOR-IN SUCCFSIUF- OVFR RFLAXTION 
— MFTHOn. - ♦/ 

FND OF-THF-SFCTION */ ^ 

_F!FAl^VFi:.T0R..n51 ARRAY (FW. NS. .ALT > -TXl — 

_J-;OMMON CTlfjC.-TnilMPJ >; 

— RFAL-UFCTOR- rXTRAI 1 TDUMP1; 

RFAL-VFCTOR r 2664 1 MY ALT 1 

— RFAl .-VFOTOR 151 ARRAY ( FW -NS .ALT ) - AD..IX « - 

.-COMMON CA ( AD JX , TOI IMP3 ) X 

—RFAl VFC.TOR CXTRA71 TDUMPS; 

—J1FAI. VFC.TOR IFW .NS.-AI TI ARRAY! 13) SSX{ - 

—COMMON . CSS ( SSX . TDUMP4 > | 

_ . RFAL-VFCTOR -t XTRA21 . TDUMPAt 

_ Rt£AI. VECTOR . C 151 ARRAY (FW .NS-ALT) -TI.AMDAl 

_COMMON PI t-TLAMDA, TDUMRl 1 ) 1 

_RFAL. VECTOR 1 XTRAl 1 TDUMP-IH 

— RFAl VFCTCR IFW.NS. ALTl S.DIV: 


-COMMON- B2. (S.DIV, TPUMPS ) ; 

RFAl VFCTOR rYTR^-^T TnuMOc:. 


■ fc M . . ** i..t ^ t i-MM n ....... 

PEAI. VKCTOR C51 AOJI 

RE Al. • VECTOR -fNI.IM. FW3 ARRAYINUM. NS. NUM_Al.TV‘Ar AMOA, Rl. AMHAt 

COMMON B3 ( Rl . AMDA , ALAMDA. TOUMPiS ) I 

RFAL VFC.TOR rXTRAOl TDHMP6I 

REAL -VECTOR CNUM..EW1 ARRAY(NI.IM«NS, NUM..AI.TPI ) 

II, V , w , n » AOJI ISTU. AD JIJSTV , 

common P4(U. V.W.n.TDUMP7)» 

REAI;- VECTOR -C XTRA41 -T0UMP7I ^ 

COMMON PS ( AO. n iSTU. An..H.ISTV, ADJI.ISTM. TDOMPS ) » 

REAL- VECTOR [XTRA51 TDUMPOl 

REAL- VECTOR CEW_NS. ALT] ARRAY ( 1 3 » - VEI » 

COMMON P<S ( VEI TOUMP9 > I 

RFAI.- VECTOR -r XTRA2 ] TnUMP9» 

RFAl.--VECTOR-CNI.IM .AI.TPM It 

INITI Al - ( 7«\0.0. 0, 0333334,0. 0333334. 0.C>333334 

,0. 1 .0.^.0.2,0.2.0.C>.0.7\); 

REAL- Al-3 . T . ERROR , 777. THETA, TH, CON, MAXD 1 V , M I NO I V,— TOl FRANCE . RLXFACt •< 

^ 1 NTEOFR-I MTNOTV.TMAXOIVt 

^TNTEOFR -NOT. CONV, r , J. K, NL, 1 1,17, 13,-14,^11 , J?. K7, 1 

/*». EEK3, I TER AT T ON, K t , I W , E I , E7, E3, F4, E5, NIJM. . T TER , LEE*/-- 

K3, ITER AT I ON , K1 , 1 W , E 1 , F? . E3, F4. F.5, NI.IM. I TER, MAX .. ITER, 

, EA . E7 , ES , E9 , E 1 0 . MMMM , J. J , NS. . I NCR ,--AL-T. . t NCR , - EW_OUT ? 

REAL-VECTOR -r3331. Wl AMOA, WL AMOAll 

REAL VECTOR --C 91 KL AMOA 1 

RFAI.-MAXI. AMnA. MINI AHOA.-ALU 

/«•«•♦«■****«•**«■«•«•#«•*«• ■»•**«•*♦**«•********♦«•**«•**«•*♦«•***«•«•**«•«•*••*«•*«•*♦«•««•**/ 

/.***«******»****»********************************»***«***************/ 

y*.**.«**«-*<HHHHHHHt-**fr*fr*«.**tt**»*»***<HH!-****«-********fr«-»-*****«'**«-«****»/ 

/ft***#*##****#*##**###***#****************************#************#*/ 

_/*«-«-«.***.«.<HH(»«.lHm.*«.«.«-****«-***«.****««******IHItt*<HH(**«-«#**«-*-«-**«-******«-/ 

— 4PA0E: 

PROCFni.lRF- - MA T N ; ; 

/.#LEE*/- 

WR I TE-LIN I-T f.A ) ♦» ' 1 BEG I N • MA I N PROCEDURE t 

REAn-UNTT (V.) 


( RI.XFAC.--TOI ERANCF-.-AL 1-,-- MAX^J-TFR, - N.3.,1 NCR.,- ALT-.:: VNCR,- 

EUI..OUT) 

43 ( E9 . 0 - /.) , -3 C 1 9-./.U—1 94 *, 

— /*l.EE*/ 

=-=-_PLR - JANUARY. 23 -1981. 

lE-EW OUT- >-NIIM..EW -THEN 

WR ITE - UN IT ( A >_ I FW. OUT-) 4.SX '’U.SER- WANTS T5tt f 

EW_OUT:«imuM_EW< 

END It 

i OPFNMAP( <F1 '.TLAMOAI 1 )m, I3M1 .FI .-n.AROFM t 

OPFNMAPt '■F2'-,3Cn. ISM3,F2, "LAROE-'-lt 

OPENMAP.l ^F.3-'‘-,Pl.AMriA( 1 , 1 ) C l ] . ISM3.E3, ^LAROF-^^) t 

OPENMAP ( 'F4 1 . U ( 1 , 1 ) C n , ISM4, F4 , "I.ARGE^-) t 

DPENMAP.LlF5L..i^nJLISTU( 1 .-D C 1.3 ..ISMS.ES.-^! ARGF U-t 

OPFNMAPt ''FiS-' . VF! ( 1 1 C 11. TSM7.FA. 'i ARGP-’ 1 1 


i: 

I 


OP’FNMAP('CnFMTRX".TX(l )n 1. ISMl ,F8. 'LARCiE" > « 
-OPFNMAP (-•' Sni.MTRX " , SSX <-l > T U » T SM 2 ♦ F9 , 'LARGE f 


. OPFNMAP ( ' An..lMTR X ' t An. IX ( 1 ) T IJ . T ft W7 , F 1 A , ' I . AROF ' > t * 


-/♦LFF. - WRTTF I lNIT<f.) (Ft »F^» F3- F4. E5, E6. F8. E9 - 
,FIOH»'OPFIM'-, 9(3XI?MMI FF»/ 


-WRTTF -HNIT( 6 ) (Fl.F.7>.F3.F4,F5iF/S.FSvF9 — 


-,Ein)H' OPFN'.-^9(3XI^i:H*J 
NI.IM^TTFR !•- 0« 


-/ *l F.F — WR T TF- MN 1 T ( h V 


•Rl.XFAATOR 


F6.34(LER*/ 


-WRTTF UNIT(<S) 



( nVY> Nl IM_ FW , NI.IM..NS, NL1M_ Al-T , Al. 1 , TOLFRANCF i RI'XFAC ) ♦» ' PW 

FR.7t 

— -^NIIM. FW - I3v'NUM..NS ■ " , 13. 'Nl IM_ALT 'T I3r 
' Al. t — F7i 3. ' TOL FRANCE , F8. 4, 


( nvv, Nl IN . FW . Nl IN-NS . NUM.. ALT . AL 1 . TOL FRANCE . Rl. XFAC , - 

MAX.TTER)#'- nVY » '-.F5.7. 

^--NHM. FW . 13, '-NIJM. NS ■ r3,'^-NIIM.. Al T “ -^ , T3, 

.AL1-- '-,F7.3,<— TOLERANCE ■ ',F8.4, 


-^Bl XFACTfiR-= ' , FA-3, 

4:.flN « » PTV-1 30. Ol 

_K 2 »• NIIM.Al T +1 t 


MAX. TTFR 13*1 1 - 


/». INTTrALIZ INC, ALAMDA AND BLAMDA- ARRAYS-OT VINO — 

FTR3T 0LIF3;5 VALI.IE3 */ 

.7* -I AMBDA VAl IIFS ARF FOUND ON PAOF FILE F3 */ 


-7.* WHFN -THF FOl.I .OWTNO IMPLICIT- I /O IS 

— T-S- COMMFNTFn OUT, 01. P- VAl UFS-FROM SAVED 
-/* -FILE- F3 ARE NOT OVER-WRITTEN 






*/- 

*/. 

r-.— */— 


+ DYY* 

TO NUN.. NS no.. 


THETA»»90,0 

FOR .1 1 

theta..! a THETA-DYY! /*NOT JMEEDED 


PL.R-FEB 17,31 .*/ 

A* I.IR-THF.-.INI.T.IAL . VALUES- PY-1 . OE+OA — */ 

FOR K - Ja 1_T0-NUM .Al T DO 


_Bl AMDA(^l,KXn :NUM_EW1 


. ALAMDA ( . I , K >*C1 1 NUM F.W 1 


-ia 3.0F+0A* 

I ■ -7. OE+OA! 


-ENDF! 

ENDFl 

.EOR_<l_:a_l._TO. 


NUM. NS no 

FOR -K_x« -1 -TiX-NUM .ALT DO _ 

FOR -I.i»_l BY 2 TO NUM_FW DO 

m.AMDAl J.K) Cn -X»-l.OE+OAt 

ALAMDA (J,K) m «*_1 .OF+OAt 

ENDFt 

£NDF! 




-ENDF! 


. END OF_I /O.^— -INSERT COMMENT - HERE . -•:=■--*/- 

---.FLOAT APRTI 7TH — */ 

BLAMDA(5.f.> C3A3 -» ■ 1 .5F+0At 

_ ALAMDA(5,5).C3A1 t 1 ,5E+0A* - *✓- 

FOR .1 :a 1 TO NUM W«; r.O 


I 


Wl.AMnAC..n I- BI.AMDA(J,!5)C3<SH 

FWnFl 

— •'♦I FF*/WATTF -IINTT<<>H»*'1^WI AMHA ARFi 'Ft /«! FF*/ 

.4RTTF IJNIT<(S) (WLAMDACtlNUM .NSl )F5F14. 7F« 

FOR- K- »■- I- TO -NUM. Al T-PO 

Kl . AMOA r K-1 - I «- Rl . AMOA < 5 . K ) C t 

-V*IFF- WRITF IINTT(6KKI AMPAn t Nl IH„Al T 1 Hl'Kl.AMDA - ? 5 F 1 4 . 7Fl 

WRTTF-»tNTT(6HBI.AMOA(f5.5)C.'^Oi403>#'l.AMriA--— ‘T!5Ft4.7F! 

WRTTF UNIT(/snPlAMnA(5,5)ri 3.W-AMIiA(5»5)rNl.lM. FW3 ) 

#"LANDA ■ •^,2F14,7Ftl.FF*/ 

— /HFF*/WRrTF IINIT(6> ( Kl AMDAr 1 » NMM.ALT3 >«♦' KI.AMDA ■ ' ,?iF14.-7Ft 

WRTTF I.INTT(«.)<RI.AMnA(5.5>C30i401>#'-I.AMDA ■ 'iSFlA.TFt 

WR I TF--I IN T T ( M ( Bl AMOA ( ?. 5 ) f I 3 . PI AMPA ( 5 » 5 1 f NUM..FW 3 ) 

-l.AMDA -■ -i.,7F\4,7F» 

FOR -I - l« I TO t.3 no - 

VFL< T ) n tFW. NS-Al.n t»- O.Ot 

FNPFt 

/»_ I NT T T Al r 7 I NO W 1 Nr F T FI P -4 U V, W ANP ■ PFNSI TV- P- ARRAYS «7 

THFTA *■ -90.0-DYVt 

for-vI-i« l-Tn.NllM. ^'^ -no 

THFTA I ■ THFTA+riYYt : 

TH I ■ -THFTA* nONt — 

FOR K-i« I TO K2 DO 

IMJ.KItl SNIIM. FWl 

-SJ ( ..I , K ) r t : NI.IM. FUn - -3 . 0 1 

WC.UKIC 1 INIIM. FW3-J- O.Ot^ 

nt.I.KK liNUM FW3 1 --T.O 1 

FNPFt /*! OOP ON K — LFF»/ 

F.NDFt /*I.OOP ON -.J I.FF*/ — — 

T.AI l.--PIVRONCRt 

CA|.l..--SORt 

/* — AI.IO -77TH.4 981 «/ 

A*_pnTNT OI.IT-1 0 — l.AROFST ANP -♦/ 

4*~f;MAl.l FST-PTV8 PFR I AYFR. . . OFl FO RFTiIRNS *7 

/*_f.-of^noi.iccf!?.sfmi. comparisons I - */ 

4ri3R_Kl- .1 TO-NIIM .Al T-PO 

-NL I ■Ot — — ' ■ — 

WR I-TF-HN I T I ^. ) (KltttHl . ///, 40X , 'Al TITIIPF ♦ 13.7/ /#t 

FOR ...It- I -JO NIIM..NS DO 

FOR-I.t«t -TO NMM.FM PO 

— 

MVAl.TCNL3t»PIVr (K-1 )*77'+(J-rl )*A4R-t-I 


FNPFt 

ENPFt- 

MY ALT t.l t 2464 3 t « APS ( MY ALT f 1 1 24A4 lit.. 

FOR JJ t-.UT0-20 DO-- 

MAXPIV: =MAX(MYALTU 

MTNDTV;»MTN(HVA!.7)t . 


IMAXniVi-?^Fl FO<MYAI.T,MAXnTV)-H « 

IMlNniVl-SFLFA(MVAkT»MlNntV)-fn 

i,i»i-iMAxnTv.niv.7^-m 

ut»ihAxnTv-( (.11-1 >*?:>>» 

j?i«TMrNnrv.niv.7?^i» — 

17i-1MIN0TV-( (.J?-1 >*7C>H 

WRITF IINIT(^) (K.MAVDTV, IMAXOTV, T I * J1 « MTNni V, IMINDIV. I2»J2» 

^fSX^TJS.FlS.Pl.^t-i.FlfS.S.^irSIM 

MVAi.Tr TMAxnrvji -0. ot 

MYALTCIMTNDtVll-0.0* — — 

( f^ifFi 


FNi.iF» 

CL0SF('F1^.F1 )( 

Cl 0«F<''F2^.e2>i- ’ 

CI.OSFl -^Fa'.FaM 

CL0:>F<-'F4".F4)I 

:_Cl OSF( <FI^'»Ff.M 

CI.O!aF.C^F6.^.F(l>H 

a OSF(-<COrMTRX '»FftM 

<;|.OSF( ■'Sni t1TRX<' ,E9)l 


a nSF( 'An..(MTRX''-,F 10 )l 

/*LFE WRITE IINTT(<S) (FI .R2.E3,F4»Er5.E6.Fa 

,F9.FJ0H»^Cl CjSF", 9(3XI?)4tLFF*/ 

WRITE UNIT(-SMFl,R2.F3.F4.F5.F6»E8 

»F9.F10H*" CUCI3F', -9(3Xr?>R4 

4jH.FF*/‘ WRITE UNIT(^)« ■ IL.EAVR MAIN PRnCRni.«RR '*M /*LRR*/l 

FNnPt /*MAIN -I FF*/ ' 


/♦*****♦*#♦**♦***•«•♦«•*#»*«•*#**«•♦♦#»#»**♦***«.**»♦«.*«■*♦♦#♦#*.♦♦♦**♦*»****/ 

/*«**«»»»«««*«**»**»«»«*»»**»«*»«»««»«»««**«•*«*«*.>«»»««»**»»**»«.«*««/ . — 


_»PAOE» ^ 


PROrpni IRF - l.'i I VRONC F » — 

-_4* RRVtSIO — !J1.1NR ft. 1931- -RI./RARIO #/ 

/♦ TRFAT- FnOF KOnFft FOR VRI (4.3.10) «•/ 

/-» K»l CPOlINn ROOF - */ 

/-* K«9- TROPOPAIIF.F FDOE «7 


/«i.F.E*y-.WRlTF_UNITlA.)*‘''lBFfiIN PROrFOURF niVRONrr''«» /«! FF*/ 

THIS PROCEDURE ESTIMATES THE VAI.IIR.OF RTOHT HAND Sir»E 

OF_Fi*<IIATrONS, - AT_T»:r FNH I T CAI Cl II. ATFS THF HI VFRGFWCF 

OF OBSERVED OR AD. iU'" TED WIND FIEI.D.. */ 

NL-t.»-Ol 

EDR-.*I- I ■ .1 JTCl>NUM_NS -00 

; FOR K -i«.i - to-num.-ai.t do 

; FOR T0-NUM_EW DO 

NL :» JMI.+ U 

1 1 - t ■ ■ )<f 1 ( 

12 .1- .1-11 

.. IF TIbNIIM Pul^^t TUPM r« tii 1tCMI%Tt 
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IF i;>*0 THFN 12 tm NUM.FWtFNntI 

vFi (i )rNi ) i» ».M.j,K)cin*oivi.K)r int 

VF«.(2»rNI.1-l--U<J.K)CI21*0<J.K>CI21l 

VFlOtfNI J-i- U(J.K+1 )m*DI,I.K-H if I3t 

tF-K-t-THFN 

VFL(4>tNl.3-i— <2«-UIJrl >r r 3-U<J.2> r I 3 »*l2»n(.J, 1 >ri 3 

D(.».2)cn>» 

-T r- K-9 THFN 

A^El. (4>CNL3- J.©)f I3-U(.l.9> r n J»fi) M3 

D<J.‘9)Cn>» 

ENDII 

— -TF< (K>1 ) . AWn. (K<9nTHFN 

VFl.<4>CNl.3-»«-U(.J.K-l)CI3#D(.J.K“l )Cnt 

FNnri 

FNDFt/*l.nnP ON I- -I.EF#/ 

FNriF»/#I.OOP-ON K I.FF#/ 

ENDFl /*LOOP ON-^ J— I FF*/ 

NI.-4M 0* 1 

FOR-J !■ -l_TO-NMM.Nft-nO 

for. «_i_to nijm .Ai.:r-Dn ^ 

FOR I-» *- l - TO -NI IM.F.M 00 

Nl^i --MI. + 1 : 

I- 

-J-U 

IF- OJ «NMM. NS+I -THFN Jl-l ■ NIIM_NS-1 iFNDI » 

L T F- ..12*0 ^THFN - J2-t - -2 1 FND T » • 

^vFi (5»tNi 1 :« v(jt.K)m*n(vii»K»rr3« — 

VFI.(-S)CNI.l-i» V(..I2.K) Ml*D<..l2.K)Cnt 

^VFl (7)rNL3-»»-V( J.K+1 )tr3«0( J.K+J >1131 

IF-K'-I, 

^THFN -VFI ( 8 ) r Nl . 3 -t • • -I 2*V ( ) r f 3 -V 1 ..U 2 ) r I 3 ) « ( ? « TK ..I ..1 ) T I 3 — 

-.0 ( .J, 2 >-U 1> I 



_ I 


I 


FNPir t 

lF.K-9. 

THFNi VFI. (8> TNI -3-4 « (2*V(J.8>T 1 1-V( J.-9-) f l-l ) * ( 2« 01 8 ) f I 3- 

jrri(..l. 9.) C 13) I 

FNnil 


JF.UKVt >.,AND.-<k'<9> ) -THFN 

VFI (8)rWI..3-i«-VUI.K'l )M3*n(.I.K-l >ri3« 

END II 

FWOFt /*l finp. fiN -I 1 FF*/ 

FNHFt /*( OOP ON K-.- I.F.F*/ 

ENOFt /«L.fjOP ON J -1,FF»/ 

Nl— *« .01 


-FOR -.1-4 ■ - 1 _T 0 Nl IM..N8 OO 







vEuiocNi.i t= (?#M(j, nrn-w(J.2>rn)*(7*n( j. nrn 


IF K»9 THEN 


-IF( (K>1 ).ANn. (k'<9) >THFN 

-VFL ( 1 0 > C NL ]-* * -W ( J rK~ l-K 13 *n ( d , K- 1 ) C I l 1 

-FNDH 

VFL( 11 ) ENL3 i» U( J»K) r I3»n( J,K)n 3 t 

VEM 1 2 > CNL3 : » V < K ) C 13*D < J. K ) T 1 3 1 

VFI. n:3>CNL3 s»-W(J»K)n3*rU J,K)r 131 


_ . __ 

— r.ixM.ir i — - — 







/» -f:AI.CUl ATING RIGHT- HAND SIDE OF EQUATIONS 

i^riiFU TT tK rt fSt,, - - 


y* 

1 - RFnTlsl rriMMFMT HFRF it-/ 






i • ru Ti • — 1 • ru ai Tii- 



•• iT 1 'JI* — * O l. 1 •• TT. ITIO»«.rnL . 1 

SSX < I > C-t : EW-.NS_ALT3*VEL( 1 ) C 1 : EW„NS_ALT3 



F.NDF; : 

. / *_^^-_Fr4D -COMMENT HFRF */ 

4*-CALCULATIN0 DIVERGFNCE-OF- WIND FIELD-*/ 

AL?—: 0*AL 1 *AU 1-; 

D I V C 1 : EW _NS...ALT 3 : = ^9 C I : EW..N3 .Al T 3 / AL3 ; = 

MAXDIV.-i=_MAX(DIV); 

MINOJU_:-=-MIN(DIV).‘». 

/*l FF WRITE- LINIT(^) (MAXDTV.MINniV)41'’MAXDIV -= -'-,E14.7, 

Oil ND I V- E 14.7#? I .EE* / 

WR I TF-UN I-T ( 6 ) ( MAXD J V , M I ND I V ) # <-MAXD I V- - - F 1-4 . 7, 

i_.MIN0tV-^=-^,El4.7#; 

BBB: WRITE- UNIT,( 6 )#cn.FAVF- PROCEDURE -PI VRONCF^’# 

ENDP ; / *D I VRONCE-LEE* / 

__Z****************************************************«*********««****/-- 

—.y.***********************.*********************************************/- 
__/*********************************************** *********************/ - 
__/-*******♦*********************************♦**************»*********»*/ - 

/♦.♦•»:**•****#*.*»*******#***♦****«•*********#***«**«•*»*****♦****«****«***/..... 

*PAGE? - 

PROCEDLIRF,- SOR? 

INTEGER-JMOD, KMOD.-IMOD? 

/.*LEE* /. . WR I TE . UN I T ( 6 ) tt ^ 1 BEG I N PROCEDURE-SOR " # ; / *LEF* / 

! A*^ THIS PROCEDURE IS A PROGRAM FOR A SUCCESSIVE -OVER - 

! RELAXATION (S-O-Rl . ITFRATION_METHOD TO SOLVE SIMULTANEOUS 

I EQUATIONS. FIRST VELOCITIES- AND - LAMDA -TERMS ARE ARRANGED 

j IN-VECTOR FORMS, -AND THFN S-O-R EQUATION IN GENERAL 

i F0RM_LS -U3ED FOR -ITERATIONS.- */ 



28 


NUM.ITFR l« 0! I 

rNTTIALIZI Nn Tl AMDA- ARRAYS • */ ' 

FOR T I- I TO EW..NS„ALT DO 

^TiAMnAn>riii5] I- o.oi 

^/*l.FF NOT. rONV I « It*/ — ; 

NOT..CONV l« Ol — i 

FOR -t.1 - 1 - - 1 TO Nl IM . NS • 00 

FOR K' i» -t TO NI.IM . AI.T DO 

FOR r «“ 1 TO NIIM_FW 00 - " 

Nl. -!■ Nl. + t» 

II ,« i + i t - 

i3..,._r+?i 

TR.-r-l=NI.IM_FW+l -THFN I l--t « -1 t FNOn | 

IF i;>=«0 TMFN \2 t» NUM..FWSF.N0I5 !• I 

- IF. <I1=NI.IM,.NS+I THFN J1 t* Nl IM..NS-J ; l-NDH 

IF-»I2«0 - ..-THFN .J2 :« 2?FNntt - 

IF Kl=NilM_ALT'tl --THFN Kt t« NUM-.AI.T; FNDIt 

IF -K3*0 THFN K3 : = It FNDT I - 

IF -I3=NIIM. FlJ+2 THFN IS-t* 7tFNni» 

I F . I ;-s=NI.IM .F W+ 1 THFN 1 3 : » 1 1 FNO 1 1 

IF I4«=0 THFN 14 ;= NHM..FW; PNOI t 

IF-I4=-l — THFN 14 t» NUN FW-ltFNDIt 

TlAMOACNl ..) r 1 3-1* W AMOA(J.K)l lit 

TI.AMDA(NI.K21 : « Pl.AMOAt .J.k > C 1 11 1 

TL AMOA ( NL ) l<i 3 . -* « .01. ANOA 1 J . k ) 1 1 ? 3 t 

T I.AMOA ( Nl . >141 -S = -DLAMOA ( ..U . k U 13 « 

TLAMOAINL 1X51. -01 AMOA ( k ) 1 I 3 t 

TLAMnA(NL)C61 : =» .DLAMDA(..l,kl Ullt 

.ILAMOAINI.) 173 - PI AMDAl J.kS) 13 3 t — 

TLAMDA(Nl.)C81 ta RLiVlOAl J, k U X 1 13 1 

-Tl AMOA (Nl.) 193 I 01. AMOA (J,k3 > til 3 t 

- -Tl . AMOA (Nl)t 101 :a Ol.AMDAl J.kl )CI 21t 

-Tl AMOA I NL ) C 113 i a Pl.AMDA ( J.kO) ! 1231 

-TI.AMnAlNI )rr?l :« BLAMnAl.Jl ,k1 )CI3t 

TL AMIiA ( NL ) r I 3 J » B 01 AMOA UH , k3 ) III t - 


TI.AMDA(NI.>C141-i= ni,AMnA(.l2,kl ) 1 1 Jl 

Tl AMOA ( Nl ) r 1 5 1 -X » 01 AMOA ( J7 . k3 11131 

/* FNniNO OF- VFOTORIZATION */ 

/« .PFr.lNl NO OF- S-O-R HFTHOO */ 

-T.:= TI.AMOAIND.nOT.TXlNI )t - 

PI AMHAUl.k ) r I 3 -xa Al.AMOAl J.kir n + (Rl XFAC*(SrNl 3-T)/ 

TX.tNDCl])! 


/♦ DIAONOSTIC FVFRY lOTH GUY */ 


0WODIALW65J 

OF POOR QUAUn 
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MARCH -SV. PLR */ 

OMODi-.I.MOn. tOl 

IMODi«I.MOn.tO» 

KMCnt-K,MOn, lOl - — 

1F( (JMOn-t ) ,ANn. (KMOn-1 ) .and. ( TMOD-I > > THFN 

WRITR »IN!T(6» 

—i 1 , K' . Al AMHA ( J , K' > r n . PI AMDA ( J, K ) CI3 . SC NL 1 . 

T.TX<NI.>m)1*3T5,5Elf5.51M 

IF( (vl»1 > ,riR.~( J-NUM.NS) > THFN 

THnN ..U t- 3tFNnXt 

TF( J"»NnM..NS> THFN J J I « NUM ,NS-I « FNCIt 

. amDA ( .J, K ) C I 3 I ■ SUM < BLAMPA ( .J.J . K ) C t « NUM.FW 1 ) /Nl IM ..FW t - - 

/# UP HY KOF-HXS MAR 31 -«•/ 

^TRY A"FlOAT" */ ' - 

/♦ PI AMnA(5,5) C3/S3 •« 1.5F+06« */ - - 

TF(AI-AMnA<..UK)r n-0.0) THFN ALAMDAX.I, K ) C T 1 l- O.OOOttFNHTi 

FRROR I- 1 .0-(BLAMDA(vl.K>rn/AI AMDA< J, K) f T 3 > » - 

TF APS (FRROR )>TOl.F.RANCF THFN NOT_CONV I » NOT..CONV+ 1 ; FNOn 

/♦ XPMPORARY RH3 PRINTOUT - 

PLR MARPH \2 1931 */ 

. - for 0 FROM 1 RY 100 TO FW. NS-ALT DO - - 

WRITE- UNIT (A) (0, SCOT )#SX» 15, 5X,F.15^iM 

FOR .1 :«-l-TO NUM NS HO > - 

Wl AMDA CUT s» PLAMDA<..I,5) C3/>3 1 - - 

/*l.FFfr/WRITF IINITIMtt 'OWI AMDA -ARFi 'PI /*l.FF*/ 

WRITE- UNIT(*1 (Wl AMDACl :NLIM-NS3 )P5F14.7P? 

for K »= 1 TO NUM. .AIT DO 

KLAMDACkl t«-BLAMDA<5,K)C361*- — - 

„./*l FF WRITF IINIT(A) (kl AMDAC 1 : NUM, Al T 1 » P '’Kl . AMDA , 5FM . 7P| 

WRITE IINIT(4> (Bl.AMDA(5,5>C5O;403)P I.AMDA =» ^5El4.7Pl - 

WRITF -UNIT (M ( PI. AMDA ( 5,5 1 1 13 , PI. AMDA( 5, 5 1 1 NUM. FW3 > - 

P I. AMDA » :.?EI4.7PUEE*/ - 

! WRITE UNITIM (Kl AMDAC 1» NUM. ALTl IP/. KT AMDA * '1, 5F 1 4 . 7P I . 

WRITF UNIT(6) (BI.AMDA(5,5) C30I403 )P • l.AMDA » . 5F14. 7Pt — ■ — j 

WRITF IINITIM (BI.AMDA(5, sun. PI AMDA(5. 5) CNIIM, FWI ) i 

P'. l.AMDA «_Lv2F.l4.7P» ] 

MMMM..t».(XJ 1 

FOR J_J«=^..1T0 NUM.NS DO. ] 

F.OR TO NUM..AI.T.DO — . 

WLAMDACNl.3 t ■ MAX ( PI AMDA(. I, K II ; J 

Wl.ANDAirNM «■ MtWIRI AMnfti.i.n »» ] 



30 


• •«••• tl «Wt * * • I' 

IF (MMMMO) THFN 

TF (Wl AMPAJ rNU<0.0)- THFN * 

FF WRTTF »INTT(6) (..UK.WI AMOAt tNLl ■ 

s'W( amPAI •- 'SFlA.TMtl FF*/ , 

WRTTF IINinf » UNK.WI AMDAl fNl 3 . 13. ' K 53. 

-'-UI.AMDAl----',F14.7Hl 

MMMM «■ MMMM+1 « ^ 

FNHU 

TF( (il«37» . AMD, (k-n ) THFN 1 

/H FF WRTTF UNTT(/S) (.I.K.WLAMPAICNI.IHJ'.I « I3. K » •",13. : 

AHPAI - -'.F14, 7AU.FF*/ ' 

WRITF IIN1T(#.H.I.K,WI AMHAl TNU )4' J - '.13,^ K ■ '.33. j 

' WLAHDAl =-'"-,Ft4.7*3l 

Al AMnA(J.k)ritNUM.EW3 -l» PI.AMHAT J.kM 1 s NUM.FWl ? ' ^ 

MAXI AMPA J» MAX(W1AMPAM 

— MINI AMPA ■*» mN(WI,AMPAl )% ; ; 

A*l FFWRTTF-.I,INI-T(<%> < MAXI AMPA. MINI. AMPA M^^^kAMPA s-^-r^'MAX- ■ -^.ElA,-?. 

_...-MTN » •■'.Ft4.7«*.I.FF«/ • — — 

WRITF UNITIf.) (MAXI, AMPA. MINI AMPAM»<-LANPA *s - "-MAX - , F 1 4. 7 | 

NIIM.-TTFR NIIM,. ITFR-'*- 11 ^ ii 

y«l FF WRITF-IINIT(vk) ( N0T..rpNV ) »-'NUMPFR OF Fl.FMFNTS — NOT. rONWH- c— ,j 

WRITF UNTTIM (NOT CONV.NUM .ITFR) I 

’4-1 , . TO. ■ FLFMFNTS HAVF.-NOT CONVKROFP AFTER-' . ] 

IB, ITFRATTONS. -'I*? - - | 

/*! FF - WRTTF lINTTIM INIIM. ITFRIA- COUNTING OF I TFRAT IONS . lOtt; I FF.* / i 

/*| FF IFTNIIM TTFRIMAX.. ITFRTTHFN 00 TO RPP: FNPT ? I FF*/ - - 1 

^I.F(NUM .-ITrR>MAX ITFRlTHFIvl GO-TO-BBP* FNPI 5 ' j 

>:»! FF -IF NOT t'ONVM THFN 00 TO I ClFNOU#/- | 

IF..NOT nON'y > 0 THEN 00 TO LCt-ENPH 

/«! FF- WRTTF UNIT(6M ITFRATT0N)*»1H0, 'NUMPFR OF- ITERATIONS » | 

WRITF I.INITfM INIIM.ITFRIH'O-"-.— 13.— -'ITFRATIDN'o iM-.FORF TF'RM I NAT ION' , ] 

£jc\f s-o -r sou .it I on -' 4 » 

/* FNIITNO OF S-O-R NFTMOP « / 

/« .PRINTING OF- l.AMDA TERMS */ | 

I J* FOR J-»« 1 TO NI IM NS PO •« / ; 

— /*! FF*/WRITF-UNIT-<^)4 OFTNAI.. I AMBPA ARF« ■4t/*l.F.F.*/ 

-....FOR J :«_1 BY NS..INC.R -TO-NIIM .NS no 

FOR k_: «- \ - BY Al T TNGR -TO NI.IM ALT -DO 

WRITF UNITIB) (.I.K )4)HO. r..l-= '‘.13, 'k » ''.I34I ; 

- WR I TF- J.IN I T ( O ) ( Dl .AMPA < ..I. K U 1 : 1 0 1 > «5F 1 4 , 7 4 « j 

__/* APJl ISTMFNT-OF OBSERVED -WIND-FIEkD-TO- MAKF-J-I-MAaS-- I 

1 CONSERVATIVE */ ! 








Nl, l» 0» 

FfiR-tl-t •-1-Tfi- NMM-NS HO 

— FOR K l«-t TO NI.IM Al.T 00 

FOR T t«-l -TO NUM.FW PO - 


NIr-i- -Nl>n 

APjr H5i t» ADJXtNi. )r H5)t 

— — . It- — .1 .1 I. ■ 

IF ( n-NMM.FM+n THFN T1 l—lJFNnil 

IF (T?»0 ) THFN 12 t» NI.IM .FWt FNDU 

IF THFN ^11 NMM_NS-1 »»!?«» NUM_NSiS FNPI » j 

IF (.I2»0 ) THFN »■ UFNDH - - I 

IF H 1»NIIM..ALT+J ' -THFN K1 l« NnM..Al T« FNPI I - | 

IF (K:^»0 -) THFN KO *» llFNDl? - — j 

AD.iusTiK^i.K) t n -«= IK j.k) cn+APjr 1 3*«Pi AMPA( vKk )r I n ~ 

-Bl AMOA<..UK)CI23)-ADJr 21*(BI-AHnA(J,Kl )Cn ^ 

m. AMPA ( . I, K?s ) r 1 .1 M 

AP..iusTV( j.k) r n v(v».k)rn+An.ic3HMPi.AhPA(ji.K)tn - - 

^-RI.AMnA(J2-KUn)-AD.JC41*(nLAMDA(.I.Kl>Cn — > — 

PI.AMPA(J,K3)[m» - 

An.lUf5TW(^Kk‘)r M I- W( J.K)m+APJC51*(PL AMDA(vKkl ) r n 

PI AMDA(.J,K';'3>cn>? - - 


__fndf< 

FNOF? - 

PR I NT I NO OF AP.IUSTFP WTNP FIFI..P «•/ 

/* FOR ,1 s=-t NUM_NS DO -- */ 

FOR J iB„t. PV-NB, INf;R TO NIIM_NS PO - 

F0R--IC~:= 1 BY-Al-T-.INCR-TO NUN.Al.TPl PO 

/*l FF WRITF MNTT(6» ( J.k»41H0, -'J = '•.IB. -'K - 

UR ITF 1 IN I T < O ) ( l.l ( . I - k > C 1 : 1 0 1 , AD..IUSTU ( J , K ) C 1 : 1 01 . 

V'(J.k)f1t 101 , AnjMOTV(.I.K) r U 10], • 

u ( . . 1 , k 1 r 1 : 1 o 1 , AO. II loT w ( . I , k ) r u 1 o 1 ) 1 4 . 7 4 ? ' FF* / 

WRITF- llNTT (/SHJ.k )»' OJ « - < , 1 3, J‘k » 'VI3IM - 

URTTF I.INIT(/>)(l.l(..I.K)Cl:FW..0UT1)#' UNAn.JIJOTFn II ARF‘ , 5F1 4. 7«1 

WRITF IINITIM ( AD.III:-.TU( J,k ) C 1 tFW OUT] >4" AMJMOTFn U ARF ' , BF 1 4 . 7lt I 

WRITF -I.INTT(6) (V(.l,k )r 1 :FW 0UT1 >4 UNAP.IUSTFO V ARF ‘ 1 4. 74l 

WRITF IINITIM (AP.IIISTOI.I.Kir 1 !FW 0UTH4' AP.IUOTFP V ARf ' , f<F 1 4 . 741 

URTTF I.INTTIA) <W(..I,KK 1 :FW.0IIT1)4' MNAO.IUSTFn W ARF .r:.K14.74l - 

WRITF LINIT(6) t APJUSTWI J,K) C 1 1FW..0UT] )4i'- APJIISTFn U ARF ' , BF 1 4 . 7« 1 

FNOFl 





i?.AU nivRnNinFj 

-/♦I Ft- -RRR: FNDPU FF*/ 

— FF*/ 

-BRBl-WRrTF FAVF PROnFIIURF- SOR'«M 

FNnP? /*<5i1R LFF»/ 

— /♦IFF*/ 

FNIim 

.. F»T»T-FNCni ll\|TFRFn. 




FILE WINOP 


ive-vhii ndp 



IWrNOxllINnP 

o pv 5 h F » 111 1 n »i p 



I 


« 


-31,. 


MOrmi F WTNn(I.INITl»ril.PTRnP,I.INrT4«C0FMTRX.UNIT7aS0l MTRX 

r«.IN I TB« AD..IMTRX . UN I T6-RFSI H.T’. UN I T9-S0LMTR ) 

PFOIN MAINI i 

-y*cr: PROORAMMF WIND I FADS TO AN AnjI.ISTEO MARS -CONSERVATIVE — CCCC 0002- 

•C WTNO FIELD FROM RF.AL-T IMF* ODSFRVFD WIND FIELD DATA.-THIS C 0003 * 

-ft PROGRAMMF - 1 NVOI. VFS "NUMFR I CAL SOI I IT I ONS - OF 'THF PARTI Al. C'— 0004 - 

0 DIFFERENTIAL- FOUATIDNS RFSUI.TINO FROM -A CONSTRAINED C * 0005 - 

n VAR TATI ONAI. -• FORM! U AT I ON - ( FI U. FR-I AORANGF TYPE ) TH I S PROGRAMME C * - *0006 • - 

-CCCn IS- WR I TTFN -IN STAR LANGUAGF.-l < SL I )■ OP CYGFR * COMPUTER . CCCC 0007 - 

CCCC OFF I N I T I ON - OF- TFPMS — CCCC 0008 — 

-e ^X I "EASTWARD DIRECTION C 0009 - 

-e Y!«NORTHWARD niRFCTION C 00)0 — 

-e -7l«VFRTICAL -DIRECT ION C 001 1 - 

G XX:=GRin-SI7F-IN X-PIRECTION, DFGRFFS RAPIANS RESPECTIVELY C 00)2— 

-G VYiaGRID-SI7E-IN Y-D I RECTI ON, DEGREES -«<- RAD IANS RESPECTIVELY C 0013- 

-G 7 Z I aGR I D S I 7F-I N Z -D T RFCT I ON 0 - 00 1 4 -! 

-G -IL:-aNUMBERS-OF GRIDS IN X- DIRECTION C- — 0015 — 

-G ..ll.taNUMPERS OF GRIDS -IN Y—PIRFCTION C OOIG - 

-G KI.*aNUMPERS-OF GRIDS IN -Z— DIRECTION C — 0017- < 

-G -THFTA-«/-TH -t aANGl.E-IN DEGREFS RAD IANS- RESPECTIVE!. Y FROM SOUTH C 0018 — 

-G PO! E ( “90. 0 > -TO- NORTH- POI-.E < +90 . 0 > C 0019 - 

x:: Nl AMDAt aDIIMMY LAGRANGE MULTIPLIER FORCOUNTING PURPOSE C - 0020 ■ 

T ALPHA : =THREE D I MENS I ONAI .COEFFI C I ENT MATR IX 

.CCCC SSJaONF DIMENSIONAL SOLUTION MATRIX 0023 --J 

THE FOLLOWING -SECTION -NFEPS-TO- BE -CHANGED FOR 1 

TUFFFRFNT-GR ID -SYSTEM. XTRAl- T0-ISM3 -ARF-PAGING 1 


PARAMFTRFS— TH I FR VALUES -SHOULD -BE ESTIMATED 

ACCORD T NGl -Y. -AI ) ,Al 3 AND • AL2 ARE -WEIGHTING FACTORS 

FOR -U,V-AND_W .VELOCITY- COMPONENETS- RESPECTIVELY.-.*/- 

PEG I N I NG OF-THF^ SFCT I ON * / 

.l.TTFRAl .LY- 0 '• -DX X 1 — 

-LITERAL! Y.-'i5. 0"- DYYl 1 ^ 

ajt.TERAI L« 

1_ITFRA! )-V— "37 "- JU 

-L4TFRALLY~i'9.» -Kl t 


I^ITFRA! l-Y—t'-) 0''k KL!U -/*k'Kl !.aKL+1 */ 

L4 TFRAI . LY-.'12r-!974. "- -I .IK’.1 

laTFRALLY—l’0 . 0 J-I!-Al,l-» 

UP[FRAU.X-Ilx-).01 i!_Alv3t 

L-l-TERA!..l.-Y_.'ll ^ fi " .. Al. 2 1 

l_LTERALLY.JL3.357/s"_XTRAl ? 

l-lTFRALLY-.L')fi992''— XTRA2I 

LITERALLY-IU 1 192.'1_ XTRA3; 

i-J.TERALLY— IL7SP." . I .SMl I 

LJTERAUJY_11640!'_JSM2t 

i_IX£RALLY_!l2f.S.".. ISM3! 

y* -END OF _THF -SECT I ON -*/ 

REAL. VECTOR [ kia_L3_7.t 

1 N l.TJ AL f Z a \ 0 . f)33.3.?.3A , 0 . 0-33.3334 , 0 . 03-33.334 , 0 . 0333.334 

, 0 ,-U 0 . 2 , 0 . 2 . 0 . 2 , 0 . 2 , 0 .-2\ ) ; 



LITFRALl.V "ZOFCTriRrfc'T' 77, 


ORIGINAL PAGE IS 
OF POOR QUALITY 
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.TTFRAI.LY "Z7trKl" 71 1 

rl TFRAI I V-« 7 Z ^»r K 3 •*- 7 ? 1 

-TTRRAI. lv -"ZZrjtKV "73t 

^FAI— VFHTOR fKl-3-771 T 772r773i ZVFCTORt 

ir ITRRAI-.L-Y— 59v.o535R9793‘* - PI » 

irTTFRAU Y-"MRITF“- PRINT* 

»-T TFR Al . IrY-" RTRR C-t-) H I t T I vV* - T3 1 » 

i -TTFRALL Y-«»RTRF: ( ? > M « 1 1 L 3 " - TS3 1 

l-l-TFRAI. ItY— ''RTRR (-3 > F 1 1 F I »_*3 " -TS3* 

l.-TTFRA«J.Y— "RTRS ( 4 ) 1 1 1 T TL3 'TS4t 00 

t-1-TFRAI.I.Y— ''RTRS(5) r J t iri. 3" 

FrTTFRAULY-"R(.J*K)C t« TTL3"R0l 0040 

y* OFF TNIT IONS 

1 «l AHOA ( , I, K ; C n 1 T?»L AMDA ( J . K > r I + 1 3 , T3-LAMDA < J ^ K ) 1 1 - 1 3 

T4-I.AMDA( J+l . K ) C n . T5-LAM0A (il-t-rK) C l 1 , T6«I.AM0A( .J, Kf D C I 3 

— _T7»L AMriA(..l,K-l ) r T 3 ,T8«LAMDA( J.K-M )CI + 1 3«-T9»l.AMDA( J.K-l )f I + l 3 

T-10*LAMOA(J,K-t-l )CI-13,TU»LAMnA(.J,K-UCI-l3 

—^1 ?»l.AMDA ( ,1+ 1 , K+ 1 ) r I 3 T T 1 3-L ANOA < J+ 1 . K-1 > f-I 3 » f 

T 1 4=1 . AMOA < , K+ 1-) m , T 1 !f.-sl AMOA ( ..I- 1 » K-1 ITT 1 

RTRR»niSTANCF- RFTUFFN -SI IRF-ACF- ANO TROPOPAI ISF -AT- ANY - CiRID. 

OAMAI l-^AMAV -AND-OAMAU' -ARR -COEFFICIENTS -FOR -THE 

ADJUSTING FOI lAT I ONS. « , ^ 

l.TTFRAU.V -I'-TALPHAC 1 3^T14 ^ 

11 TERAUl. V--" T Al PHA C 

L-I TFRAI.LY -"-T Al PHA C 3 3 " — T3 1 

L I TFRALUY-y TAl PHA F 4 3 -T4 » 

>_T_TPRAI..IrY -"-TALPHAC 5 3 T51-: — 

_4_j tFRAI LY "-T Al PHA T t> 3 TS * 

LT TFRAI I Y- "T ALPHA C 7 3 T7 F— — 


_<_I_TFRAI ;L Y - " TAI..PH A r S 3 "-T8 1 — 

I-TTFRAI I.Y— "TALPHAC93 "-T91 

LTTERALL Y— "TAL-PHAf 103 "-T10» 

1 T TFRAI :LY~"T ALPHA! 1 1 I"— TIH 

I^I-TFRAl |.;Y -"TAl.PHA F 1 2 3 " -T 1 2 1 

J..I-TERAU.Y— '\tAI.PHAC1 3 3 T-1 3» 

l,I-TFRAI I -Y- -''TALPHAt 1 43 'L-T14* 

U-TERAIJ _Y..." TAL PHA T 1 .=L3 •-'_T-1 5 ; 

LU'FRAI.LY-J.'GAMAUU-) F J : I Il.3i'-GUl t 

I.ITERALUY—S'GAMAIJ ( Z) C4 : I-II.3 ’i-GI 

Ul-TFRAl. UY_"GAMAV (-1 ) F I : I-I I.-3 i' -GVl * 

LITERALLY— '-'0AM AV (-2.) C I *-I ILl-'l-0V2t 

J ITFRAI LY- -'JGAMAWI4 LF-t » Ml -l-!!-GWM 

_LT-TFRAI.I,Y— "OAMALU 1 ) C T l"OriUl * 

-LITERALI-Y-'IOAMAU ( 2 ) F I 3 "GGI 12 J 

.JJTFRAI.14'_';0AMAV( l-> C I 3^'OGVU 

_XJTrRAl.Ly _"OAMAV ( 21 C I 3^'GGV2* 

_L LTER Al .LY _'iO AM AUl < 1 1 C-I 3 " GOU t ? 

-EXTERNAL- PROCEDI.IRFI INTFGFR.BI T-VFCT0RC64 3 . REAL-, 

I NTFOFR . I NTFGFR ) - Q30PNMAP » 

EXTERNAl- PRCiCFDUREU NTEOER. BI.T_ VECTORCi.4 3 l-Q^CUCiSE-L 
RFAI. VFCTOR C153 ARRAY! UK) TXl 
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COMMON CT(TX, THUMP 1)1 

REAl— VICTOR- C XTRA11 - TOUMPl I 

RFAl- VECTOR -FI 3 5 - ARRAY < I JK)-SX I 

COMMON CR ( sx ♦ thump:? > I 

RFAL -VECTOR C XTRA?3 •TnUMP?» 

RFAI- VECTOR C53 ARRAY( lUK)- AHJXf 

COMMON CA ( AO.JX . THI IMP3 ) 1 

RFAI, VFCTOR-CXTRA31 -THUMP3f 

RFAI— VECTOR -r T JKl - ARPAY( 1 31-SSXl 

COMMON CSS(^:SXrTHUMPA)t 

RFAL VF.CTOR rXTRA^J THUMP4I 

RFAI . VECTOR CUT ARRAY < - OAMAI J , OAMAV . 0AM AW X 

RFAI- VECTOR C53 AHJt 

TNTEOER -lERl, IFR:?, IER3, IER4I 

INTFOFR-f.^UK-,NI , I I T 4, cM K3» KKvJJ. K?* I II.. COUNTr“- 

REAL -VECTOR TM RRO» 

RFAI - AI..4.CON,7Z7, XX, THFTA,TH,YY I * 

RFAI - VECTOR r ILL ARRAY(O) ALPHA, PFTAl 

RFAI . - VECTOR - C-T I . T ARRAY ( 5 ) RTRS I 

RFAI.-- VECTOR - r 15 1 -TAl PHA X 

RF Air VECTOR f lOl-SSl 

__rFAU.-VFGTOR C II-3-A, P, C» Ft 

RF Al:r-VFCTOR C U -3 -ARR A Y ( Jl. , KL. ) R I 

RFAI . -VEGTOR - E - II 3 TRAP I 

INTT4AI-(TR0P«V9. a,P.3,9,4,9.5,9.7,9v9, 10. 1, 10.3, 10.5, 10.7, 

40.9,1-1. 1, U.a, 16.8, 17.8, 10. 0, 18.0, 18.0, 18. ( 

i!^0 ,-18,0, 18. 0, 17-rR, 46. B, 1 Werl Iv l r lO. 9, lOi l 

1 10.3, 10. l,9.9,9.7,9.5,9ir4r9-.-3,9.9.\)l 

RFAU VFCT0RCU.3RR1 *RR2,RR3,RR4,RR5,RR6l 

PROCKOI IRE-MAT Nl 

WR I TF— I IN I T ( 6 ) ( HX X . HYY , Kl . , T ,.iK , AL 1 , AI..2 , AI.8 ) 4 ' OX X - ^ , F5 . 1 

<-HY Y« •' , F5 . 2 , < KL»''-, 1 2 , ^ I«.IK« < , 1 7 AU 1 - f 

■_sA1.2='‘., F7,3. 'AL3«<-, F7 . 34 1 

7VECT0RI: UKL-3: *< 7C2t Kk'U 3+7C I *KL1 )/2.0l 

27 1 r 1 -» Kl -3 t » 1 . 0 /-7 1? : KKl. 1-3 1 

7-Z3ri:KI.3t = t #0/7C ItKLlI 

7.72f1 iKl 3s*=7ri n ;KL3+773ri»Kl.3t 

C-Al.l. -Q3nPNMAP( TFRl , X\434E464n54525820\, TX ( I ) E 11 , ISI 

^CAU.-03nPNMAP I I FR2 . X V534F4C4ri545?2020\ , pX I 1 ) C 1 1 , 1 

CALL- 030PNMAP( TFR3, X\41 444A4054525820\ , AH..IX ( 1 ) C 11 , 

WRITE-UNI T(6). LI FRl ,-IFR2, IFR3)4' JFRle^^, 716, 4X, 

<A ER2*--'-, 7 1 6 , 4X , I ER3« ^,71641 

COl.INT t »0 ? 

CON S «P.I / 1 80 . 0 1 

y»_INT.TIAl I7IN(i. TALPHA,An,.l,8S- ARRAYS */ 

^TAI PHACl ! 153:«O.Ot 

ADJLl : 5>l»0.0l 

S8CJ : I31t=0.0; 

XX : «nx X«C0N; 

yY:=DYY*CQN; 

K2:=KI +1 : 


0043- 


0, 

7, lO.Sr- 


2 , 

F7.3, - 


0051 


Ml , I ) t — 
M2, 1 ) I — 
!SM3,l)t 


OWGtSXtVAOBC 

pch)R 


Al. 4t -(AM /Al .?)»»?« 

H(. I "n-t 

-V*CC — INITIAL-f7rN0 Fl.FMFNTS- AF R( ri.'3 

— FOR *vji - 1-To di. no 

RFAD-HNITd »-IR(^l*-I)M» II. .i >- 3. <%(♦! 

R(J. l>Ctt II.-ll»<S371000-0-+R(.J, t)Clt ILlI 

R(.:.KL >t It IL 3 I -6371 000.0 --+TR0PC J3*1000. Ol 

7721-0.01 

FOR KI-2 T0 Kl-l -nO 

-Z7Z»-277+7tK3l 

TSli-R(.I.KI Kli II.3-RCJ. 1 ITU II. 3 I — 

ROi-RM. I )Cll II. 1+ZZZ»TSU 

FNOFl 

FNDFl 

_/.#CC THF-FOl.LOWINO-DO-l OOP-ltFFINFS-THF FI.FMFNTS OF COFFFICIFNT CCCC— 

-C MATRIX "T ALPHA "-AND SOLUTION MATRI X "SS'MT HAS FIVE C 

-C PROCFPURFS WHICH OF AL WITH FQUATIONS SUITABI.F FOR PARTICULAR— C 

-C CONOT T TONS. PROCFOURF - "SNPOI .S'J- STANDS FOR- SOUTH NORTH POLAR 0-- 

-C SURFACF.SIMILIARI Y ''SNPOUP'- STANDS FOR SOUTH- NORTH POLAR C- 

_C INTFRIORt SNPOLT FOR SOUTH NORTH POl AR TROPOPAUSFl i C - 

-n ?•' INTFPS'*- FOR- INTFRIOR-PCU NTS AND -SURFACE •ANn~"-TROPOI-”-FOR G- 

-CCCC™ TROPOPAUSF- INTERIOR. CC*/ - 

THFTAl —90. O^DYYl 

FOR-<ll«.l.-TO ..IL-DO 1 

THFTA* «...THFTA+nYY< 

-TH: cTHFTA-« CON \ 

FOR t - TO -KU-DO 

TSl :»RU.Kl >r H II 3-RI J, IK 1 » II 3« r. 

IF-LUI-1 ).0R. (J-.IM I.AND. IK'«1 )— THFN-CAI L— SNPOI S» -FNDT » 

IF < ( J»l ) . OR^I. I-JL ) > . AND. ( (KM-4-. AND.-(K<KL->->— THEN CAl L— SNPOLI I- 

ENDI; ^ 

IF- (1. 1-1 ),0R. (U-JL- > ) .AND. (k-KL) THEN- CAL I-- SNPOI. T I ENDIl 

IF— < (.J>1 L.-AND. (JCJL) ).ANn. (K<KL) -THEN CALL- INTEPSl-ENDK 

^IF— ( (U>1 V.AND. I.K.IM ). AND. ( k-KI » -THFN-CALl.-TROPOI «- ENDH 

ENOFl 

FNDF »- 

CALl.-03r.LOSE.-C IFR 1 . X \434F464DS4525S^:0\ >-l 

CALL- 03CI n.SF-( IER3. X \4 1 444A4nf.45.J!f.S90\ ) J 

Ci^Ll...-O.30PNMAP ( TER4. X\S34F4C4nS4f59!:.S20\-.SSX ( L->Cn. ISM9:r-l ) % 

WRI-TF HNITI6) ( IFR4 Ht ' IER4-'' . 7 1 AIM 

FOR I:»UTO -13 DO 

<11 -1 1 

i^f_s:3X( nr<n*»sx(..i)ci4« 

IF- |..K=I JIO-THFN -CO-TO-l Sl-ENDH 

FNDF; 

CALU. 03CI .OSF I I FR9-.-X \ f.34F4r.4riS452?i'>90\H 

CALL J03CL0SE ( I ERA , X \.'5.34F4C4nS4S25R20\ ) I 

WRITE LINI T I A)R'^r.l.nRF''IM 

WRJ IF. -UmiX6U J ER 1 . I ER2 . IER3.-IFR4 > I* HER Jjb^7 1 6. 4X 

•'.TFR2«^. 71A. .*X. TFR.'^a''. 71 A. AV. •'TPPd«^.7lA*ki 


PNOPt 0M3 • 

PROGPOMPP -RNPOl. S I 0 1 44 

-/#CC PROCFnMRF" "SNPni.S" PFAl S WTTH THORF FOMATIPNR WHICH ARF CCCC - 014R- 

.G SlItTABI.R FCiR SOUTH AND NORTH PO' FS SURFACFS- IT FSTIMATFS C 0144- 

.G COFFFTCIFNTS-OF FAOH CORRESPONPINO “LAMPA" TFRMS ANP TRANRFFRS ~C 0147- 

CCCC THFM -TO THF WORKING MATRTCFS-''TAl.PHA'* , "SS" . * CC*/ * 01 40 * 

IF J.I -THFN FNPn 01 RO 

IF U-..IL— THFN .Ui»Jt-2l ENOTl OlfH “ 

RRini II 5i»R(.u-n .1 >ti j ii.-n-R(u. nriMLit — oin? ~ 

Al;PHA(2)rii (Al.l*Al .l/< liO*AI.2*AI 2«RTRS( I > C I * tl;5 > > + ( ( (RRt 

rn riM/(Riu. I >r 1 1 n. 3 «vv) )** 2 )»Ai.A/n .0 «rtrs< nn< ili > >»7i t 

AI.PHA(3>Clt TL3i-<Al.4*RRUll IL1>/(R(.J. I > C 1 « ILI#R< J. 1 ) Cl l TUI* 

PFTA<?)r It-FI 3»-?i O*Al.l*AU*(RRl f 1 l U 3 » / ( R ( . I. 1 ) f 1 « IS. 3 »YY ) 1 0157- 

-BFTASaiCli I».1»-7.0*AI.1*AL1 1 0160" 

OVli-1 .0/(2,0*Al.3*Al 3*R0«YVM 

0V2»«RR1C1 » II.3/<2.0*AL:?»AI..1*TS1*R0*YY*( 1,0/71)11 

C4Jt:-1.0/(2,0*AI.?*AU?*TSl*(1.0/71))l . 

FriR_U--l-TO II.- PO — ’ 0163 - 

T-l-t -Al. PHA ( 2 » r I 3 +ABS ( AI.PHA (3)113)1 

IF<..I«1> THFN T4I— AF»S(ALPHA(3)Cn)tFNPIl ^ 

JF( J«Jl ) -THFN -T5J —ABS ( AI.PHA( 3 ) f DM FNPH 

T6s —Al .PHA ( 2 > C n I - 

SSI 1 23 I — BF TA ( 2 ) f I 3 1 

S5Sri33*»PFTA(.3>CTH- — 

AO .t! ii3 : -OAMAV II ) 1) 31 

AD-.ll. 4 1 : -OAMAV ( 2 ) C T 3 1 ;; • 

AP.iC.53: -OAMAWI 1 )ri3t 

i?.ALL-TRANSFFR» 

FNIiFt 01 70 - 

FNPP» * 0171- 

PROCFniJRF-SNPOl I 72 

_/*Cf>.^TPROCFDMRF— "RNPOH" PFAI.S WTTH THOSF Fi^llATIONS WHICH ARF CCCC - 0 1 7.3 - 

_C SUITABI F FOR -SOUTH ANP NORTH POLFS INTFRIOR POINTS. -IT C 0174 

_C FST.INATFS COFFFTCTF.NTS OF FACH CORRESPOND I NO "LANOA" TERMS AN)3 -c-— 0175- 

_CCCC^.s.-=IRANr.FFRS THFM TO THE WORKl NO -MATR ir:ES-"TAl.pHA” . "SS” . CC*/ 0176 

IF-,.1»1- THFN ..Ui»l-ENni» 01 7R 

L.IF_wl«..ll . -THFN... 1,1 1 ». M.-2 ENPri 0177 

RR2C 1 : Tl. 3 1 »R ( ,1. 1-M . K ) C 1 1 II. 3-R ( ..I. K ) Cl : 1 1 1 1 01 SO 

ALPHA! mi: II. 1:b(AI 4*RR2f 1 : Tl 3t7l )/( l.O«RTRSI 1 >f 1 : IL 3*RO*R0*YY M 

Al.PHA(2) Cl t Tl^l:«AI 4/(R0*R0*YY) 1 

Al PHA(3)ri: n,j:«ABS(AI PHA( 1 )fi: II. 34 ♦Al.PHA( 2 ) C 1 1 II. 3t 

BETA(t)Ctt TU3:«-2.0*AI.1*AL1/R0» 

CiVI :«1 .0/ (?.0*AI.3«AI.3«R0« YY) : 

OV2:»RR.2C i: 11.3/ < 4.0*AL3nAI.3*R0*TSl *YY*77. > I 

OWl (4.0*AL2«AL?*TS1 *774 « 

FOR l:nJ.._TO TU. DO ’US? 

Tl:«+Al PHA(.3)ri3t 

T6: —AIHS! AI.PHA ( 1 )Cn)« 

TF(.I«1) THFN T 4 : c-Ai pua: r T 1 • CMAT . 




- 39 - 


• • ft * 


TF(.J«.« ) THFN T5i*-ALPHAC>)Cn«FNDn 

~ 1 F ( R < J..«^ ivK ) C I KR t . K ) f n 1* THFN 

— SSC 1 2 >« — RFT A n > C n t 


-FI RF 

— FiRn ? 3 »' 
-FNOt » 


•PFTA( ivri;> 


-ArMr3U-AAMAV( nttii 


— AnJC41t*nAMAV(2)CIH - 
— -An.ir5H«AAMAWM » r Mt- 


-CAI.L— TRANSFFRt- 


-FNDFI 


-FNDPi 


OIP ? - 

0\9^ 

01^>4 


PRfiCFDI IRE RNPOl.H 

-V#CC— — PPOr.FDURF- ''SNPiXT*'- riFAl.S-UITM THORF FOl.lATTONS WHICH -ARE CCCC 01^7 

4: SlITTARI.E-FnR -SOI ITH AND NORTH POl FR- TROPOPAI.IRFS. T T -FST IMATFS C 0» V«, -■ 

-C rOFFFtCIENTS-tTF EACH CORRFSPONOING! 'M.AWnA"TeRHS AND TRANSFER:? C 0197 

-CCCC THFM TO THF -WORkTINO HATR TCRS-"T ALPHA" , "SS"-^ CC*/ 019R 

IF_,>V THFN FNOII ; Ot9S>- 

IF THFN ..Ml ■.« -2 FL; •/ 1 0200 

Al ,PHA( I ) n * Tl.T I -AI.A/ ( RO*RO*YY ) I 


.FU=:jrA( I > M i4L5«-2.0*AI- »*AL 1 /ROt- 


-OV t * ■ 1 . 0 / ( 2 . 0* Al S* Al ,0*R0»Y Y ) « 
t1_.T0 -ii^no 


«+Al.PHA( I ) Cl 31 


.4^0R in 
Tin 

^IFUI«^U-THFN,,TAI «.-AI PHA(4 ) f T llFNCiIl— 

IF < .!« - THEN. T5 z «-AL PHA < U-C T1 X FNDI X- 

RSt-1 2.1 1 -+BFTA ( M-t LI l 

Anjcsi X *0AMAV( 1 X c nx 


.0201 


CALI TRANSFERl. 


—FNDFX — 


J^NOPX : 

PROnFnilRR INTFPSl.— 



-COCO «’> 2 1 0 — 

r. 0211- 

0212— 

„G 0217— 

OF SHRFACF -POINTS, THF C 021 A— 


_4»CC"---.^PROr.Fril IRF .‘iTNTFPS"— ORAL:? WITH -THOSE FPI IAT I0N:7- UHTCH ARF 

—C SUITARl.F FOR-,Al.L -INTERIOR POINTS-ANO SURFArF POINTS FXCFPT 

_C SOUTH AND NORTH POI..RS SURFACRS.-.IT - ESTIMATES r.ORFFICTFNTS -OF. 

_ C FAOH CORRESPONOTNO -"-LAMDA!*-TFRMS- AND TRANSFFRS-T+IFM -^TO- THF 

_C WORKING MATRir.FS^".TAl.PHA",.'iSS"... IN -CASE 

_C:r.CC-.-.-FIRST-.'lIF."_STATFMFNT ACCESS- TO ANOTHER PROC-Fni.lRF !‘1NSUR'1. 02L*;-. 

I;-,2: «R C.J+ 1 . Kl.) Cl ML I-R < ..H- 1 . U C 1 ML-Jx 02 1-7- 

RTRSlSIf 1 l-IL3:«Rt,l-l .KL in : U-3-R(..UlM U I Mlrjt 0210- 

RIRS ( A ) C 1 M L T L3 X «RTR:=5 ( 1 ) C 2 1 1 1 02 19- 

RTRSlAiCTl ll«RTRS(-l ICMl 0220- 


.RTRS(.«!U2t U.IX»RTRS( l-I-C 1 X U.-IAX- 

.RTRS151C 1 llsRTRSl 1 1 C Il.lX — 

.RRSr.-U:»R(.J.K) C2sl-R( J.K) CTLIX — 




— RR.stiLit «*R( j.K 1 r 1-1-Ri J.K1 r ii.-nx -- 


-02 


-A- 


_RR3t2xlL-13J-^R(.J.K» C.Sx U 1-R( .J.K > C-1 1 II -21x 

_BRAriML3:«RI.Hl.Kir 1: II. 3-R( J-1 .K> 1 1 Ml. 3 X 

,.RR5r l.lx«:R( J.K> C23r-2.0#R(.I,K)C1-1 + R(.J.K)CIL1X 

_RRr.LIL3x«RlJ.iail3-2.0*RlJ.Kir-U.3+R( J.K>r IL -1 31 

RR.SC2* U.-l3X«R(J.K)f-Sx U.1-2.0*R<.I.KU 2t II -1 3.*.R<.l.k ir 


X j Xl -^-Xt 


- 4 !»? 27 — H 


1 






-I|0- 


RR6tli n 5«»R(J^t,K)nUL3-?,0*R( J.K)f II IL3+RUI-1 .KJCI t II )l 
-AMI (Al;t/(AL^i*TSt > ( <RR3CH U.1)/<2.0*XX#* 

•RUUKVn I n l*TR1*CnS(TH) ) »♦♦?»♦<•< ( <RR4f n n.))7(?,0*VY*- 
R<J.K>Cli TI.3*TSU l♦♦:i>♦A^.4)^ -- 

lF-(vKvH > .AND. (K«-l )-THFN CAU -INSHRl ENDII - 


0>91 


IF“(J<JI ) .AND.-<K»1)-THFN 00 TO LA » FNDIt 

-Rt 1 1 ILll-MALl /Al ?)**?)*?.0/IR(.I.K»r 1 1 n-l*TSl )•♦•( (RR^f I* IL3 

-♦(TS4-TS5)-)/<2.D*( (R<..I,K> C 1 1 IL3*TS\ -02?!?S -* 

— *COR(TH)*XX)**?> > »-(RR5r 11 H. 1 / < TSHm RO*XX#COS< TH» > **2 

—)>>♦( ( <AI.4*RR4Cli tl.l)*<TS2-TS3) )/<2.0»( <RO*TSl 0 — 

♦YY>**2> > »-l (Al 4*RRA[ 1 1 IL3 »/(TSl*( ( O ' 

-R0*YV)#*2) )) + < <TAN(TH)#(AI.4*RR4Cli IL1) )/<2.0*YY^»Tfii*RO*RO) )1 0 * 

-CCll Il.ll— CRR3C 1 1 II 1 )/(XX*TSliK <RO*rOS(TH)»**2) M 0740 

- FCli TLJI— <AL4*RR4Cli TLl)/(YY*TSl*Rn»RO)t 0 — 


-ALFHAI 1 »rn U-3i-7fO»l ( 1 1 0/ ( (RO#COS< TH) *XX )**7» )♦( AL 4/ 1 (RO 

*YV»**2)>-t'(72*ACliin/(2.0*Z2)))1 - ‘ 

Al PHA(7mt IL3:«l.0/( (RO*C05ilTH>*XX>**2M — - -- 0744 

-ALPHA(3)C1» ILlt-( < l.O/(RO*RO*YY*YY) )-(TAN<TH)/<2.0*R0*R0*YY) ) — 

-Al PHA(4>C-M Tl 3»-( ( 1.0/-(R0*R0*YY*YY) > + (TANITH)/(7.0*R0 *rIi*YYM 


— AlPHA(5»tUII 3i-(71*Aril H 3/(77) )-MPr 1 1 II. 3/ 17.0*77 Ml • - 

— AI.PHA(/>)ClML3t«(7r^*ACtt IL1/C77) )-(BClt IL3/(7.0»ZZ) )l • 

— Al PHA(7)f )lIl.3i-CCl I Il.3/(4.0*77*XX)I =- — 0749 

— ALPHA(8)Cli im-f-Cli IL1/(4.0*Z7*YY)I 07!50 - 

DFTA(-I >r II II Jt-AI 1*Al.l /(Rn*CO«ITH)*XX)H 07S1 

— BFTA (-2 ) r-l i-Il. 1 1 «AL 1 ♦Al. 1 / ( RO*YY H * 0257 

— PFTA(3)f II II 31-AI. 1#AU /(TS1*77)I 0253 — 

-..BFTA<4>i:ii-IL3i-ALt*ALl#(RR3C1 1 TL1>/(2.0*XX*ZZ*TSI*RO#COS(TH) )l " 0254 -- - 

RFTA(5)C 1 1 IL)l-Al.1*AU *(RR4r I t IL3 ) / ( 2. 0*YY*Z7*RO«TP.l ) I 0255 

— BFTA(MCli II 1i»7.0*Al.t*AI-l*TAN(TH)/R0l 025<S 

_—PFTA(7) t It II 3i»4.n*AU *01. 1 /ROl : 0757 

r.l J 1 4 « 1 . 0 / ( 4 . 0* Al A ♦AL I ♦R0*C08 ( TH ) *X X > I 

r.ll7i-RR3ril II 3/IP.0*AU*Al 1 ♦RO*TSl ♦fOS( TH)*XX*77 ) 1 "" 

OVl X »l . 0/ ( 4. 0*AI..3*AI.3*R0*YY) 


T.O?l «RR4rj HI. 3 / ( P. 0*AI 3*AL3*R0*TP1 ♦YY*77 ) 1 - 

QU 1 X ■ 1 . O/ ( 4*. 0*AL2* AL 2* T3 1 ♦ Z 7 > I 

f OR-4 X - i—To -I u no 

-Lt X ■ T ♦ t4 ■ ■ ' ... ■ ■ . — 

IF 1 1 -Tl.+ l -THFN .1 1 X »1 X FNni X 

I F— I ?»0 - THFN 1 7 X ■ 1 1 -I FMD I X 

T 4 X «4-alpha ( 1 ) rn X 

.T7X —ALPHA! 7) IT 31 

J3x «sAl PHA < 7 1 1 13 X 

T4X—AI PHA(3)C IlX 

JSI —ALPHA ( 4 » [ n X 

I6X— ALPHA(5)C 111 

^T7t— Al PHA(f.»tnt 

^TSx— ALPHA(7)CI3X 

T9X-+AI PHA(7ir Tlx 


0760 - 

07 A 1 — 

- 

07 6 

0?64 . 

— — — - 

^ 

0??4 

0?7A 

o:>76 . 

0?7j? 




T10:»+ALPHA(7)CIlt 



02S4 

. TM t ( 7 ) r 1 V t‘ 



07 

-1 r HurnH i-j i 

-1 t rif - 




—TV4s =+ALPHA( 8) L i J ! 

-T1 5? =-ALPHA ( B > M 3 1 

— 

— 

... r 

02 - 


-- ssrtu»+BETA(t>i:ni 1 

— G.S(?]t=-BFTA(l)rr]ll ■ 

— ssi:!=:i t »+bi-:ta<-^2 vc i 3 1 — — — 

— SSf M i «-PFTA ( 2 > M 

— ^BSC-93 ; =+B€TA ( 3 ) C IU 

— SStl03t =-BFTA(3)r I3» ' 

— ^SSC33:=-BFTA<4>m? 

— SSE4 3!=+BFTA(4> ri3? 

— SSC7-3 : «-BF ' 3 5 ) C T-3 

— SStA3:*+BF; ■ >5>rT3{ — — 

— S;rir-121 : =-BP.3 A < 6. / !! 1 3 ; 

— SBr 1 33 : =+BFTA ( 7 ) T-1 3 ? 

ADOt-l 3 : »OAMALI ( I ) C T 1 

An.JC23 » «GAHAI.I( 2) r 13 « - 

AB. IC43-1 : *GAMAV < -1-) l-I-3 y 

ADJt 4 3 : »GAMAy ( 2 > C 1 3 ; - 

AB.JC5-3 : =GAMAi4 < 1 > C 1-3 ? - 

GAl L— TRANSFER? 


^MBF-?- 0307 

4,A ENDP ? i 

RROCFriURF— -I NSl IR ' > j 


_/*CG— FROCFfcilJRR — N3UR " DEALS W I TH -THOSE EQUAT T ONS WH T OH- ARR CCCC 03 10 ' 

.X SLI I TAPLF-FOR- ALL- SURFACF-PO I NTS-EXCLUO I NO - SOUTH- AND NORTH 03 14 

_C POl.F-S -ONE'-S,.— T-T-EST T MATES -GOFFF T C I ENTS -OF- EACH -CORRESPONDING C 03 1-2 

_C ?il AMDAII-- TERMS- AND -TRANSFERS -THEM-TO 4HE-W0RK; ) NO -M ATR I CES G 03 1 3 

_C =i'-TAl.PHA-"->-:iSSi:,: -♦/ 

ALPHA LW C-1-: IL-3 : =APS( RR3C-l-t-IL-3 / (-2-^0itXX« XX«RO*RO*- 

COS ( TH > *COS ( TH ) ) ) ? 

ALPHA L2-) f l : 1 1. 3 : =ABS ( RR4 f U-I L 3 *AL4/ ( 2. 0*R0«-R0*W*VY )-) ? 

ALPHA ( 3 y C 1-: I? .-3 : =A C 1 : 1 L 1 *TS t 

PF.-TA < 1 4 f -I s IL!I s =ALJ *AL 1 «-RR3r 1 : Il=-a-74R0*C0S( TH>-*XX >? 

BFTA (-S-J U s I1..1 : =AL1*AL1*RR4C4 : IL34-< R0*VY->-? 



PETAiSi r-l-? H..-3t=2.0*AU *A14 t 

I3U 1 :-= 1 . 0 X- ( 4 , 0* AL I *AL 1 *RO*TS 1 4C0S ( TH ) *X X > ? 

G 1 .I 2 -? =RR3r i-t.U.3*71 /(4.0*Al.l *At m-RO«TSl *COS( TH ) *XX )+ 

OVl ; =1 .-0/44. 0*AL-;3»AL3»R0»YY-> ?— 

GV2 : *RR4 r 1 : 1 1. 3 *7 1 / ( 4 , 0*AI 3*Al 3#R0*TS 1 *Y Y-)-? 

GUI ? »Z-t/4 2, 0«AL-2*AI.2*TSl ) ? 

FCiR-I :-=-1_T0-U— no 

I-l:=I + 1? 

^I2!=T--U 

IF-I-t = IL+t -THEN I I : =l ?-ENDI ? 

^IF— 12=0 -THFN - I2s =IL? - F-NDI 

I4!.-ALPHAC4-) C-I 3 +A1..RHA (-2-)-C I J+AI.PHA ( 34 C4-3 ? 

. . TA:=-AI PHA(3irTlt , . . 


ORIGINAL PAtifi il.' 
Of' PC/OR QUALITY 


'i 


-0319'- 
-0320 • 


-0321 ■] 

.0:1124 j 



ssn n s “-BETA n u n t 

SS r l ? I t “-BET A ( ? ) r i-i » 

BSC 133 1“BETA ( 3> C U ^ 

— IFHR(-clrl )Cn 3>“R< J» 1 > f I?3 )tAND. (R( J+l * 1 )Cn “ 

^>R<*J-lrl)Cn) > THEN 

T3 1— ALPHA H-ye-int 

T5i— ALPHA(?> Cnt 

—END! I 

— IFdRUI. 1 ) r I I JC«R(c»«-l > Cr;>3-K AND. ( R( 1-3 — 

CR(r.l-l rl > C I 3 ) ) -THEN- 

_T2! —ALPHA ( 1 ) C T 3 t 

-T4 1 “-ALPHA ( 2 ) C I 3 J 

ENDH 

-I-F-(-<R(J, 1 )C ri3<R(Jvl )ri23».AND. (R(J+1, 1 )CI3 

^>“R( J-1 , 1 ) CT 3 ) )--THEN 

T2:“-ALPHA(l)CIlt 

-TF.:=-Al PHA(2)ri3t 

PNOI ; V 

^IF.( (R( J,-l >ri1 3>R(Jrl )ri23)*AND.-(R(..l+l, 1 )M3 

<“R(.J-V, 1 >CI3> ) -THEN . 

T3:»-ALPHA( L)CI3t 

^T4; =-alPHA ( 2) r I 3 5 

ENBI? 

AD.IC 1 3 ;“GAMAI.H J ) C I 3t 

AD..IC;^-3 : “GAM AIK 2) C l 1 % ^ 

^An..ir-33 : “GAMAV ( 1 > C I 3-t 

ATiJC 4 3 s “0AM AV; ( 2 ) C I 3 ? 

AD. Iff. 3 : “OAMAW (1 ) f I 3! 

CALL— TRANSFER! 

FNOF! ^ 

ENDPt 

FROCFniJPE-TROPO I ! 

Vi*CC.- — PROCEDURE — “-TROPOi'i. j:i|rA|..s. wxTH THOSE EQUATIONS -WHICH ARE 

C SUITABIrF FOR •TROPOPAUSE INTERIOR POINTS. IT ESTIMATES 

C COEFF-ICIENTS-*OF-FACH CORRESPONDING "LAHDA"-TERMS AND 

CfXtC— TRANSFERS -THEM- TO THE WORKING- MATRICES-"TAl.PHA" , “SS" .- 

AUTOMAT I C -REAI -UEC-TOR C 1 1. 3 RR 1 0 , RR 1 1 ! 

RR 1 0 r U t “R ( U » K ). r-2 3 -R ( J , K ) r 1 1 -3 ? 

RR to C 14. -3 : =R ( . J , K > C n -R ( . J . K ) C I L- 1-3 ! 

— RR10C2:-Il.--l-3: =R(...I»K) rs: Il-3-R(...LK> C 1 1 II. -23! 

RR 1 1 r,l : .1 1 . 1 ! =R ( . J+ 1 , K ) C 1 s I L 3 -R ( .J-4 » K ) C-t : T L 1 ! 

ALPHAU )Cl: IU3t“l . 0/ ( ( RO«-r:OS ( TH) « XX ) <h(-2 ) ! 

ALPHA ( 2 U-1 ! I L-3 ! = (( 1 . 0- < TAN ( TH ) *YY/2 . 0 > ) / (-( RO^^YY ) > ) *AI.4 ! 

^ ALPHA(2) r.l s IL 3I=( ( 1-.-0+ ( TAN (TH > *YY/2 . 0 )) / ( ( ROfrYY ) <Hi-2 )-) «-AL 4! 

A C t: I L 3 : = (-( AL \ / ( AL2*TS I > ) **2 ) + (- < ( RR 1 0 C4 : I L 3 ) / ( 2 . 0* X X *R0 

*TSt*COS(TH) ).)«-*2» + ( ( ( (RRl 1 C 1 : IL 3 ) / ( 2, 0*YYfrR0«-TSl ) > **2 > «-AI4 ) ! 

ALPHA ( 4 ) C-i-i I L 3 : “+ C2 . 0 / ( ( RO*COS ( TH ) *X X ) ) > + ( 2 . 0»AI. 4 / ( ( RO*YY )- 

** 7 .) .) + L72*A r. J-! 1 1 . 3 / ( 7 7 ) ) s 

ALPHA ( .E. > C t ! 1 1. 1 : “ 7 2 * A Cl : I L 3 / ( 7 7 ) ! 

BETA(-1 ) El : IL3i“Al4 «Al t/ (RO*COS(TH)*XX > ! — : 

. PFTA(2iriXll 3!=AI 1*AI 1/<R0*YYIt 


I 


v*« ♦. • « 


‘ BFTAOKU n 3J-AI l*AU/(TSl*7?>l 

AHTA<4>rtin 3l-AI.1#Al.l*(RR10CttTl.1>/(2.0»XX*77»Tfil»R0*cnS(TM) )l ■ * ~ 

PFTA ( 71 ) tl f I », 3 1 - Al. 1 *Al 1 * < RR 1 1 n I 1 1 3 ) / ( ;■* ,0« YY« 7 7 »T$ 1« RO M 0S:^4 - 

PFTA(-SK1,1I 3l»2.0*AI.1*AI t*TAN(TH) /ROl OS'SSS • 

FlFTAI7>m n.ll-4.0*AUfrAI 1 /HOt - - 

r,l u I - 1 . 0/ ( 4 . 0<» At. t *AI 1 #Rn*rOS ( TH ) X M 

OVll-l .0/(4,0*^Al.;'^*Al 3#R0*YYM - - 

FOR T«-1 VO n no - - 

l U-T-t-n 07. v; - 

IP n«n.+ t THFN Tt«-UFNnn - - -053?- 

; IP i^«o thfn i?t«iMFNnn osao - 

— • — Tll"AI PHA(4>m« - “ • V'54.‘ - 

; T?«— A» PHAM >n 1« - 05 A4 

j T3 i i.-AI,PHA( t K 1 11 0! .4^- 

L T4«»-Al PHA(?)rnt - - 0540 

I T5»— Al PHA(3>tni — — - , - OSSvl 

T7j—Al.PHA(SVrn« '• OMr..v 

:>3rtl:— ►BFTAUUni 

S5t ?U«-PFTA( Ur Ml - - 

Aors 1 1 «-t-HFTA < ,■> ) r 1 1 % 

— 33rM»«-HrTA(?Ktl< - ^ - - -- 

i 33C‘?l:«»BFTA(3Un» — — 

ssr ioi:=-F«FTA(3>rn» 

.S3[3U=-BF.TA(4Kn; : • -• 

?.3i 7ij»-nFTA(f.> r m; — — — .. 

33r41:=i-HFTA<4)r ns 

scroll =+nrTA(5) r MS - 

SSCl?lsa-BFrA<.->>rMS — - 


131s et + OFTA ( V ) t M s 

An.lC 1 1 : -OAMAIK DM 1 S 

ADJl 31 s =OAMAV ( 1-H 1 3 1 


CALL -TRAN^FF.Rs 

PROcrnnRF tramsffrs . - - 

/.* THIO PROl FnilRr TRANsFFR f-ALFl « ATF O VAl UF.3 CtF 

FI.FNFN V3. OF-. COFFF 1C 1 FNT NA3 R I X ( Ci.'ir NTRX ) . 501 .LIT T ON 

MAVRIX (SOI MTRX) v-AOvinSTlNO MATRIX . (An.lMTRX K-TO 

-THR CORr.SPiTNOlMO FILF:S NANFO AS rOFHTRX . SOI N VRX - 

APJMIRX. -THFSF- FILLS WII I -Kl' Ai.r.FSSFTt IN WINOAtUUST 

COUNTS -COMNT+1 S . — - 

TX (t.uUNI ) 1 1 ' 151s «TAI.PHAS 

SX * f.OI INT-) LI ' 1 ,31 s «SSS 

ADJX(CnUNT)CI 'STs^AUU- 

TAl PHAT 1 S J53S =0.0t - - 

AOJCISf.ls^O.Ot 

' SSt IS I33ssi0.0s.._ 

; ENDPt 

1 F.NOMt 


0SS4-- 


FILE LAMDA 



— nrt •-> 9mrt» ^ 

*RFWTNn,LAMOA. 

- /rr>f»Y? hf-* 1 ■ — — -— - — 

1 Monui, F W I NriAnjS ( UN I T l -COFMTR X » I in I T?=F0I MTRX . UN r T3= 

— An..lMTRX,UNIT/S=RF:=.Ul.T)- RFI'TN MATN; — 

=/»-THIS FRDriRAM AH. JUSTS OBRFRVFD WIND FIKI P TO MAKF 

ITL MASS-COMSFRVATTVF wind FTFJ.n. AP.IU3TING- TFCHNIQIJF — 

IS PASFH ON OARITTONAI, FORMUI ATION. SIMI U TANFOIIS 

FQUATIONS- FOR I.AMDA -TFRHS ARF--SOI VFO IJSTNO -SI ICCFSStOF 

OOFR RFIAXATION (S-O-R) ITFRATTON HFTHOO. FIl.FS NANFO AS 

r.OFMTRX . SOI MTRX . AND AD.JMTRX — UAVF- NFSSASARY 

rSTJ MATFn f.OFFFrC.IFNTS^ -THFSF-FIl.FS-WII.l.-.RF-RFAn -I.N....rMIS^ 

— F'RiTlTRAM. #/ . - — — . . 

/.* -THE-FOU.niJING . SFCTION NFFDS TO -PE.-CHANOF.n .FCiR . niFFFRFNT 

GRID -SYSTEM. _XTRAt-,TO XTRA7-AND TSM1 TO-TSM7- ARE 

RASINfi PARAMFTFRS. -AL 1 -IS -WF IGHT ING FACTOR -FOR -U -OFl OCITY- 

C.OMPONF.NT */. 

/♦...PF.GlNrNG.nF-THF-SFrTrCiN «•/ 

1 . T TFR Al I , Y S . 0 'i- 0 YY ? 

I-ITFRAILY- "T?." 11. ♦. 

L T TRRAU.Y '-l ■ Jl. : 

L-IITFRAII.Y 

l.ITF;RAI.I.Y-l'\0''kl.U.; /* -KUJ.= KI. + 1 -*/ 

1 ITFRAi I V- .1'^S•?7^•• UK? 

U T TER ALLY - " S3:"<7(S ", -X TR A 1 ? 

LI TFRAL LY . .'U 5992." X.TRA2 ? 

LTTFRALLY ‘'17SS4'' XTRAS? . 


I 


-i,5- 


LITFRAI.IY ''7:4512" XTRA4I 

It T TRR ALI ; V - " 5 n 5 ? '• -X TR A5 1 

ItJ-TFRAU V-*‘ H 1 9?“ -XTRA7I 

L-. TTRRAl hY ••7A8" ISMH - 

UTFRAl LY ''6.40” tRM?t 

1.-T TFRAI.LY-" 1 ?8 - T SM3I 

LITFRAI.I Y’ "?56" -ISM4I 

IJTFRAU.Y "256" ISM5t 

ItI TFRA l-l-Y ?56 1 FM7t 

LITFRAI l.Y "0.01"Al.ll 

1 IXFRaU.Y "0. 001" -TOI .RNCFI 

ItI TRR At. I . V - *• 1 . 75 "RL XF A^ t 

I.ITERAl.l.Y "3. 14159?65;‘-!SR'=>793"Pn 

/* -IL«NI.IMBFR OF ARTOS IN FAST-WRST DIRRCTTON. 

^ xll.»NIIMPFR OF ORinS IN NORTH-SOUTH HIRFCTION. 

KT “NMMBFR OF ORTAS IN Al.TITUDR BTREGTION. 

ijkbtotai. numbfr of grids. - 

toi.rncr«toi.rancr in ITFRATIONS. - - 

RLXFAC=RFl AXATION FACTOR IN SUiX.FSTVR OVFR RFL AXTTON 

METHOD. ♦/ 

V.*-4^Nn-0F--THR- SECT TON */ 

RFAl VECTOR 1155 ARRAY( UK) TX I ’ 

COMMON <5T ( TX . TDI IMPl M 

RFAL VFCTOR rXTRAll TDUMPlI 

RF.AI vector C5V ARRAY (UK) AD.JXl 

COMMON CA(An.,IX,Tni.lMP3)» ^ 

REAL - VECTOR CXTRA71 TDI.IMP3: 

RFAl VECTOR r T.IKl ARRAY ( 13» SSXl 

COMMON C.SS<SSX.TDL1MP4)I 

RFAl -VECTOR r XTRA25 -TrillMP4! 

RFAL VECTOR- C 151 ARRAY! UK) TLAMDAl 

— COMMON Bl (TI.AMnA.TDUMPtl » ; 

RFAl, VECTOR- C XTRA1 1 -TOUMPl I ; 

RFAl VFf.TiSR riJKl S.KIVl 

j COMMON fi2(S.DW,TnUMP5) ! 

i RFAl -VECTOR TXTRA31 TDUMP51 

RF Al. -VECTOR C51 AD..I; — 

RFAl VFCIOR nil ARRAY (JL.KI.) ALAMDA , PLAMPA; — 

COMMON - B3 ( PI AMDA , Al . AMDA . TOI IMP6 > I — 

I real VECTOR TXTRA33 TOUMPA? 

j RFAL- VECTOR TIL] -ARRAVUL.Kl LI )-U. V. W. P, APJl ISTIU ADUUSl V, 

AO..IUSTW? 

COMMON R4 ( I I , V , W . n , TDMMP7 ) ; 

, REAL VECTOR T XTRA4 1 . J PUMP71 

' COMMON P5( An.tl.ISTI I, An.lUSTV, ADUUSTW, TOUMPO) ; 

; RFAl VECTOR I XTRA53 TOlIMPSt 

, REAL. VECTOR CIJKl ARRAY ( 13 )...VEL; 

• COMMON PA I VEL . TPI IMP’? ) ; 

i REAL VECTOR rXTRA?'l T0UMP9; 

s REAL VECTOR Ckl.l.l 3 7? 

, TNT TTOI ( 7=\o. n. n ri^ . a 


*. 


♦ 0# I I 

RFP\i. Al:3 . T . FRROR * 7 Z 7 • THFTA TH , CON , MA XP t V . M I NO! V t 

INTFOFR CHFC.Ki t , . K. Nlvrl 1 rT;?. 13. T4, .J1 , .J'S, K2. 

K3 . r TFRAT I ON . K J i T W , F I . F? . F3 . F4 . FS , COUNT I NO i 

E6 . E7 , F3 , E9 , F 1 0 . NMMM . . .U » 

— REAL VFCTOR- r3331 Wl AMDA. Wl AhinAH 

PFAI. VECTOR [93 KLAMOAl 

RFAU MAXI AMDA.MINl.AMDAJ 

PPOCFOI IRE MAIN! 

OPFNMAP( ■fFl '-vTl AMOA( 1 K t J. FRMl .FI . 'I.AROF'’ > » 

OPFNM AP ( • F2 , S C n . 1 3M3 , E 3 . ‘ I. AROE ‘ > » 

OPFNMAP( '■F3'‘,PI AMnA( 1 , I ) f I 1 . T3M3. F3. 'I.AROF')! 

— OPFNM AP ( ' F4 , U ( 1 , I ) [ n , ISM4 . E4 , 'LAROE' ) t 

OPFNMAP( 'FR'’*. AOJIISTMC 1 . 1 U 1 3. ISM5.F5. AROF'' ) « 

-OPFNMAPF- VFL< 1 > C 1 1 . TSM3,F-^, 'LAROE I 

OPFMMAP( 'COFMTRX", TX( 1 >[13. TSM1 ,FP. 'I AROF' M •. 

OPENMAP ( 'SOL MTRX ' , S3X < I ) C t 3 . TSM2, E9, 'LAROE' ) 1 

-OPFNMAP( < AO..IMTRX . AOJX ( t) Cl 3 . FSM7, FJ 0. ^I.AROF'-) X 

-WRITE UNI T ( O ) ( F I , F3 . E3. E4. F5. FO. E3, E9 — 

,FJ0>4''CiPFN'', -9<3Xr3>#; . 



COUNT T NO : =0 X ^ 

WRITF LINTT-(6) (DYY. |T, , Jl , KL. AU , TOLRNCF. Rl XFAC )» 'nYY='' , F5, 2, - 

'IL=:’ . 13, 13. •'KT.=-'-.-I3, ' Al..t = -',F7.3,.''T0LRNi:;E=--.F3.4,- 

:irRl.,XFACT0R=<..FA.34? 

430N:=PI/-tOO.O? 

k'-?! *kl . + l X 

/-«=_ INTT-IAl T7TN0 ALAMPA ANP PLAMPA-ARRAYS-OTVING 

F T ROT 01.IE33 VAI .IJES-*/ 

-THFTA; =+90. 0+PYY; 




-1 


j 

i 



_FOR ..I:=:t TO-^.IL-DO 

-THFTA; =thfta-pyy; 

..for k:=l jo kl no 

PI AHDAt.i.kU i; ri 3;=?.o;— 

ALAMDA<..i,l<) Cl ; Il.3;=3.0;— - 

.F.NPF; 


FOR -J; = i-..TD Ji no 

FOR K : =1 _TQ _KL DO 

-FOR-J-; = i-iHY-:r;_TO-iL._nQ 

—J^l . AMD A ( J . k ) r 1-3 i =» 1 . 0 S 

_Al.AMPA( J.k) C I 3;=l .0: 

-ENDF! FNOF; FNPF ; */ 

-.PL AM DA ( r. . f. ) r 3/S 1 ; = 1 . 5 ; 

. Al.AMDA(5.S) C3/Sl;=l ,5? 

FOR -J;=) -TO-.il. PO 

-Wl.AMPACJl: =01 AMnA(..L5) [363; 

ENPF; 

WRITE LINITIt) (Wl AMOAf l; Jl 3 )4SE14.7#; - 

FOR K;=l TO KL DO ■ 

-KLAMDA [ K.li=BL AMD A C5...1C1 C 36 3.; i 

FNOF; 



-Cl Oi>F( ■ 


- Cl.n;>F( "F3-,F3)» 

—Cl OSF ( •• F 4 , F 4 > ; 

-CKOCFl 'FS '-,F5) ; 

— Cl nSFl F/S'-.F^^M 

- Cl OSF ( CCFMTRX ' , FCO ? 

Cl OSF ( ' RCil MTRX , F9 ) 5 

—Cl OSF ( AHv IMTR X - - F 1 0 ) : 

WRITF -UNTT(6.) ( FI , F? , F3. F4. F5, F6. FR 

,F9,F10)»- CLOCF'i -9(3X12)11: - 

FNDP : 

« 


PROCFnilRF n I VRCNCF : - — 

_/*._TMTR PRnCFnURF-FRTTMATFC THF- VAl UF of RIGHT HANP RinF 

OF FOUATTONR.. AT THF FND IT CAI.CIJI.ATF3--THF-DIVF.ROFNCF 

OE-DBRFRVFn-nR An^iuRTFn wiNn FiFi n */ 

FOR J: = t-TO Jl. no — 

FOR K:=sl-JrO KL-DO 

FOR I:>-1-T0 II - no 

X2:=>J-l; — 

IF Il = Il.+ l THRM Il:=i:fcNni: — 

JF-12=0 THFIM I2:=Il.?ENDl; — 

vn.n > CNL] ! =i.i( j.K) r i n»n( j,k> c in ; 

VFI.(2)[NL3;=U(.I,K)CI21*D(J.K)CI2D: 

VELI3) IHLIl =U(J.K+I ) f-I 3*P( vU k + 1 ) 1 1 3 < 




IF ka1 THFN 


vFL( 4 )rNi.i«-(?*u(j, 1 )rn-u(.,i.2)rn>« (;'«-n(wU i > m 


n<ij*2>cn>» 


FOR- .J: = t -TO- ►IL— DO 
POR-K: = l_TO-Klr-nCi 
FOR- l: = l— TO-TL--DO 


VFL ( 9 ) r NL 1 : =W (a. 1^ K+1 ) C I 3*0 ( J. K+l -VC I -3 ; 

IF--K>1 -then 

-VFl . (1 0 > r NL 1 : = ( ( J , 1 ) C IJ -W « 2 » C-I 3 ) * ( ( J , l-> C T-3 

J. D<.J.2)ri-3>» 


VFl-( 10) CNl 3»=W( J,K-I ) r I-3frDUUK-'l-) 11-3 


VFl < 1 J ) CNI.3 i =U(0-K> C T 3 «-D ( K >^C- 1-3 
VFM-1 2 ) r NL 3 : =V ( .J . K > C T 3 *D ( vU K ) C I 3 
VFl.-( 1 3 ) r NL-3 s =W UU K > r-I 3 *0 UU K > r-i-3 


y-*_CAl.C;UL.AT I NO -R J OHT-HAND STDE OF-EQUATT i.iNS -«•/ 

SCl.MJK3:=0.0*, — 

FOR-i s =t_Tn-i 3- no 

SU LT Jk'3: =sr 1 : I JK3+33X ( I-) C U IJK3*VFI ( I-> 1 1 1 MK-3 

ENDF 

/-»-r.ALr:UI. ATTNO niVFROENr.F-OF UINP FIFI D -* / 

AL3! =-r-2, 0»AL 1 *AI. 1 1 

ni.VC.U.I JK3 : =SL1 1 IJK.3-4AL3: 

MAXniV:*MAX(DIV) ; 







M!NOTVl«MIN(niV)l 

— WRTTE-IJNI‘T(6)(MAXPIV.HTNDTVH»^MAXDTV«^»F14.7« 

•‘MINDIV«^>ei4.-7#» ^ 



PROCEDURE ROR» — 

/*.THIR PROCFnt.lRF--IS A-PROORAM-FOR-A--SUCCFSSIVE- OVER 

RELAXATION (S-O-R) ITERATION METHOD TO SOLVE SIMULTANEOUS 

FQMATIONS. FIRST- VELOCITIES-AND LAMDA TERMS-ARE ARRANGED 

IN VECTOR FORMS.— AND-THEN -S-O-R -EQUATION -IN- GENERAL- 

FORM -IS USFD- FOR -ITERATIONS.—*/ 



LDI NL« 

/♦-INITIALIZING TLAMOA ARRAYS -*/ 

FOR -I I » 1 -TO I . »K- 00 — 

^Tl.AMDA ( H r-l I 1 l «0. Ol ^ 



CHFCKt - t t — 

FOR- Ji - 1 TO -tJU-00 h 

FOR K J » I -TO - Kl.— no 

FOR— I «»1 -TO- Ib-DO 

— : NI.-i-»NLr*l I 

Ill - I + l I 





1 4 1 » I -21 ^ 

— Jl : ».J+ 1 

U2J »J-l » — — 

Kl : »K+l » ; 



1 F T l » I L-H- THEN -M i » 1 » END 1 1 

^IFi_T2»0 -THEN I2i»TLlEN0I1 

-IF -U1 «..l| -THFN-JI t «.JI -1 1 ENDI % 

IF-J2*0 THEN .J2« “2» ENOI I 

L_l F- K 1 =KL.+ 1 •— THEN K 1 1 -Kl. 1 END 1 1 

^ ^IF- K3=a~THEN -K3» »l s ENOI 

F 1 3= 1 1. +2 THF.N 1 3 « -2 » END I 

I F-T 3- T L+ 1— THEN -1 3 1 = I S END 

IF—1 4=0 THEN 1 4 1 » T I .« FNO I » 


_4F-l4=.^-l THEN -1 4* =11.-4-; ENDI ; 

^TLAMOA ( Nl. > r I 3 1 »Pl AMOA ( K > I T 3 1 

Tl. AMDA ( Nl . ) C 2 1 ; =01 .AMOA < .J . K ) C II 3 1 — 

-TL AMD A ( Nl.) r 3 3 » =PI. AMOA ( J , K > C 1 2 3 X 

TLAMOA ( NL ) i: 4 3 1 =BL AMOA ( ..U , K ) Cl 3 X 

TLAMOA ( NL ) C 5 3 1 =PI AMOA ( vl? . K ) C I 3 ; 

TLAMOA ( NL ) C A 3 t «BL AMOA ( . J , K I > 1 13 1 

TL AMD A ( Nl, ) r 7 3 » «Bl AMDA ( J . K3 ) C I 3 1 — 

ILAMDA(NL) C83i =BLAMDA( J.Kl ) Cl 14»— 

^TLAMDA ( NL ) r 9 3 » »BI. AMDA I J . K3 > C II 3 «— 

TLAM0A(NI.)C103s=BI AMDA(J.Kl)CI23» 

TL AMDA ( NL > C 11 3 l =PI ..AMD A C J. K3 ) C 1 2 3 1 


Tl AMnftfwi 1 r t-jitasW .fc'i \mt . 


Tl.AHDA ( NL ) n 3 ] t «BU AHDA ( J 1 « K3 ) C T 1 « 

ThAHOAf NL-) Cl 4 3 1 ■BlAMDACJ^J , Kf) C 13 1 

Tl AMHA ( Nl . » C 1 S 3 1 -Al. AMOA I U2 i K3 1 1 1 3 r 

/#-F.MOINO OF VECTORlZATION *A 

/«-PFGINtNO OF-S-O-R MFTHOO */ 

T I «TI .AMOA ( NU . DOT .*TX ( NL > 1 

PL AMDA ( cl» K > n 3 « -ALAMDA I K ITT 3 ♦ ( Rlr.XFAC* ( S f Nl 3 -T »7 - 

TX<NL)r-13)f 

!F< KI«H . OR. (vl-^«rH -THEN 

lF<J»n-THEN JJl»?tENDIt 

TF< J»JL> - THEN .I.Ji -JU-l lENDM 

PLAMPA U, K ) r n « -SUM ( Bl AMDA ( JJ. K) C 1 1 U.3 ) /II .1 

-ENOn 

PLAMOA ( 5, 5 ) r 1 1 • 1 1 

^IF ( Al AMDA I J, K4 f T J-0. OI-THFN -Al AMPAI K> t I-3 1 -0.0001 » ENDH — 

ERROR I ■ 1 . 0- < BL AMD A < J » K > C 1 3 / Al . AMDA ( d « K ) Cl 3 > I 

IF- ABS< ERROR ) >TOl RNCE— THFN CHECK* -CHF.CK+H ENDI t 


» t— no 

IJL AMriAr* 11 s AMDA(*J« 5> C361I 


NDFt 

- 


WR I TF- IIN I T ( ^ L I WL AMDA r 1 tOI. 3 ) #5F 1 4 . 741 

POR-K:»V~Tn KL-nO 

Kl .AMDA r K 3 I -Bl. AMDA 1 5L K > r 36 Tt : 

FNDF 

WRME-IINITI6) (KI.AMOAr UKL3 )«‘'KLAMnA»''»5E14. 7#« ~ 

UR I-TE.-I IN T T ( 6 ) < Bl. AMD A ( 55 , .*S ) C 30 * 40 3 H» ' I. AMDA«^ 55E 1 4 . ; 

WR I TF -UN I T (7v> ( PI AMDA ( f=i . 5 1 f 13 . BLAMDA (5.5)11 1, 3 ) 

4^I.AMnA»-'.2El4,7<M 

MMMMr *0 1 

NL.t--0 * — ... — — . - ... 

FOR -..I : - 1 TO JL no 

FOR -Kl «l-Vo KL DO 

Ml l-Nl +11 

Ul AMDA r NLrl I «M AX ( PL AMOA ( J . K4 M 

WLAMDA t C NUl I «M T N ( BLAMDA <J» K )44 

IF— ( MMMM<3 ) THFN 

^IF-lULAMDAl CNLKO.Ol- THEN 

WRITE-J.INI T (6 ) ( J.K. WLAMDA UNL3 • 

iWLAMDAl-'',E14..7«M 

MMMMt «MMMM+t L 

ENDI t FNDI « 

IF( (.l-ST) .AND, (K-U.)-.THFN 

WRITE- UNIT(6>(.J»K,.WLAMDAltNL1)4-' J»''-»-I3,-''K-‘S 

: iWLAMriA)«^F14.7IM 

ENDI I 

ALAMDA (J.K) Ml II -34 «Pl AMDAIcUk >-C) i IL3 1 


FNDF I 


A 
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MAXI AMDAt «MAX ( ULAHnA ) I 

HTNLAMOAl«MIN<W».AHOAt>« 

— WRT TF H IN r T ( 6 H HAXl .AHAA, M I Nl. AMPA ) #' l;AHOA«''*f ' MAX-^-» El 4 ; 7i 

^tN--'-iE14-.71M 

ITERATIONl-tTERATIONi-H 

WR I TF - UNI T 1 6 H CHECK > 4 ^ NUMBER- OF- FI.FMENTS~NOT-CHECKED«' • 



COUNTlNOi -COUNTINO+l I 

WRITE UNIT(6HC0UNTINI3)#'’C0UNTTN0 0F-TTERATI0NB-*,-I8<M — 

^IF(COIJNTTNO>tO)THEN OO-TO-BBBlENOn 

IF-CHFCK>1 THEN 00 TO I.CiF.NDH 

WRITE UNIT(6><ITERAT10N)#1H0.-'NIIMBER 0F-ITERAT10NR-' 



/♦- ENOINO- OF- S-O-R METHOD-*/ 

PRINT-INO -OF-UAMDA - TERMS- */ 

— / »- — FOR -l-TO -^Ib- 00-—-*/ 

FOR -JI--1- PY- 6- TO JI..- PO — 

FOR Kl-l -^TO KL-DO — 

— MRITE-UNIT(6M»1.KH»1H0, , 1341 

WRITE-UNIT( A) (BLAMOA( J» K> Cl * 101 )4PE14. 74l 

ENOFI 

FNOF« ^ 

/* ADJUSTMENT OF -OBSFRVED-WIND FIELD TO MAKE- IT-MASS — 

r^CONSERVATIVE— */ . 

hH.t -Ol 

FOR. Jt»l-TO-Jl.-00 

FOR K I ■ 1 -TO KL-nO 

FOR-n-l- TO -lU DO 

NlJ>NL4-l.t 

ADJC I » 51 1 ■AO.JX < NL-> C 1 1 5-1-1 

M»«r + u 

I2»«I-4» 

Jlt».H»l.« 

.J7I «.J-l « 

K1 X -K+l » 

■■ - - - — 

lltl ( 1 1 » u 4^ u -THFN r 1 1 ■ n FNO I X 

IF-U2-0 UTHFN I2i*IL»ENDH 

IF-J Jl-vIL+l )-THFN J1-»«^IL-1» J2»»Jl.» ENDII 

IF- < J2«0 _ -> ..THFN J2* -1 1 END 1 1 

— IF— IKl -KL+1 J -THFN K1 « -KL t ENOI X 

IF_(K3»0 >-THEN K3*»UENDT» 


AD JUSTU C J . K ) r U l -UU .KU-I J + AD J tJ 3 * C BL AMDA CJ . K > f 1 1 3 

-BL AMDA < J , K > C.T 2 3 ) -AD J C 2 3 * < BL AMDA ( J «.K 1 ) C M 

^-rBl AMDA(J»K3iri3)l 

AD.JUSTV1.J. K) C 1 3 1 ( J, KU 1 3+ADJC31* ( BLAMDAL J1 , K) Cl 3 

^Bl.AMDA ( J2, K ) r n ) -ADJf 4 3 * I PL AMDA ( J. K1 ) 1 1 3 

^BLAMDA ( J, K3 ) C 1 1 ) » 

AOJUSTW ( J, K ) n 3 » -W IJ, K VC 13 ♦ADJC 53* ( BLAMDA ( J. Kl ) Cl 3 

-Rl AMDAl.UKSirTIlt 






— /♦-PRTNTINO OF-AOJUSTFD WIND FFR-D */ 

—FOR ,it»\ TO vlt OO- — 

FOR Ji«l RV Ji no 

— FOR Kl-1 TO KU“00 

— WRITF-UNTT(<v» r3«>^K*S f3«M 

WRtTF UNtT<6)<U<J.K)CtU01,ADJIJSTU<J.K)CU 101» 

V( J»K>ril 101.AnJUSTV(U»K>Clt 103, 

WC.J.K)CJi t01,AnvlUSTW<drK>Cli 103)R5ei4^ 7#t 



«-*CAI;CI.H.ATINO DIVFROFNCF-OF AnJUSTFD WINO FIFI.D -♦/ 

FOR di«l TO-JL-DO 

FOR~Ki«l~TO K\. no r 

U<.lrKKt • TL1»-A0.JUSTU(d,K>C-ti Tblf < 

V(vl,K>n i rL3««AD.MJRTV(J,K>C-n Il..3t 

W< .J, K ) C 1 1 TL 3 1 -AO.JliRTW( .J, K ) C 1 1 IL3 ^ 

FNDF I 

FNDFI 

— CAl.L-niVRGNCF» 


BBRt ENDP t - 

-ENDMI- 

— FO t~ FNCOMNTFRFO^- 
/I DLF 



m 



5J- 


FiLE HASSC 


\ 


i 


1 

i 


K«< nddma fit — "■ 

— — 

♦RFWTMn-MAI^Rn,, 

■ ■■ ....— ■ — ■■■■ 

JSTRJC FU TlOO.CMt ?0000» 

-AT-T/M INT . 

_XHAR0Fv.ii):>?F:3,i.Rr.. : 

«DEU.IVER. R2 1 1 AWWASRFiM . AKHTAR 

-OR T ( MARJ5 1 

—CRT-C OI BTROP) 

-ATTAr.H^Rl. J -0 IN«L TRRARy — 

RFTMirR.- 





I 


RFI.,|?0000. 

Stt-WMASS» W 

ReWACF* miMMY: 

RFWTNn.Si.1BTN. 

TOSTARn*. SMPfN«PI //U.0LRTH0P-/7U1 

DAYFTI F.OAYi>* 

RFPlACFinAY?! 

FXTT# 

RFFI -ACF. I DUMMY 1 

nAVFll.F.nAY^>» 

RFPI ACF I DAY? 1 

— I STORF - <»R9«.50 -1 A??5R MASS P U f 

MASS.T70. ! 

DFI rVFR(R?l I4> 

l-OAD tSI 1 P IN. I~(«SI-I I 1 B »CN«On. OU»l CDF f I #200 \ 

OROl.««C1 .OROI •»C2.nROI.>«C.S) 

ATTACH. CDFF . VFL OCI TY 

RFOIIFSTrFIV-l^.-i ^ 

00. 

DAYF 1 1 F.^nAY4-.. 

XOAST J&«AS9A50F^RFSI IUT4_: 

FXTT, ^ 

DAYF ILF.. DAY 1. 

IDAS (Z>AS9ASCMURFSIJ1.T 1 > 

DSMPnMR.. 

FOUJ^NCOI.INTFRFD — ■ - • 
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file UINDC 


/ r*>«rt l h q » ><« ri>dc- 

- HiRHUTNn..wiNnr,. 

/.copy?ihf.vniindc 

4STRvinB»T400* CMt ?OOCr.. 

Ai:;C:OUNT » 

CHAROR-^ 1 OS^SSU.RC, 

OEUI VER. R?4.1 4 VVUASFEM - AKHTAR 

GET (W T NnP ) 

GEJ4GLRTROR1 : 

___AT-TACH. Sl_17UN=I.IBRARY 

reduce:., r-. 

RFLi 120000 

-SLl.JJINDP, DUMMy..i.SU BIN 

REPLACE. DUMMY 

REWIND. SLl BIN, 

TOSTARU ..SL1RIN:?BT7/U. GLBTR0Ps/7U4 

DAYFTLE.DAY2 

REPLACE (.n AY2 1 

EX 

REPLACE I DUMMY ) ! 

DA YE T LE uD AY2 

REPLACE.I DAY2J 

1 STORE 6S9.^50.._102?5£! .WINDP 

WIND 0.70. 1 

DEI IVF.R(R21 141 
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file WIND 


— P<ki«t) r««( « wi nH — 
«REUTNn.UTND. 
-/t*^w i »- 


♦REMIND, MA;5R. 

/copvshf-, ma jm 

-\ MODI H E- CONT I NMI T V ( UN I T 1 


DLBTROP, UN I T?=COFF . UN I T<>=RFSUl T 

U.INIT.'3t=VEl.OCITY>- BEGIN MAIN? — 

XftCr.^-PROGR AM CONT-I NU T TV - 1 EADS -TO - AN AD^tt lSTED MASS — , «.-C- 

CGNSFRVATIWF-WINO-F-IFI.D-FROM-REAL-TIME-OBSFRVED-OATA. 

^THE OBSFRVFD MIND -COMPONENTS U AND V ARE SUPSITUTEP 

IN THF- CONTINUITY EOUAT ION TO ESTIMATE- W COMPONENT. : 

THE CONTINUITY FOUATION-IS-WRITTEN -IN SPHERICAL 

4-:00RDINATFS with -a -vertical- coordinate system WHICH-TAkFS 

^I NTO ACCOUNT -THE -VAR I ABLE -TROPOGRAPHY-OF— THF EARTH ^ SURFACE 

AND VARIABI.F--TROPOPAUSF HEIGHT.— THIS -PROGRAMME IS ^ 

UR I TTFN- 1 N STAR -I .ANGUAGE -1—1 SL l-»- FOR -CYBER MACH I NE . CC*V- 




-/ *CC---rDEF I N I T I ONS— 

X=FASTWAARn D I RFCT I ON 

-Y=N0RTHWARD D I RECT I ON - 

Z= VFRTI CAU— Ti I RECT T ON- - 


— XX=GRTD-S.T 7F -IN -X -direct I on 

— VY=GR I D S T 7 E-I N -V -D I RECT I ON 

_Z7=GRID SITE- IN -Z-OIRECTION 

—I I =NI.IMRER - OF- GRI D - PO I NTS-I N X -B I RFCT I ON 

, sNl IMBER-OF-GR ID PO I NTS -I N-Y-D I RECTI ON 

_.KL=NUMRFR-0F— GRID POINTS- IN -7 -PI RECT I ON — 

JTHETA -?<_TH=ANGLE -T N-DEGREFS ^ - RAD I ANS RESPECT I VFL V- FROM- 

SOI ITH POI E^ -90 ) -TO NORTH POl F ( +90 ) 

-Ra_D I STANCE -FROM EARTH -CENTER-TO THE POINT OF -INTFREST— 


-IN -TROPOSPHFRF- 


_TROP= -TROPOPAI ISE HE T GHT 

JLl. V_&._W=ll I NDl COMPONENTS - UM-X-. Y- S<-Z-ni RECT-I ONS -RESPECT I VFLY- 
END OF DEFINITIONS CC*/ 


/» THF FOU OMTNr, SFCTTON NFFDS TO PF OHANOFD FOR A GfVFN 1 

ORTO-RYSTFM^FOR RXAMPI.FI-FOR -A ORIO SYSTFii 5- DEGRFF -I.ONOI 1 

TMOF-X— R-RFORFF I.ATTTMOF-K-^ VFRTICAIr-tFVFl S WE-HAVF | j 

«37,-KI -9 ANO TOTAL GRID POINTS I.JK-23976. - XTRAl - TO XTRA3 AND i 

ISMI-T0-ISM3 ARF PAGING PARAMFTFRS WHOSE VAI MFS SHOULD BF ESTIMATED ^ 

FOR A -01 YEN SYSTFM*--IN Tl. IS PROGRAMME VARIARLES-R - »U. V.W AND 

and AI.PHA ARF paged SFPARFTEI Y ON D I FFERENT-L AROF - PAGES . 

FXAMPI.E* - ALPHA VARTARl.E HAS-TOTAL- 13A080-RIT W0RDS-(6X31RX72> . 

THTS -fS -FOUIVAl.FNT -TO- SX < Jl.--?) X (KL— H X IL-. 

A SINGLE- LARGE PAGE CAN HAVE 6RS3A BIT WORDS OR 12S- SMALL PAGES. 

THFRFFORF-WE NEFD THREE LARGE PAGES WHICH WILL BE FII.LFD PARTIALLY 

BY -ALPHA -VARIABLE.- EXTRA PLACE IN-THF SECOND PAGE WILL- BE 

Fil l FD BY A DUMMY PARAMFTFR DUMPS, THE DIMENSION OF THIS PARAMETER 

CAN BE ESTIMATED- AS XTRA3=3X65S3«.-13S0S0«S05-2S AND PARAMETER I SM3 

< IN -OPEN - MAP STATEMENT I «3X 1 2S«3SA. ^ CC*/ 

LITERALLY "S.O" DXXl 

1. ITFRAI L Y -’SS. 0"- DYY? . 

L4 TERALL Y 72 ?* T L » i 

L.ITFRALLY -•137!'_J| t ' 

I^FRALLV^?‘9“ KL» 

M-TFRAUvV-S-' l-O'iKk'LLl /*KKLL=KL-»-l--»/ ^ 

l.ITFRALLY--“-23976" I.JKl 

IrrITFRAI I Y "SIR" I.JK1 X /«- IJKi = < Jl.-2)X-KT -*/ 

L-ITERALLY -••4lR60"XTRAl t 

l^ITERALl Y -yRJI52"XTRA?» ^ 

i L^ITERAI LY ^G0R23"XTRA3L— — 

W-TFRALI V- ••-12S"ISMl X 

LITERAIJ.Y- ••25-k'*ISM25 

I-JTERAI.I.V- "3R4"ISM3« — 

/♦ — END OF- THE SFCTION TO BE CHANGED ♦/ 


— LI TERAI L Y-" Z VECTOR CK 3 7 7 * 

l.I.TFRAL LY- «R (vl, K ) f 1 1 I IL ] " ROl 

— M TFRALI rY -''RTRSC t * I IL-1 *' -TSl X 

— I ITFRALLV- '!Al PHAM ,Nl.)tU-ML3" Al» 

— I .ITFRALLY--!'ALPHA<2»NL.)Ct : Ilt.I*' A2? 

— L ITFRALLV "At. PHA( 3. NL I Cl t I II - A3l 

— i..ITERAI I Y— "ALPHA(4,NL>n t IIL3"-A4? 

— LITERAl LY -''ALPHAIR.NI If 1 : IIL1“-AS» 

— 1,4 TERAI. I Y-i; Al ,PH A ( A . Nl. ) U : 1 1 L 3 A6 1 

LITERAL I. Y- "3. 1 4 1 S92AR3589793" -PI t 

— REAU-VECTOR CKKU.3 7t ^ 

RE Al- VECTOR CKL 3 • 7VFCT0R1 

REAL- VECTOR CJLl TROPt 

REAL. -VECTOR - 1 U. 3- RTRS . RR 1 , RR2 , SS« 


REAL- VECTOR -Cl L 1 -ARRAY ( A > - VEL 1 

REAL- VECTOR r UK 3 WVFl.t 

Z^ft-PAGING OE-VARIABLFS R.U.V.W- AND ALPHA— 4/. 

REAI-. VECTOR - C I L 3 ARRAY ( JL . KU-R » 

COMMON CUR. DUMP 1 > » - 

^ — REAL -VECTOR -X-XTRAI 3 DUMPl 5 

. REAL VECTOR TTI 1 ARPAVf.H .fcvt i ^ ii. «i w« 
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» «• • • « # •NtXS.W • VJi » V » W » 


COMMON C?IU.V.M.nUMP?>t 

Pf- Al - VECTOR C XTRA^ 1 DUMP^t 

PFAI. VFCTOR-MI.J-ARRAY(A. I JKI ) Al.PHAt 

common C3<AI.PHA, DUMP3) I 

RF AL VECTOR t X TR A3 1 HI IMPS i 

-/♦-FND OF -PAOINO */ 

INTFGFR r.vKK.NL ,EI .F?,F3.KHKK. ni;l 

PFAI. MAXW. MINW. 17.1 . THETA, TH, XX . YV, CONt 

TNFTT Al I 7-VO. -0333334 , 0-. 0333334 ,0.0333334 ,0.0333334 

vO. 1 , 0. 2, 0. 3, Oi 2v0. 2, 0r2\ > I 

-TNITrAI--<TR0P“\9.?,9.3.9.4,9.5,9.7,9.-9. 10. 1, 10.3, JO. 5. 10.7, 

10.9, 1t< 1 , It .3, 16.8, t7.-ft, 13.0, 18.0,18.0, 18.0, 

18.0,18.0,18.0, 17r8rl6.8, 1 Jv8, 11 . 1 , 10*9, 10.7, 10.5, 

10.3, 10.:t^9.9,9.7,9.5,9.4,9.3,9.2\)l 

PROCFOURF MAIN* 

OPFNMAP ( <F 1 ' . R ( 1 , 1 ) r 13 , r 8M1 , F 1 , ''l.ARCE-''-) 1 

^ lOPFNMAP ( ■' VFLOC I TY < , U < 1 , 1 ) T t 3 , 1 SM2 , E2 , ' I . AROF. < ) ; 

OPFNMAP( ‘'rOFF'’-.AI,PHA(1 ,1 )f l 3, I8M3,E3,'^l.ARr,F'’-)-l < 

WRI-TE-IJNIT ( 6 > (El-, E2, E3 > H-iOPEN-VS (3X 12 ) 41 i 

M( t«ri 1 

CON: »P t /-1 80, 01 : 

-/♦ 1 N I T I Al . I ? r NG-U - V AND W VEl.OC t T Y - VAR I ABL FS * / 

THFTA: =-90. 0-DYY-l 

TH:=THETA*r.ONt ^ 

FOR . Js at_TO .11. -GO 

FOR-- k ! al— TO - K kU.-nn 

UC.J, k ) r-t-t4l.3 •! » -1 0.0? 

^VlvI.Kin : n 3 J=3.O*C08(TH>l 

U.(,l, K ) r 1 : TUj! = 0.-01 

FNGFlFNnFl- — ^ 

THE FOI I nuING-,STATFMFNT -i'CAU— COFFFICrFNT''— l-FAGS THE 

EXECUTION TO PROCEGURF-COEFFICTENT-^-THIS- PROCEDURE 

F8T r MATES -THE- NECESSARY COFFF I C f FNTS- FOR -THE- CONT I NU T TY 

EQUATION -FOR A GIVEN COND I T I ONS • SUCH 

AS -01. OBAI..— TOPOGRAPHY- OF- EARTH "S-SURFACF-.-TROPOPAUSF- HEIGHT - 

ANG -GRID -SI 7ES. -THEREFORE FOR- A- GIVEN CONGTTIONS-THIS 

PROCFOURE-SHOUl. n - PE ACCESSED -ONI. Y-ONCE.—THF- COEFFICIENTS 

THUS ESTTMATED-WTU- -BE STORED PERMANENTI.Y IN A FILE 

CAM ED-i'COFFii.-IN ORGER--TO AVOIG-EXECI.ITION OE THIS 

STATEMENT -TYPE -/* AND ♦/ ON -I EFT-AND R I OHT STGE OE 

l'jr.AU_ COFFEICIENT'l- STATEMENT, 


CAI I COEFF I C T ENT l 


-THE FOI I. ON I NO GO- 1 OOP- i!-F I RST -CONVERTS VELOCITY 

. GATA IN THE VECTOR -FORM AND THEN USING THE CORRESPONDING • 
-COFFF I C I ENTS FROM FILE- " COEF " - FST I MATES VERT I C-Al —VFLOC I-T Y 

-COMPONENT “Wi -WHICH SATISFIES CONTINUITY -EQUATION.- */- 

— Nl =0 1 

FOR -J? =?-T0- II -1 - DO 

FOR k'l=1 TO kl GO 




K-t * «K- t « 

^IF-(K1«0) -THFN KlfttFNDlr 

VFI;( l)Cl* TL-3t«W(.J,K)Cti IL1« 

vEi 7( ? ) r n I «i M j > K > t ? 3-iM j . K > c n. 3 f 

VFL- (-2 K T I ; 1 1 «l J ( J , K ) C 1 3 -LK . JtK ) C I L- 1 3 1 i 

^VFl.(?»r?« II--13*«U( J.K)f3s IL 3-UUI.K)Cl t IL-21t 1 

^VFl.(3)C-tMI.H=V<J+l»K)Cll ILl-V(.J-l,K)Clt IL3I i 

VF.l (4>M» Ilr3i»VK».K>Cl « 11-3 » 

VFI_ < 5 > C 1 « 1 1 -3 » «U ( J ♦ K+ 1 > 1 1 1 1 L IHJ ( .J , K I ) T I f 1 1. 3 ? 

^vFL(^v)Mi ii 3»=v(.j,K+i >rn ri.3-v(j,Ki>ri j II 3t 

S3C II IL-3s»0.0l 

FOR KKt«J-TO -^.-DO 

SRf 1 HI 3 1 =SSr t HL3-M VFI.(KK) C 1 ! 1 1-3 *ALPHA ( KK . NL Kit 11-3 ) I 



1 F ( K» 1 ) -THFN W ( J , K-M > r 1 1 I L 3 1 «SS mi 1-3 » FNH It \ ^ 

IF(K>n THEN UI(.J.K+l )Clt TL3i*W(..l»Kt >Clt TL3+SSC1: TLltENOn 

FNOFt 

FWOF; - 

^*_TfM-ST4MATF-MAX I MUM -ANO- MINIMUM-VALUE OF-''W"- VELOCITY- -♦/ 

Nl.s^O? : 

F0R-U:=1 -TO Jl- DO 

^OR K t «-! -TO ia-430 

FOR-I ! *I -TO -IL- DO ^ ^ 

NL:*NI„+1-t 

WVFl r Nlv-3 : =W ( J ) r I- 3 1 

ENOF t FNDF t F.NOF » 

MAXU:=MAX(WVFl. ) t 

MI NW : =M T N ( WVFL ) t 

4JR I TE-UN T T 1 ( MA X W . M INW ) # 'MAX W « ' rF, 1 4 .-7 » M I NW= , F 1 4 . 7# t 

Cl . OSF -( ' F l-i., E 1 ) t : 

Cl ORE*. l-iVFI OC I TY E7 ) t 

Cl. ORE (-^CnFF'-tEaM 

— WR I-TF-UN I-T-l A44E I . E3 » 4 'CLOSF»< < 3X 12 ) 4 1 

ENDPt 1 

PROCEDURE-COEFF I C-I ENT 

J.* SEE -DETAIl.ED COMMENT-ABOUT-THUS-PROCEDURR. IN THE 

ILPROCEDURE-MAINiL^ */ 

7. VECTOR C 1 : KL 3 X « < 7 C 2 X KKLL 3 + Z C I : KL 3 ) /2 . 0 1 

XXteDXXfrCONt 

yY:=OVY«CONX 

^CC^T-_I NI.T.I ALI 7 1 NO ELFMFNTS-OF- R ( JL, KU I II J --~CC*/ 

FOR-JX-Sl-IO -JL..DO 

READ-UNIT(4J (R( J. l» C U II.3J4AEl^. A4t 

R( J...l)Clx IL3X«6371000.0+R(J. UClt IL3 x 

RLJj KL .) C.LUUJ =637.1000. 0+-TR0P-CJ3 *1000. 01 

TS1X=R(.J,KL>C1X IL3-R(J, UClUL.3.t 

7.ZZ.x.=0.0t 

FOR_K;-=2_TO_KL-DO 

. . 7ZZte777+7rk'ii 
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ROi»R( J, I >C1 1 TI.]-«-7Z?#TSU 

FNHPt- FWnFt — ^ 

/♦ FRT I MAT I NG - THE COEFF I C I ENTS -OF’ CON I TUI NY EQUAT I ON T/ 

Nl.i«Ot — 



FOR vH -2- -TO Jl.-t- DO 

THETAl 



FOR- K» «1 TO K'l -no — 

«NI 

TSl ^»R^J.k•l.)t l t II ]-R(J. 1 irii H.-3t 

RRt Cl-l:=«R<rJ,K-> CZ1-R( J»KK-IL3t 

RRi r iiri»»Ri..ivK)cn-R( J*K)f iL-nt 

RRir-Zi Tl.-n» »R(..UK>C3i TL1-R<.J»K>C1* TI.-2II — 

RR2f-U I Lr3 8 =R ( -.1+ 1 -K > 1 1 M I. ) -R I K > Cl « II. 3 I 



A1 J =IR0-?.0«TS1 *Z7> /ROl — 

A2:=-T:51*ZZ/(2.0*RO*COS(TH)*XX)I * — 

A3« =-TS1 *77/(2.0*R0«-VYM 

A4 1 =-•? . 0*Y Y» A3*T AN ( TH >-l 

FNDi t — : 

TF(K>1 )- THEN 

A1 4»-4,0*TSHi77/ROt 

A2 * =-T3 1 *7 Z7 ( Rn*COS < TH ) *X X ) » 

> A3: =-fSl */7 /(RO«VY) V 

A4 , 0*V Y*A3*1 AN (-TH7-? — 
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Each NWS data field is prefixed with a 384-bit label. The first 192 
bits contain values which indicate the quantity, the surface or layer, and 
the time of the data being labeled. The next 64 bits contain the initial 
time, the number of data points, and generating code information. !lhe re- 
maining bits hold, among other things, scale factors partinent to packed data. 
Appendix C shows each entry and its location within the label. 


All data fields can be thought of as some quantity, ^ on a horizontal 
or quasi-horizontal surface S at some level of value L. If the quantity 
is by its nature, a layer-defined quantity, the limits of the layer are 
given Si at In and §2 at 1^. ' Code figure for Q and S are both given in 
Table 1. 

The numerical value of the level 1^ for the corresponding surface S_ is 
composed of two numbers : an integer £ and a power of ten with the exponent ' 

E such that 

L = C x 10®, 

Both £ and E are signed integers (sign and magnitude). The space allocated 
for £ in the identifier is 20 bits which is large enou£^ to accomodate the 
binary equivalent of any 3-digit decimal integer. !!he convention which must 
be followed in order to provide a unique bit configuration for any given 
Icvci is as follows: 

“N 

(a) ■ C must be a 5-dirit decimal integer whose leading digit 
(the highest order digit) must be nonzero, unless L is 
zero. 



Cnee £ is determined, the value of the exponent E naturally follows such 
that toe resulting L will have its true value in the standard units given in 
Table 1. For axanq)le, if £ is 500 mb 

c = 50000 

E = -2 

50000 x 10-2 = 500.00 

# 

or if L is 0.7 mb 


c = 70000 

E - *5 

70000 X 10-5 = .70000 


Appendix B contains the hexadecimal equivalents of commcnly used values of £. 

Table la gives the values for the marker R which is used 
t for identifying different spectral quantities. 

I Table 2 giv».s the values for the marker £ which is used 

t ■ ' for identi''ying fields involvini: time and explains how 

■ • and Fjj arc to be used for each value of £. 



Tabic 3 gives the values for the marker m which is used for 
identify ng fields involving layers* ** 

Table U gives value to be used for the marker X* In 
general, this marker will be non«4ero only for fields in 
a guess file. 

Table Ua gives the values for the markers ^ which 

identify climatological data. 

Table 5 gives the values for the markerd K which identifies 
the grid to which these data apply. 

Table 5a gives the values for the markers which ixidi- 
cates if the grid was derived from spectral or other 
special methods. . . ' 

fihe initial hour I of the forecast (or the observation time on'which 
the analysis is based) is entered in accordance with the 2U-hour clock GMF 
to the nearest hours. For example: I » CO for midnight GMT and I ■ 12 for 

noon GMT. Y is the year within the century. For exan^le, for the year 1973, 
Y_ - 73. M is the month of the year, and D is the day of the month. 

Table 6 gives the values for the marker R which identifies 
the rm within the cycle. ' , 


Tabic 7 gives the values for the marker G assigned to the 
program which generated the field. v 

Each datian point value on a surface S (Qs) is scaled according 
to the equation . 




A 

QS 


• V 

A)/2U. * 




Qs must be in ^he standard units given in Table 1 . 


The procedui'e followed for scaling is to scan the data to find the 
ma^mum end Mnixaim values of Qs. The value of A is set equal to [(Qs) max. 

* tQs) Qln. ]/2. An integer n, is then chosen such that 2^ is the smedlest value 
whi^satisfies the condition (I(Qs) max-A3/2’*)<l. This procedure allow max- 
imum accuracy to be maintained in packing the data. 


t ■ , a 

The mid-range value A is placed in the identifiers as a IB(4 360 rehl 

number, i.e., 8-bit sign and exponent and 2U-bit fraction. The scaling value 
n is inserted as an IBM 360 half-word integer, i.e., l6-bit integer, 2's 
complement if negative. 

A , 

The scaled quantity Qs, on unnormalized fraction . is packed into 16 bits 

using 2*3 complement if negative. 
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Procedurea for Creating or Unpacking FGCE Level III Data Seta 
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General 

Tha docuaentation for FGCE Level III fotnata is described in Appendix 11* 

• • • 

FGCE International Data Managenent Plan, "Format for the International 
Exchange of Level III Data Seta During the FCGE." 


■ • 


• . • • ,■ 


'i-V. 

• I . 


Purpose '■ 

. i‘ . • . • • 

■ a ■ 

Thia note describee the basic ateps needed for creating or unpacking Level 
III data sets with the use of attached subroutine listings. These sub- 
routines were written for use on IBM 360 machines that have 32-bit words. 
Some modifications may be necessary to utlllxe these routine on other 
machines. 


Level III Tape Files 


V- ->• V 



File I — TEST file 

. . • , 

File 2* — TAPE HEADER file 

File 3 ’ ' — GRID DESCRIPTOR file 


_■ ■■■■ : 



% ...^ 




f f* I 


Fu.. - tEFFL XII MIA fii.(.) 

To Create Level III Mafntetic Tape 

The basic sequence of events in s computer program should be as follows, 

I 

assuming that the required analysis fields are available In an analysis 

file and are ready for processing: 

A. Write the required TEST file on output tape (File 1). 

B. Screen available analyses to determine first and last dates to be 
processed and written on the output tape. Prepare the documenta- 
tion for the TAPE HEADER file with appropriate dates (File 2), 
and write it on the output tape (See example. Attachment 1) . 






p;.- 'I*. « •.*<VV: *1 •• •;.-, ^ ^ ’ll 
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. C. Describe the grldCs) used for the Level III data (See example, 

i * 

Attachment 1). Write GRID DESCRIPTOR file on output tape. 

J De Select an analysis field from the analysis flle^ prepare the 




J' necessary unique Identification words ( see subroutine W3FI32,l 

Attachment 2), scale and pack the data (see subroutine W3AI00, / . Sit/ 

MY 

Attachment 3), and write the packed field In a LEVEL III DATA I (f 
file on the output tape. , 

E. Repeat D above for other fields for the given date/time. 

F. Write an end-of-flle tape mark on the output tape. I 

G. Repeat steps D, E, and F above for any remaining date/times. / . - 

H. Terminate the output tape with at least two (2) consecutive / 


• - ••u. • 




end-of-flle tape marks. 


To Unpack Level III -Data 


:. r:c ' 


The format required to read FGGE Level III data sets/ls written In clear 


text In the second file (TAPE HEALER) of each tape. The third file of 






each tape will contain information concerning the arrangement of the data 
for the grld(s) used In the DATA fUe(s) that start with the fourth file. 

. A. Read a record from selected DATA file. . ' • •',/ ‘ r 

B. Determine if you desire to process and unpack the record. 

' Subroutine W3FI33 (Attachment A) could be used to convert the 

first 256 bits of the record identification to Indi-.-ldual field 
Identifiers. 

C. If the field Is desired, t!.a packed data (beginning at bit 385) 
can be unpacked and rescaled with the use of subroutine W3AI01 
(see Attachment 5). If not desired, repeat A and B above. 


O 



Jh*'* •^.•-*1' 




k sample FORTRAN program Co locate and unpack N. Hemisphere 700 mb heights 


and to print out some values over the United States, Is given In Attachnent 
6* , Descriptions of buUc-ln and Intrinsic functions used In the various 


subroutines are given In Attachment 7 
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APPENDIX 1 



FCGE 3 A 00317601020078010812 - ■ • ' 

10780 FIXED LENGTH RECORDS (BLOCKSIZE - RECORD SIZE) 7V *^ Q 

2CA4 USE ( 100(27A4) ) FOR FULL RECORD ' ‘ u » 

0123456789-:> /STUVWXXZ.(-JKL4N0PQR*);+ABCDEFGH.) (< 

9-TRACK, 800 BPI : • : • ' 

BINARY, ODD PARITY 

IBM 360 / 195 , 32 BITS/wORD, 8 BITS/dYTE. • •* (1 ' ' • 

NATIOIWL METEOROLOGICAL CENTER, NWS, NQAA, WASHINGTON, D.C., USA 
■ ( 100(27A4) ) - . • . • ' 

HOUGH FUNCTION ANALYSIS METHOD, GLOBAL ' (fXATTERY) 

FIRST GUESS COEFFICIENTS DERIVED FROM A 9-IAYER PRIMITIVE EQUATION 
FORECAST MODEL ON A 2.5 X 2.5 LATITUDE /LONGITUDE GRID. 

100. • (42640000) HEXADECIMAL FLOATIHO POINT REPRESENTATION 

-100. ■ (C2640000) HEXADECIMAL FLOATING POINT REPRESENTATION 


(C2640000) HEXADECIMAL 


FLOATING POINT REPRESENT.\TION 


• ANALYSIS FIEIDS NORMALLY PROVIDED: 

.12 lEVEIS OF U- and V-WIND COMPONENTS, TEMPERATURES, AND HEIGHTS 

(AT 1000, 850 , 700 , 500 , 400, 300 , 250 , 200 , 150 , 100 , 70 , AND 50 KB); 

.. 6 LEVELS OF RELATIVE HUI-HDITIES (AT 1000 , 85 O, 700 , 500 , 400, AND 300 ^S) ; 
TROPOPAUSE PRESSURE AND TEMPERATURE (MODEIED) ; SNOW FIEID) 

SEA SURFACE TEMPERATURE (FROM SATELLITE DATA) ; SEA lEVEL PRESSURE. 

HISSING ANALYSES: NONE .. ' ; • - • 


. ... END OF- TEXT 
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EXAMPLE “ GRID DESCRIPTOR FIIE 


GRID 029 RECTANGIE 

GEO(SUU>HIC 

RIGHT-HANDED 

GRIDLINES o 037 , DELTAJ - 2.50, FIRSTPT-(O.OON, O.OOE) 
JdVLL, NI*»145, DELEAI - 2.50 
... END OF TEXT 


GRID 030 RECTANGIE 

GEOGRAPHIC 

RIGHT-HA1H)ED 

CRIDLIIJES - 037 , DELIAJ ■ 2.50, FIRSTPT»( 90»OOS, O.OOE) 
JsALL, NI»145, DELEAI - 2.50 
... END OF TEXT 
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W3FI32 
Pack ID 

(Modified for FGGE Level III) 


FORTRAN H EXT + 
IBM 360 
NMC, DCA 



;v PURPOSE; : ‘ . 

To convert an array of the 27 data field Identifiers Into an array of the 
flrat 8 identification words. (See Appendix 11 » FGGE International Data 
Management Plan, "Format for the International Exchange of Level III Data 
Sets During the FGGE." 


. USAGE: 


Callliie Statement 


CALL W3FI32(KARRAT,PKDMT) 


where. 


KARRAY ■ array containing the 27 data field Identifiers 
(nrrEGER*4) . . 

PKDNT " array to receive the 8 Identification words. 


NOTE 1' 


(a) If any number n In (KARRAY (I) ,1-1,27) is erroneously large, 
W3FI32 will print: 

'VALUE IN KARRAY(I)=n ID TOO LARGE TO PACK' 

(b) If any number n In KARRAY (I) Is erroneously negative, 
W3FI32 will print: 

'VALUE IN KARRAY(I)»n SHOULD NOT BE NEGATIVE' 

(c) If either (a) or (b) occurs, that portion of the packed 
word corresponding to KARRAY (I) will be set to binary 


Examples 

1. Suppose KARRAY(4)-30. Then W3FI32 will print 'VALUE IN KARRAY (4)-30 

IS TOO LARGE TO PACK' and will place ones In the first 4 bits of 
PKDNT (2) . * 

2. To form PKDNT from the 27 data field Identlfers for 500 mb height from 
FMANL (limited-area fine mesh analysis file) dated OOZ March 15, 1974, 
KARRAY should be Initialized as follows: 


i 

% 

\ 



IDENTIFIER 


KARRAT 


KARRAY 


IDENTIFIER 


KARFAT 


IDENTIFIER 


( 2)-8 • - 

(3) -0 < 

(4) -0 

(5) -50000 

( 6 ) — 2 . - 

(7) -0 : . 

( 8 ) -0 

(9)-o r ; 


Unused 


Umtsed 


Then the resulting FKDNT array would be as follows (In Kexadecimal) 


(5) 00000005 

(6) 000005F3 

(7) 4A030F00 

(8) 00130BCD 


'i) 00100800 

(2) 00C35082 

(3) 00000000 

(4) 00000000 






.*? ; 


" 7*.*^V ^l<' 

• . • ■ 'X **y 

4. ' t%S 


U3PI32.' 



SUBROUTINE 

SUBROUTINE 


SiHSi 


nti.run°nLo identifiers. 
ciltlbSIlSiiitSgSS^Il.ZOOICO^OE 

3«7OOU0dO3«7000B0CU3«ZO0OO0BO3*Z 


DIMENSION KIONT (B) * I TABLE (27) 1 1 ARRAY i^Zi 
Data iTABLE/zooiAOCoi tzoooeocoi tZoooooBO 
lZoloOO8O2.ZO0ICOAO3tZO0U0dO3fZO096OCO3, 
220l08lA0AfZ010008OA,Z0Ol80dO|tZ0Oi0pd05t 
3200 1C0A06* ZOO 100C06» 20000 1006* ZOO 180807 1 
4Z00000807* ZOO I 8080 b* ZOO lOOdOdfZOOOO 1008/ 


Z010S1402* 

oIcu<»o^* 


A ITABLE/Z001AOCOI*ZOOOBOC01 *^OOOOUHUi*^OQli.uHuc*4Viy?*7'' 

000802 * zoo iCOAO3*ZO0 1^0803 *Z00060C03*ZOOOpOtt03*ZOOlCU40A, 

081A0A *2010008 OAtZ0Ol8p8o5*ZOOlOO8O5*ZO00 80 805*ZOOOOpflO>* 
lCOA06lzO0lO0C06;2OOg0l0O6*ZOOl808O7 *200100807 #20008080/# 


DATA KX/2FFFFFFFF/^ 
DATA MASK/2000000FF/ 
MAKE KIOENTS*0 
DO 20 1=1*8 
KlDNT(I)sO 

continue 

IS16N=0 

DO 40 1=1*27 

ISC=ITA0LE(I) 


n=LANO(ISC.MASK)^ 
i2=LANO( SHFTR(ISC# 8) 
13=LAND( SHFTR(ISC#16 
I4=LAND (SHFTR < ISC*24) 


8) *MASK) 
16) *MASK) 


5 IGN TEST. 
Iv= 1ARRAY(I) 
IFdV.GE.OlGO TO 12 
1FU4.NE.0) GO TO 1 
WRITE (6*22) .... 

> format (///* 1X#» VALUE 
/TIVEV//I 
GO TO 30 
1 IVTlABsaV) 


I MASK) 


KARRAY(*#I2*Ms»#Ill#? SHOULD NOT BE NEGA '.) 


YSI6N=RSI6N . . : , 

• K=I2/4 

DO 25 M=l*K . . ' . • ^ ^ 

25 ISlGN=SMFTL(ISiGN*AJ - . . .... u, 

. IS1GN=SHFTR(ISIGN*1) 

IV=L0R(IV*ISIGN) . 

12 continue 

in SH?TR(IV*I2) .EOrOIGO TO lOO 

21 FORMATt///*{x* • VALUE IN KARRAY ( • # 12* • ) =• # I 11 # • 
/CK»///) 

30 IV=KX 

IA=32-I2 . i 


LARGE 


I 


IV=^HFTR(IV*IA) 


KIi)NT(I1)=LOR(KIDNT(I1)# 

continue 

RETURN 

END 


SHFTL(IV#I3)) 





W3AI0j)^ 


Vi V V/f^v v: 


•7V.S 1 


APPENDIX 3 



SUBROUTINE W3AI00 <f)EAL4«PACK«LAQEU 


c scale and convert a field of Rf al*a data to malf-worq integers and 

C PACK TOGETHER WITH LABEL DATA IN FCGE LEVEL III DATA FORMAT* - 

( CALL H3AI00(REALA*PACK*LABEL) ‘ 

C WHERE REALA a ARRAY OF REAL»A DATA TO BE PACKED. „ V 

C PACK a OUTPUT ARHAY. WORDS 1-8 COPY LABEL WORQS i-8* WORDS C ' 

C 9-12 are GENERATED* WOR05 13*14,... AHt PACKEO^OATA. 

C allow l2^(J»l)/2 FULLWORDS. WHERE J IS NUMUER^OF POINTS*.., 

C LABEL alNPUT LABEL DATA. WORD 8* BITS 16-J31* MUST CONTAIN THE 

C COUNT OF THE NUMBER OF POINTS IN ARRAY REAL4 FOR THIS 

C . ROUTINE TO WORK. ■ -/ju 

c generated data to be found in Pack ^ • Vk- 

C WORD 9 BITS 0-15: NUMBER OF BYTES IN WHOLE RECORD. • 

C BITS 16-31* EAClUSIVE-OR CHECKSU;! OF- WMOlE RECORD lEXCEPT t':f 

C checksum ITSELF) • BY HALFWORDS* ^ .. n'* 

C WORD 10 BITS 0-31* RtAL*4 CENTERING VALUE A a THE MEAN OF THE !HAX ' 

C ‘ ^ ANO min values OF ARRAY REAL4. v- 

C WORD 11 8JTS 22“§?J tNT?GER*2 SHIFT VALUE N, THE LEAST.INTEGER SUCH 1 


;C WORD 10 BITS 

C WORD 11 BITS 
C BITS 


C ^ BITS 24-31* iNTEGER*2 SHIFT VALUE N, THE LEAST INTEGER SUCH 

C that ABS<X-A)/2«*N is less than i FOR ALL X IN 

C - , array REAL4. N IF.'NEGATIVE IS IN 2»S COMPLEMENT 

C FORM. VALUE OF N WILL NOT EXCEED 127. 

C WORD 12 BITS 0-31* ZERO 
C WORD 13 BITS 0-15* FIRST PACKED DATUM 

C BITS 16-31* second PACKED DATUM , • ?'. 

C WORD 14 BITS 0-15* THIRD PACKED DATUM ' ... • 

c • • • 

C THE binary point OF THE PACKED DATA IS CCNSIOERED^TO BE TO THE RIGHT 
C OF THE LEFTMOST, OR SIGN BIT. DATA IS .Itt 2*S COMPLEMENT FORM IF 
C NEGATIVE. IF PACK IS EQUIVALENCEO TO AN INTECER<»2 AWRAY HLF, THE 


C DATA MAY BE RECONSTRUCT EO.,AS FOLLOWS* 

C REAL4<J) a HLF(J*24)*2.**(N-TS)^;*^A^^ , 

C WHERE A ANO N ARE TO Bt FOUND IN WORDS .10 AND 11. 


. . REAL*4 REAL4(1) : 

INTEGER*2 LABELa)?PACK(I)*IA(2) 
equivalence (A* IAU) ) • (X*1X) 

? TRANSFER LABEL DATA TO WORDS 1-8.- GET 4i£H0C0UNT* COMPUTE BYTES. 

^ OO^IO 1^1*16 

PACK(I) = LABEL (1>* 

10 CONTINUE 

3 = LABEL(16) 

M a J*24 
PACK (17) a m*M 

PACK (18) a 0 • • 

c •• 

C FIND MAX# MlN OF DATA* COMPUTE A AND N« , ' 

RMAX a PEAL4(l) 

RMIN a ,,MAX 
DO 20 Ia2,J 




<<;«•?••* • rrs - j ** 


r4^> 


f 




• 1 ‘* 


■/ 

W3AI00 


. ^ r . 
f y 




- 90 - 



20 


30 


RMAX s AMAXl (RMAX*REAL4(I) ) 
RMIN a AMlNKRMlNtHEALAd) ) 
CONTINUE ^ ^ 

A = 0.5^(RMAX»RMIN) 

X = RMAX-A 

N * LAND(SHFTR(IX*24» *1271 
N a 4<MN-64) . 

IF (TBIT(1X«8)) GO TO 30 
N - N-l 

IF tTBIT(IX*9)) GO TO 30 
N » N-1 
IF_(TOlT 

N = MAX0(-127*MIN0(127*N) » 
PACK (19) a lA(l) 

PACK(20) = IA(2) 

PACK(21) a 0 
PACK(2?) s N 


(lX,10n GO TO 30 


PACK (23) a 0 
PACK (24) a 0 


i - 





C NOW PACK UP the DATA 
C 

' . TWON a 2.««(15-N) 

DO 60 1=1. J 
X * (REAL4(I)-A)*TW0N 

60 TO *0 

. . K = MIN0(32767.K) 

• . GO TO 50 ^ 

. s8 PACkT?? 24? a K 

60 CONTINUE - •. 

COMPUTE CHECKSUM AND STORE 

IXOR a 0 
DO 70 lal.M 

' ?x5R"^^l!iii(IX0R.K) 

70 CONTINUE 

PACK (18) a IXOR 
RETURN 
. END 


f 

C 


■J 

:d 

••■•1 

.1 


A3 - 2 
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. . • W3PI33 
Unpack ID 

(Modified for FCCE Level III) 


ju-ruu^xA •» 


, • • ♦ r ■ • • ■ * #■ 

purpose':’ ’ ■ ’ '■"•• 


FORTRAN H EXT + 
IBM 360 
NMC.DCA . 


To convert an array of the first 8 ‘Identification words Into an array of 
"27 data field Identifiers. (Sec Appendix 11, FGGB International Data 
Management Plan, "Format for the International Exchange of Level III Data 
Sets During the FGGE.") 

. • * 
Calling Statement • . • « .. 

. • CALL W3FI33(PKDNT,KARRAY) 



where, f - • ‘ 

’ PKDNT array containing the 8 Identification words. ' 

KARRAY " array to receive the 27 data field identifiers (INTEGER*4) • 

Example ' * . ■ • • • •. 

Suppose the 8 identification words for 500 mb height from FMANL ‘ \.‘t. 

(limited-area fine mesh analysis file) dated OOZ March 15, 1974 are 
, given (In hexadecimal) ..... »•..•• 


PKDNT (1) 00100800 
(2) OOC35032 
. (3) 00000000 
(4) OOOOOOOO 


(5) 00000005 

(6) 000005F3 

(7) 4AO30F00 

(8) .00130BCD 


Then the resulting KARRAY array would be as follows: 


KARRAY • 1DE^^TIFIER KARRAY IDENTIFIER KARRAY IDENTIFIER 


(1)-1 . 

Q 

(30)-0 

(2)-8 . 


(ll)-O 

(3)-0 


(12)“0 

(4)-0 

t^ 

(13)-0 

(5)-50000 


(14)-0 

(6)— 2 


(15)-0 

(7)-0 

m 

(16)-0 

<S)-0 

X 

• (17)-5 

<9)-0 

S2 

(18)»0 


^2 

(19)-0 

Unused 


(20)»1523 

NU 

C- 

(2D-74 

JJ 

E? 

(22)*»3 

MM 

c2 

(23)»15 

YY 

CM 

(24)-0 

GG 

Ks 

(25)-0 

R 

R 

(26X9 

G 

Unused 

(27)=3021 

J 



■ . . . ' . 


• T 


►*.; K' 1 ••*■ 
..•> . ?:» 


. 0 • . t 


W3FI33^ 



SUBHOUTINE W3FI33 ( IONT*LARRAY) , . . f 

,* OIrlENSION LARRAY (27) , IDNT (8) t ITA0LE127) I J<6) . JS(6) ‘ , 

^ SUdHOUTiNE W3FI33 UNPACKS THE 8 IOENUFICATION WORDS INTO AN ARRAY;. 1 
■ OF 27 DATA TIELO IDENTIFIERS. „ ^ ^ _ ‘v’ 

; DATA J/ZOOOOOOOF*2000000FF#ZOOOOOFFF»20000FFFr fZOOOFFFFF* 

^SaTA^JS/Z00000007*20000007F,Z000007FF*200007FFF*Z0007FFEF* . : . K 

‘ . ^SX?A^MASK/Z000000FF/ • • I’i*;; ' 

’ - DATA ITA8LE/ZOQ1AOCOUZP0080C01*ZOPOQ0801*ZOOICOA02*Z01081A02* ’ 'v;" 
l2Ol0O08O2t2OOiCOAO3,ZOOlAO8O3.ZOOp0CO3*ZOOOQ08O3*ZOOlCyAOA, • . 

~«Z001006Ub«Z00080805«Z0000080b« 


IZO 1000802 1 ZOO iC0A03, ZOO 1A0803. ZOOQ80CO3*ZOOOQ08O3» ZOO ICUAOA, 
2Z0l0el40Af Z0100080A, ZOO 180805, ZOOi006Ub«Z00080805«Z0000080b, 
3Z0U 1C0A06* ZOOl 00C06,Z0000l 006, ZOOl 80807, Z0010080Z«Z00060807, 
AZ00000807, ZOO 180808, 2001 00808,20000 1008/ 

00 50 1=1*27 
ISC=ITABlE(I) 


* I1»LANQ(ISC,MASK) , . . 

I2=LAND(SHFTR(ISC*B) ,MASK) 

.I3=1.AN0( SHFTR(ISC*16) *MASK) 

-H I4=lAN0(SMFTR(ISC*2A) ,MASK) • ~ * 

• * IX=I2/4 

LAHRAY(I)=LAND(SHFTR(I0NT(II) *13) ,j(lxn 
'■ TO SEE IS THE NUMBER IS MINUS 
IF(IA.EQ.O) GO TO 50 

' IB=LARRAY(I) . • 

IC=SHFTL(IB,1) 


K=SHFTR(!ij4) 


‘ IF(IC.EU.O) 60 TO 50 

abs. value, then minus. 

^ LAHRAY(I)=-<LAND(lb,JS 

50 continue - - . . ‘ 

RETURN- 

-,^ENO • .. • ; 


(IX) )). 


OKIGINAL PAGE IS 
(71' POOR QUALITY 
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AU . 2 
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APPENDIX $ 


SUBROUTINE W3AI01 (PACK*RE^4«LABEU 


UNPACK AND FLOAT A FIELD OF PACKED DATA IN^FGGE LEVEL^III DATA 
FORMAT ACCORDING TO bPECIFICAT IONS IN THE ID WORDS OF THE 
PACKED DATA. ALSO MOVE THE ID WORDS TO A LABEL ARRAY* 

.CALUW3AI0UPACK*REALA, LABEL! 


where pack 

HEAL4 


Array of id words and packed data to be unpackeoa • 

REAL*4 array TO RECEIVE THE UNPACKED DATA WHICH WILL 


BE DE-SCALED ACCORDING TO ID SPECIFICATIONS* _ 

LABEL s A 12-sWOHO ARRAY INTO WHICH THE 10 WORDS WILL BE COPIED* - 


C PACK MUST contain THE FOLLOWING FIELDS. .*.• 

C WORD B BITS 16-31: NUMBER OF DATA POINTS J. 

■C WORD 9 BITS 16-31: EXCLUSIVE-OR CHECKSUM OF REST OF ARRAY BY 

HALFWORDS* THUS CHECKSUM. OF ENTIRE ARRAY SHOULD wt 
BE ZERO* . ■ 

W0r 8 }l bit! iS-il: INlkoERSz^SHlFT^VALUE'^Nf IN 2»S COMPLEMENT FORM 


WORD 

WORD 


WORD 13 BITS 


C WORD lA BITS 


BITS 16-31: 
BITS 0-lS: 


IF^NEGATIVE* 
0-15: DATUM 1 


DATUM 2 
DATUM 3 



2 THERE- WILL BE 12^(J*l)/2 FULLWOROS IN ARRAY .PACK t ANQ J FULLWOROS IN 
:C ARRAY REALA. 


»FAL*a rfa 
NTE6ER*2^ 
‘QUIVALENC 




(1) »IA(2) 


C TRANSFER ID WORDS TO LABELo . 

C 

DO 10 I-1*2A 
LABEL (I) = PACK(I) 

10 CONriNuE 

8 get word COUNT# A* N* . 

J - PACK(l6i 
IA(U a PACK(19) 

IA(2) a PACK(20) 

NS PACK (22) 

TWON.s .2***(N-15) 

8 unpack* convert REAL*A DATA 

c 

do 20 IMal, J 
I = J^l-IM 

REALA(I) a paCK(I*2A)*TWON^A 

20 continue 

RETURN 

END 


O 
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0 
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z 
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mmmm 

mm 

m 
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00 00000000 00000 oooc<o*r 
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ORIGINAL PAGE IS 
' <W POOR QUAUTY 


APPENDIX 6 










i6- 


• ••• 

APPEHDIX 7 

Built-In t\inctlons * ’• . * ‘ *•••■ 

IV - LAND (a, b) • v . 

, * . ’ • • , 

■ > vhare.a,b may be a 1-, 2-, or A-byte logical Integer expression. 

The value of LAND Is obtained by AND- In g the Individual bits of the 

arguments. The resulting value, IV, will be considered to be loglcal*A, 

but may be used as an Integer. 

,IV - LOR (a, b) , ... 

where a, b aay be a 1-, 2-, or A-byte logical or Integer expression. 

The value of LOR Is obtained by OR-lng the Individual bits of the arguments 
The resulting vclue, IV, will be considered to be logical*A, but may be 

used as an Integer. ' ' ' 

•* • -• ‘ . • • • . .., •. 

IV- LXOR (a, b) 

\diere a, b may be a 1-, 2-, or A-byte logical or Integer expression. '* 

• 

The value of LXOR is obtained by exclusive OR-lng the Individual bits of 
the arguments. The resulting value, IV, will be considered to be logical*A 
but may be used as an Integer. .. 

IV - SHFTL (J, K) - • V 

T.V - SHFTR (J, K) 

>• 

where J la a A-byte variable, 

K is the actual number of bits to be shifted. 

The values of SHFTL and SHFTR are obtained by shifting the first 
argxioent, J, left or right, respectively, the nuc&er of bits specified 
by the second argument, R. The resulting value, W, will be considered 
to be *i.ogical*A, but may be used as an Integer. 
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.. 


./ i-i ' 


;r . 




W:. 


IV - TBIT (A, K) 


/ 


where A la a variable* 4-bytes or less* 

K Is the number assigned to the bit to be tested 
The value of TBIT la .TR'JE. or .FALSE, depending whether bit position 
K of the variable A Is OW or OFF (ON-1, OFF^O) . Bit 0 Is the leftmost bit 
of variable A. The result, IV, will be declared as loglcal*4. 
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J . 

) * 


Intrinsic Functions 
AMAXl - Maximum, Real 

A»AMAX1(Xi,X2 X^) 

where A and X are real *4 

AMINl - Minimum, Real 

A-AMINKXj^.X^, . . . ,X^) 
where A and X are Real *4 

MAXO • Maximum, Integer 

I^MAXOCM^, M2,...,Mj^) 

where L and M are Integer *4 

• 

MINO - Minimum, Integer 

D-MIN0(Mj^,M2,...,M^) 
where L and M and Integer *4 


lABS - Absolute Value, Integer 
L - lABS(M) 

where L and M are Integer *4 


SIGN - 


Y-SIGN(Xj^,X2) 

where X and Y are Real *4 

Y-(slgn of X 2 ).|X 3 l| s 


c 


O 


A7 - 2 
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V ' 

* ^ ' 

September 18, 1978 


RV^Mi) * Woiurnal Occnnic unJ 

1 |:NV1^U1NM^NTAL DATA SinVlCC 


(lattci'.fll CHwHc CenUT 
(chloral ruOd»n4 
A^hevUl**. H* t. 26B01 


' • ' w* h’ C' 

D523107136)/SD 


TO; M*** 

NASA , ^ ^ 

Lon^ley Rcsoerch Center 
Hampton, VA 2366b 


‘I . ' 


i.-n -oLc^, 


> n 


FROM: Cd^Richav<l^*' DaviS,\\iief > 

^ Computer Products oranch 

SUBJECT: Tape Dumps (Your Tel con 9/8/78) 

Enclosed are hex^pe dumps of the first and last five records from 
the NNC global and N. Hemisphere tapes for AprilJ-15, 1978. 


Enclosure 

V I V V.' \ 


> . I 

a 


a 


•• i. 


I ' I 


,'L , i 






rv. ; 7 - JC 7 ^ 


\ , / — * 

'j ))>.• •■* /’. **-*^V>’ 

• 





, » t.(0 • o ' 




A-i. U -■• 


1 1 * ) ,'t • • .*0 


• 

/ 

>_a •' 



' . ' V-, 1 1 


Jiv" 

( “ 

J • t 

t- 

/Nil'!' 
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DECODE Proqran Listing 


/♦ (SUDRO.,STR(;,STRZ,SrZE,FOFL,OFL) ;•/ 

nFcnnn: /♦ to iimrack kasa weather tape •/ 

PROCIIDTRE 
CPTION.T ('»MN) : 


/* 

Thi:^ pro.jran rpiuiros no sp*'clal control c^rds to ilecod** 
till FO;IE level III packed on a National Noteorological 

Center weather tape. Proper .ICL cards nust be included, 
however, SYSIJTI and SYSHT3 are the DDnames for the input 
and output datasets respectively. If only a part of the 
data is desired, some program modification will be necessary. 
The basic conversion subroutines nay renain unchanged. 

See the writeup, UNPACKING FGGE LEVEL III DATASETS IN 
INTERNATIONAL EXCHANGE FORMATS, for further information, 

♦/ 


DECPROG-HKCld’ 


Uaxmj neco(<o 


OCL 01 

RSCIN ALIENED 

9 

OS 

PARTI, 



10 

!>ATA_TYPP 

BIT 

(U’) 

10 

ty?f”si 

BIT 

(12) 

10 

T T.R kT 

BIT 

(«) , 

10 

"AHKER1 

BIT 

(<*) , 

10 

Cl 

BIT 

(20) 

10 

nsioN 

BIT 

(1) # 

10 

E 1 

BIT 

( /) , 

10 

RARKER2 

PIT 

(**) » 

10 

XNARKER 

BIT 

(H, , 

10 

TY»'F_S2 

BIT 

(12) 

10 

TT.RF2 

BIT 

(«) , 

10 

RAHKKUJ 

BIT 

CO , 

10 

C > 

BIT 

(20) 

10 

USIGN2 

BIT 

(1) » 

10 

K2 

BIT 

(0 » 

10 

CD 

BIT 

(H) , 

10 

CR 

BIT 

(0) , 

10 

NS 

BIT 

(0) , 

10 

K 

BIT 

(«) , 

10 

n 'IMS ED 

BIT 

(«) , 

10 

RN 

BIT 

(12) 

10 

NN 

BIT 

(IP) 

10 

Y 

BIT 

{^) , 

‘io 

M 

BIT 

(0) , 

10 

n 

BIT 

(0) , 

10 

I 

BIT 

(0) , 

10 

RXARKER 

BIT 

(«) , 

10 

C, 

PIT 

(1), 

10 

0 

BIT 

(1(>) 

10 

B 

BIT 

(IS) 

10 


BIT 

(If') 

10 

A • 

BIT 

(32) 

10 

llNIISE^2 

BIT 

( I'O 

10 

N 

BIT 

(16) 

10 

UNUSEO 1 

BIT 

(32) 

OS 

DATA (‘j.loS) n IT 

(16) 

OS 

XX 

BIT 

(16) 


/• STPnCTURE FOR INPUT FFCORDS ♦/ 
/♦ RFCFTVKS BINARY DIOITS FOR •/ 
/* CONVERSION TO DBCIRAL •/ 
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'Jsing Decode 


OCL 01 

TEJIORRC 

UNA LIGNSn, 



/* 

STRUCTURF TO 

RECEIVE 

♦/ 

OS 

PARTI, 





/* 

NUMBERS AFTER 

CONVERSION 

♦/ 

10 

DATA TYPE FIXED 

niNARY 

(31) 






10 

TYPE Si 

FIXED 

BINARY 

(31) 






10 

TIME? 

FIXED 

BINARY 

(31) 






10 

K ARK Em 

FIXED 

BINARY 

(31) 






10 

Cl 

FIXED 

BINARY 

(31) 






10 

nsiviN 

FIXED 

BINARY 

(31) 






10 

El 

FIXED 

BINARY 

(15) 






10 

r.AREER2 

FIXED 

BINARY 

(31) 






10 

X'lAPKER 

FIXED 

BINARY 

(31) 






10 

TYPE S2 

FIXED 

BINARY 

(31) 






10 

TT'1E2 

FIXED 

BINARY 

(31) 






10 

riARKER3 

FIXED 

BINARY 

(31) 






10 

C2 

FIXED 

BINARY 

(31) 



, 



10 

D SICK 2 

FIXED 

BINARY 

(31) 






10 

R2 

FIXED 

BINARY 

(15) 






10 

CD 

FIXED 

BINARY 

(31) 






10 

C'l 

F IXEP 

BINARY 

(31) 






10 

KS 

FIXED 

BINARY 

(-31) 






10 

K 

FIXED 

BINARY 

(31) 






10 

UMIISED 

FIXED 

BINARY 

(31) 






10 

RN 

FIXED 

BINARY 

(31) 






10 

N» 

F IXED 

BINARY 

(-31) 






10 

Y 

F TXED 

BINARY 

(31) 






10 

M 

FIXED 

BINARY 

(3 1) 






10 

D 

FIXED 

BINARY 

(31) 






10 

T 

FIXED 

BINARY 

(31) 






10 

R1ARKER 

FIXED 

BINARY 

(31) 






10 

G 

FIXED 

BINARY 

(31) 






10 

J 

FIXED 

BINARY 

(31) 






10 

D 

FIXED 

BINARY 

'31, 






10 

7. 

FIXED 

BINARY 

(15) 






10 

A 

FLOAT, 








10 

IINIISS?)2 

FIXED 

BINARY 

(31) 






10 

W 

FIXED 

BINARY 

(31) 






10 

MNIISEDO 

FIXED 

UTVARY 

(31) 






05 

DATA (13S7) FIXED 

BINARY 

(15) 






05 

XX 

F IXED 

BINARY 

(31) 

• 
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Using Deco'le 


DCL 0 1 OtlTPEC tlNALir.MFn, 
0*j PARTI, 

10 DATA_TYPE CHAR (14) 
10 Cl ~ CHAR(14) 

10 K CHAR (14), 

10 Y CHAP (14) 

10 H CHAR (14) 

10 D CHAP (14) 

10 I CHAP(14), 

05 DATA ( 133 7) CHAR (9) ; 


,/♦ STRUCTURE RETURNED BY 
/♦ SORT 


01 

W RTTSREC 

ilNALIGNED, 

/♦ 

STRUCTURE TO WRITE CUT 

*/ 


05 

PARTI, 


/♦ 

TO TAPE AFTER SORT 

♦/ 


10 

DATA TYPE 

CHAP (14; , 

/♦ 

WITH DATA CONSOLIDATED 

*/ 


10 

C 1 

C HAR (14) , 

/♦ 

FRO*1 2 RECORDS INTO 1 

*/ 


10 

Y 

CHAR ( 14) , 





10 

!1 

CHAR (14) , 

/♦ 

ONE OUTPUT RECORD IS 

*/ 


10 

D 

CHAP (14) , 

/* 

NOW TWO TIME PERIODS 

*/ 

/' 

10 

I 

CHAR (14) , 

/* 

(OWE DAY) OF DATA 

•/ 

1 

05 

DATA (1397) CHAR(9) , 





.05 DAT2(1397) CHAP (9) ; 


DECPROG-UKC332 


r 


to c cr 
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Using Decole 


DCL 

SYSUT1 

FILE 

RECORD 

INPUT, 

/• 

INPUT TAPE PILE 

*/ 


SYSUT2 

FILE 

RECORD 

OUTPUT, 

/♦ 

SORTED OUTPUT FILE 

*/ 


SYSPRINT 

FILE 

STHEAN 

OUTPUT; 

/• 

PRINT FILE 

♦/ 

DCL 

(EOF, DA DO) 

DIT (1) , 



/♦ 

FLAGS 

*/ 


(II, 

LI- IT, 

RKC NU*5DER) FIXED DIN (15) INIT(1), 

3LKS FIXED DIN (31), /♦ USED FOR COUNTERS ♦/ 

RETCODF FIXED 3IN(31,0); /* RETURN CODE FROH PLISRT*/ 


DCL (UNSPEC, 

PLI RRTC, 

PLISRTD, 

STRING, 

FLOOR) DIIILTIN; /♦ PLI FUNCTIONS ♦/ 


ON OVERFLOW DEGIN; 

PUT SKIP (4) EDIT ('***♦♦* OVERFLOW ERROR ENCOUNTERED ♦*♦♦*♦») 

(A) ; 

PUT SKIP (2) FDITC ♦♦♦♦♦* RETURN OUTPUT TO DOD CROVO ♦♦♦♦*♦») 

(A) ! 

PUT SKIP (2) EDITC****** FOR DSDUGGTNG ♦•♦♦♦♦!) 

(A) ; 

STOP; 

END; 


OF = ‘O’D; 

N F»'I)FI LE (SYSUT1) EOF = 'I'D; 

EXD FILE(SYSUTI) INTO(RECIN); /♦ READ 1ST RECORD (HEADER LADEL) INTO 

STRUCTURE RECIN TO THROW AWAY ♦/ 


DFCOROG-UKC.132 
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Using Decode 

Dn_LOOP: /• NU.'IBER 0? TTHE PERIODS TO TP\MSLATE*/ 

on DLKS»1 TO 100000 UMTI L (EOF) ;/* FOR ACTUAL PRODUCTION RUNS ♦/ 

/* DO 3LKS = 1 TO 2 ♦/ /• FOR DEBUGGING •/ 

READ_LOnP: 

/* 

Read roraining recorlij into structure RECIN and do conversion fro# 
binary to decimal by making the assignment from the binary structure 
(RECIN) to the fixed binary (31) structure (TENPREC) by name. 

Finally, convert structure to character representation for output 
to taoe in character structure OUTREC. 

♦/ 

/* PLISRTD calls GETREC 116 times before starting to sort records. 

After sort, PUTREC Is used to write to tape. Fields to bo sorted 
are YEAR, flONTH, DAY, TINE PERIOD, Q, C Descending, and K (grid), 
also Descenling. 

♦/ 

CALL PLTSRTD( 

' SORT FIELDS= (Ud.Su, CH, A, 1, 14,CH, A, 15, 14,CH, D, 29, 14,CH,D) •, 

• RECORD TYPE=F,LENC.TH=(12591) ', 

200000, RETCODC, GETREC, PUTREC, • ', 'AP') ; 

IF RETCODE-»=0 THEN DO; PUT D AT A (R E7CODE) ; /♦ BAD SORT V 

CALL PLIRFTC (RETCODE) ; 

STOP; 

END; 


/• <<<<<<< THE FOLLOWING STNT IS FOR DEBUGGING ONLY >»»» */ 

PUT SKIP EDI.T(» HAVING RETURNED FROM PLISRTD WITH RETCODE = •, RETCODE, 
• RETURN TO TOP OP LOOP AFTER READING DUMMY RECORD* ) (A, F ( 12, 3) , A) ; 
/♦ <<<<<<< REMOVE THE PRECFCING STMT FOR ACTUAL RUNS>>»»> •/ 


READ FILE(SYSUTI) INTO(RSCIN); 
END; /* 0’=’ DIG_LODP ♦/ 


/♦ READ 1ST RECORD (HEADER LABEL) INTO 
STRUCTURE RECIN TO THROW AWAY •/ 


PUT SKIP (9) SDIT(* 
PUT SKIP (9) SDIT(* 
Sl'^P; 


FORMAL FPOGRAM TERMINATION 
•,ULKS, • RECORDS PROCESSED 


*•) 

(CCL(IO) ,A) ; 
(COL (10) ,A) ; 


DECPROG-UKC332 
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Usinq Decode 

CIETPSC: /* GSTRSC REIiDS 116 RECORDS BEFORE ♦/ 

/♦ CALLING PLlRETC(fl) FOR PLTSRTD ♦/ 
PROCEDURE RETURNS (CH AR ( 1259 1) ) ; 
nCL K.'IT FIXED DIN(15) STATIC INIT (1); 

READ: IF KHT = 117 THEN DO; KNT * 1; 

RFC NUMBER = 1; 

CALL PLIRETC (9) ;/*lAST RECORD TO SORT*/ 
RETURK (STRING(OUTREC) ) ; 

END; 

BADQ = CONVERT; 

KNT = KNT 1; 

IF DADO THEN GOTO READ; /♦ SKIP THIS RECORD */ 

CALL PLTPETC(12); /* NEED MORE RECORDS FOR SORT •/ 

RRTURN(STRTNG(OUTREC) ) ; 

END; /* OP GFTREC ♦/ 


PUTREC: 

PROCEDURE (STRUCT) ; /♦ STRUCT DEF ON OUTRSC MAY BE BETTER */ 

DCL STRUCT CHAR (12591); 

DCL Q FIXED BIN (31); /♦ USED AS TEMPORARY VARIABLE */ 

DCL NEWCALL 3IT(1) STATIC INIT('I'B); 

STRING (OUTREC) = STRUCT; 

0 = OUTREC. CArA_TYPE.; 

“■ Q - ua •> (0*1) > 49 • (Q=2) ♦ 16 * (0*3) ♦ (0=4) ♦ 88 • (Q=5) ; 
(NOSTR2) :OUTREC. DATA_TYPS = Q; 

IF NEHCALL THEN DO; NEWCALL = 'O'B; 

R2C NUMBER = PEC_NUM3ER ♦ 1; 

HRITERBC = OUTREC, BY NAME; 

PND; 

ELSE 00; NEWCALL * M'B; 

WRITEPEC. DAT2 » OUTREC. DATA; 

WRITE FTLE(SYSUT2) FRCM (WRITEPEC) ; 

END; 

CALL PRINT; /♦ FOR DEBUGGING 0NLY--BE5T REMOVE ON PUN •/ 

CALI. ?LIRETC(4) ; 

END; /* OF PUTREC ♦/ 
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tlsing Decod# 

cnvvF:*??: /♦ convert tape data to decinal nosbers •/ 

PROCEDURE RET0RNS(0IT(1)) ; /♦ 1 FOR FAILHRE, 0 FOR SUCCESS •/ 

DCL 0 FIXED DIM (31); /♦ Q IS USED AS TEMPORARY VAR ♦/ 


READ FILE(SYS'JTl) INTO{RECIH); /* READ RECORD IHTO 

STRUCTURE RECIN •/ 

(N0STR2) : 

TSNPR3C.PART1*RECIM. PARTI, BY NAME;/* CONVERT FROM BINARY TO */ 

/* DECIMAL BY ASSIGNMENT •/ 

IF (TF.1PPEC.TYPE_S1 9) THEN 

RETURN (• 1 *0) /• TEST FOR DESIRED S ♦/ 

Q ~ TEMPREC.DATA_TY?E; 

Q = (Q=«8) ♦ 2 *”(Q-M9) ♦ 3 * (0=16) ♦ U • (Q=l) ♦ S * (Q=8B) ; 

IF (0 * 0) T'lCN RETURN (• 1' B) ; /* TEST FOR DESIRED O'S */ 

TEMPRSC.DATA^TYPE = 0; 

TEMPPEC.XX = RECIN. XX; 

UNSPEC (TEMPRKC. Z) = RECIN. Z; /* TRANSLATE DATA INTO */ 

UNSPECdEMPREC. A) = RECIN. A; /♦ THE CORRECT FORM. ♦/ 

IF (RECIN. B3IGN = ' 1 • B) THEN /* GIVE EXPONENT El CORRECT SIGN ♦/ 

TEMPREC.n * -(TEMPREC.E1) ; 

IF (RECIN. BSIGN2 M • B) THEN /* GIVE EXPONENT E2 CORRECT SIGN ♦/ 

TEMPREC.E2 = - (TEMPREC. E2) ; 

7SMPREC.C1 = FLOOR ( (TEMPREC, Cl * (10 TEMPREC. El)) ♦ .5); 

KK = 0; 

no JJ = 1 TO 37 BY 2; /♦ TRANSLATE ONLY ALTERNATE ARRAY */ 

HJJ = (JJ - 1 ) ♦ 145; 

DO II = 1. TO 14 5 BY 2; /* VALUES INTO DECIMAL VALDES */ 

KK = KK V 1; 

UNSPEC (TEMPREC. DATA (KK)) * R ECIN . DAT A (JJ J *■ II); 

TEMPREC. DATA (KK) * 

TEMPREC. A ♦ (TEMPREC. DATA (KK) ♦ (2** (TEMPREC. N-15) )) ; 

END; 

FED; 

(N0STP2) : 

OUTREC = TEMPREC, BY NAMS;/*CONVEPT TO CHARACTER FOR EBCDIC OUTPUT*/ 
l;'‘TUPN('0'D) ; 

END; /• OF CONVERT */ 



DECPROG-UKC3.32 


108 - 


PRINT: 

PROCEDURE; 


Using Decode 

/♦ OUTPUT CONVERTED VALUES •/ 


OUTPUT THE RESULT OP THE CONVERSION 
NECESSARY VARIA3LES AND THE VALUES 


, GIVING THE 
ASSOCIATED HITH 


THEN. 


*/ 

PUT 

PUT 


PAGE; 

SKIP EDIT 


PUT SKIP EDIT 


PUT SKIP EDIT 


PUT 

PUT 

/* 


SKIP (4) 
SKIP(2) ; 


(•RECORD NUH3ER — •) 
(RFC NtlNHER) 

CQ =• ,OUTRBC. DATA^TYPE) 
( *C1 *• ,0UTREC.C1)“ 

(•K =',OUTREC.K) 

*• ,OUTREC. Y) 
*',OUrREC.N) 

*• ,0tlTR2C.D) 

•*• ,OUTREC. I) 
> 


(•Y 
( *N 
( '0 
(•I 

EDITC 


PA'iA 


(COL(l) ,A) 

(P(7)) ; 

(A,F(6)) 

(X(24) ,A,F(C)) 
(X(24) ,A,F(6)) ; 
(^,F(6)) 

(X(24) ,A,F(6)) 
(X(24) ,A,F(6)) 

(X (24) ,A,P(6)) i 
FROM RECORD <• 


(COL (30) ,A) 


Output data entries along with the 
refjrence. The following LIMIT may 
A miltiple of 5 is preferred. 

V 

LIMIT * 225; 

no KK » 1 TO LIMIT BY 5; 

PUT SKIP EDIT (•(',KX,' - *,KK*4,*) 
OUTREC.UATA (KK) ,OUTREC.D 
OUTREC.PATA (KK»2) ,OUTREC 
(COL(l) ,A,F(5) ,A,F(5) ,A, 


array subscripts, for 
be from 1 to 1387. 


ATA (KK^I) , 

. DATA (KK^3) ,0UTPEC. DATA (KK»4) ) 
5 (A,X(2))); 


END; 


FND; 




RETURN; 




END; /* 

OF 

• 

PRINT 

•/ 

; /• 

0? 

DECODE 

•/ 
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